[论文解读] Existence and Regularity of Optimal Solution for a Dead Oil Isotherm Problem

Moulay Rchid, Delfim F. M. Torres|ArXiv.org|Aug 15, 2006

Geological Studies and Exploration参考文献 7被引用 4

一句话总结

本文通过最优控制理论与Sobolev嵌入技术，建立了描述多孔介质中死油等温流的非线性PDE系统最优解的存在性与正则性。关键贡献在于，在数据与系数满足适当的结构性和相容性条件时，证明了最优控制的存在性及其充分的Sobolev正则性。

ABSTRACT

We study a system of nonlinear partial differential equations resulting from the traditional modelling of oil engineering within the framework of the mechanics of a continuous medium. Existence and regularity of the optimal solutions for this system is established.

研究动机与目标

• 建立用于描述多孔介质中不可压缩两相流的非线性PDE系统的最优解存在性，该问题与烃类采收相关。
• 在非线性项与边界条件的结构性假设下，分析最优解在Sobolev空间中的正则性。
• 通过在控制中施加足够正则性，并确保初值与边界数据之间的相容性，保证控制到状态算子的Gâteaux可微性。
• 制定一个代价泛函，以最小化油饱和度与压力偏离期望分布的偏差，同时对控制努力及其时间导数施加惩罚。
• 为增强油藏采收中的最优井位布置与注水策略提供严格的数学基础。

提出的方法

• 建立一个双方程PDE系统：一个关于油饱和度 $ u $ 的退化抛物方程，以及一个关于压力 $ p $ 的抛物方程，其中包含非线性项 $ abla \to \nabla p $、$ \nabla \to \nabla u $，以及系数 $ \theta(u), g(u), d(u) $。
• 引入代价泛函 $ J(u,p,f) $，结合 $ u $ 与 $ p $ 的 $ L^2 $ 误差项，以及对 $ f $ 与 $ \nabla_t f $ 的 $ L^{2q_0} $ 惩罚项，其中 $ 1 < q_0 < 2 $，以确保控制的正则性。
• 应用 De Giorgi–Nash–Ladyzhenskaya 定理与嵌入定理，推导出时间导数的更高可积性与Hölder连续性。
• 对时间导数 $ \nabla_t u $、$ \nabla_t p $，以及 $ \nabla_t \nabla u $、$ \nabla_t \nabla p $ 使用能量估计与Gronwall不等式，以控制其 $ L^4 $-范数。
• 通过从 $ L^{2q_0} $-右端项出发并结合椭圆正则性，利用逐次逼近法建立 $ u \to W_{2q_0}^{2,1}(Q_T) $ 与 $ p \to W_{2q_0}^{2,1}(Q_T) $。
• 采用紧致性论证与Sobolev空间中的弱收敛，处理极小化序列的极限过程，从而证明最优解的存在性。

实验结果

研究问题

• RQ1 在给定的代价泛函与约束条件下，死油等温模型是否存在最优控制？
• RQ2 最优状态与控制变量可保证哪些正则性性质（Sobolev、Hölder、$ L^4 $-可积性）？
• RQ3 系数 $ \theta(u), g(u), d(u) $ 及其导数如何影响解的正则性？
• RQ4 在所作假设下，能否建立控制到状态算子的Gâteaux可微性？
• RQ5 时间导数惩罚项 $ \beta_2 \frac{1}{2} \norm{\nabla_t f}_{L^2(Q_T)}^2 $ 在确保控制正则性方面起什么作用？

主要发现

• 在 $ \theta, g, d $ 及初值/边界数据满足所述假设的条件下，死油等温问题在给定函数空间设定下存在最优解 $ (u^*, p^*, f^*) $。
• 时间导数 $ \nabla_t u $ 属于 $ L^4(Q_T) $，且 $ u \to W_{2q_0}^{2,1}(Q_T) $，其中 $ 1 < q_0 < 2 $，表明其正则性足以保证Gâteaux可微性。
• 压力状态 $ p $ 同样属于 $ W_{2q_0}^{2,1}(Q_T) $，这是由于其PDE右端项属于 $ L^{2q_0}(Q_T) $，且适用相同的估计。
• 正则性估计 $ \norm{\nabla u}_{L^{2q}(Q_T)} \to c(\norm{f}_{L^{2q}(Q_T)} + \norm{u_0}_{W^{1,2q}(\bar{\theta})}) $ 成立，其推导基于引理2.1。
• 压力的时间导数满足 $ \norm{\nabla_t p}_{L^2(0,T;L^2(\theta))} \to c $，且 $ \nabla_t p \to L^2(0,T;L^2(\theta)) $，其界依赖于 $ \norm{\nabla p}_{L^4(\theta)} $。
• 解满足能量型不等式 $ \frac{d}{dt} \norm{\nabla_t u}_{L^2}^2 + \norm{\nabla_t \nabla u}_{L^2}^2 \to c \norm{\nabla p}_{L^4}^4 \norm{\nabla_t u}_{L^2}^2 $，通过Gronwall引理可导出指数有界性。