×

流与地质力学耦合序列方法的稳定性和收敛性:排水和不排水劈裂。 (英语) 兹比尔1228.74106

摘要:我们对流动与地质力学耦合的序贯方法进行了稳定性和收敛性分析,其中首先解决了力学子问题。我们考虑缓慢变形,因此惯性可以忽略不计,力学问题由椭圆方程控制。我们利用Biot的自洽理论得到了经典抛物线型流动问题。我们使用了一个广义中点规则(0和1之间的参数α)时间离散化,并考虑了两种经典的序列方法:排水和不排水分裂。
冯·诺依曼方法为线性多孔弹性问题提供了精确的稳定性估计。带后向欧拉时间离散化的排水劈裂是条件稳定的,其稳定性仅取决于耦合强度,与时间步长无关。使用中点规则((alpha=0.5))的排水劈裂是无条件不稳定的。混合时间离散化,力学为(alpha=1.0\),流动为(alfa=0.5\),与反向Euler格式具有相同的稳定性。von Neumann方法表明,当(α\geq 0.5)时,不排水劈裂无条件稳定。我们通过能量法将稳定性分析扩展到非线性领域(多孔弹塑性)。众所周知,排水劈裂不继承连续介质问题的收缩性,从而排除了无条件稳定性。对于不排水劈裂,当\(alpha\geq 0.5\)时,我们表明它是B-稳定的(因此在算法水平上是无条件稳定的)。
我们还分析了排水分裂和不排水分裂的收敛性,并从矩阵代数和谱分析中导出了先验误差估计。我们证明了具有固定迭代次数的排水分裂即使在稳定的情况下也不收敛。对于可压缩系统(即有限Biot模量),具有固定迭代次数的不排水分裂是收敛的。对于近不可压缩系统(即非常大的Biot模量),不排水劈裂失去了一阶精度,并且在时间上变得不具有传染性。
我们还研究了在全迭代序列格式中使用这两个分裂时的收敛速度。当介质渗透率较高或时间步长较大(对应压力扩散较大)时,排水劈裂的误差放大率较低,因此收敛速度快于不排水劈缝。在低渗透率和小时间步长的情况下,情况正好相反。
我们提供了数值实验,支持线性和非线性情况下排水和不排水劈裂的所有稳定性和收敛性估计。我们还表明,我们的空间离散化(流动的有限体积和力学的有限元)消除了早期固结问题中记录良好的虚假不稳定性。

MSC公司:

74平方米 有限差分法在固体力学问题中的应用
74S05号 有限元方法在固体力学问题中的应用
76个M12 有限体积法在流体力学问题中的应用
74层10 流固相互作用(包括气动和水弹性、孔隙度等)
76秒05 多孔介质中的流动;过滤;渗流
65个M12 含偏微分方程初值和初边值问题数值方法的稳定性和收敛性
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Truesdell,C。;Toupin,R.A.,《经典场论》(Flügge,S.,《经典力学和场论原理》,Handbuch der Physik,第III/1卷(1960年),Springer-Verlag:Springer-Verlag Berlin),226-793
[2] 瓦夫,R。;Evans,E.A.,红细胞膜的热弹性,生物物理。J.,26,115-131(1979年)
[3] Kaczmarek,M。;Subramaniam,R.P。;Neff,S.R.,《脑积水的流体力学:圆柱几何的稳态解》,布尔。数学。《生物学》,59,2,295-323(1997)·Zbl 0904.92026号
[4] Biot,M.A.,《三维固结的一般理论》,J.Appl。物理。,12, 155-164 (1941) ·肯尼迪67.0837.01
[5] Coussy,O.,《多孔介质力学》(1995),约翰·威利和儿子:约翰·威利与儿子奇切斯特,英国
[6] 王海峰,《线性多孔弹性理论》(2000),普林斯顿大学出版社
[7] 巴顿,C.A。;医学博士Zoback。;Moos,D.,晶体岩石中潜在活动断层的流体流动,《地质学》,23,8,683-686(1995)
[8] Wiprut,D。;Zoback,M.D.,《北海北部先前休眠的正断层沿线的断层复活和流体流动》,《地质学》,28,7,595-598(2000)
[9] 怀特,J.A。;Borja,R.I.,用于耦合固体变形/流体扩散的稳定低阶有限元及其在断层带瞬变中的应用,计算。方法应用。机械。工程,197,4353-4366(2008)·Zbl 1194.74480号
[10] Jain,A.K。;Juanes,R.,沉积物中气体入侵的优先模式:多相流体流动和沉积物力学耦合的颗粒尺度机制模型,地球物理杂志。Res.,114,B08101(2009)
[11] 霍尔兹曼,R。;Juanes,R.,《可变形无序介质中从指进到破裂的交叉》,Phys。版本E,82,046305(2010)
[12] Merle,H.A。;肯蒂,C.J.P。;van Opstal,G.H.C。;Schneider,G.M.G.,《巴查奎罗研究——压实驱动/溶解气体驱动油藏行为的综合分析》,J.Petrol Technol。,1107-1114(1976年)
[13] Kosloff,D。;斯科特·R·F。;Scranton,J.,《威明顿油田沉降的有限元模拟:I.线性建模,构造物理学》,65,339-368(1980)
[14] Bruno,M.S.,补贴引起的油井故障,SPE钻井工程,148-152(1992)
[15] Fredrich,J.T。;Arguello,J.G。;Deitrick,G.L。;de Rouffignac,E.P.,《贝尔里奇硅藻土油田SPE油藏储层压实、地面沉降和套管损坏的地质力学建模》。评估。工程师,第3、4、348-359页(2000年)
[16] Settari,A。;Mourits,F.M.,《储层与地质力学耦合模拟系统》,《石油工程学报》,3,3,219-226(1998)
[17] Settari,A。;Walters,D.A.,《地质力学与储层建模耦合进展及其在储层压实中的应用》,《石油工程学报》,6,3,334-342(2001)
[18] Maingue,M。;Longuemare,P.,《耦合流体流动和岩石力学:储层和地质力学模拟器之间的部分耦合公式》,《油气科学》。技术。,57, 355-367 (2002)
[19] Minkoff,S.E。;斯通,C.M。;布莱恩特,S。;Peszynska,M。;Wheeler,M.F.,《耦合流体流动和地质力学变形建模》,《石油科学杂志》。工程,38,37-56(2003)
[20] 托马斯·L·K。;Chin,L.Y。;皮尔逊·R·G。;Sylte,J.E.,《地质力学与储层模拟耦合》,《石油工程学会期刊》,8,4,350-358(2003)
[21] Tran,D。;Settari,A。;Nghiem,L.,储层模拟器和地质力学模块之间的新迭代耦合,Soc.Petrol Eng.J.,9,362-369(2004)
[22] Dean,R.H。;盖,X。;斯通,C.M。;Minkoff,S.E.,《耦合多孔流体和地质力学技术的比较》,《石油工程学报》,11,1,132-140(2006)
[23] Jha,B。;Juanes,R.,《模拟耦合流体流动和油藏地质力学的局部保守有限元框架》,《岩土学报》。,2, 3, 139-153 (2007)
[24] 刘易斯,R.W。;Sukirman,Y.,变形饱和油藏中三相流动的有限元建模,Int.J.Numer。分析。方法地质力学。,17, 8, 577-598 (1993) ·Zbl 0785.76042号
[25] 刘易斯,R.W。;Schrefler,B.A.,多孔介质静态和动态变形及固结的有限元方法(1998),John Wiley and Sons:John Willey and Sons Chichester,英国·Zbl 0935.74004号
[26] Armero,F.,完全饱和条件下有限应变下耦合多孔塑性乘法模型的公式和有限元实现,计算。方法应用。机械。工程,171,205-241(1999)·Zbl 0968.74062号
[27] Jean,L。;Maingue,M。;马森,R。;Vidal-Gilbert,S.,《加速耦合地质力学-储层模拟的收敛》,国际期刊数值。分析。方法地质力学。,31, 1163-1181 (2007) ·Zbl 1196.74285号
[28] Park,K.C.,《孔隙流体-土壤相互作用分析的分区溶液程序的稳定性》,国际期刊数值。方法。工程,第19、11、1669-1673页(1983年)·Zbl 0519.76095号
[29] 齐恩凯维奇,O.C。;保罗·D·K。;Chan,A.H.C.,《土-孔-流体相互作用问题的无条件稳定交错解程序》,《国际数值杂志》。方法。工程,26,5,1039-1055(1988)·Zbl 0634.73110号
[30] Armero,F。;Simo,J.C.,《耦合热力问题的一种新的无条件稳定分步方法》,国际J·数值。方法。工程,35,737-766(1992)·Zbl 0784.73085号
[31] Prevost,J.H.,耦合场问题同时积分的分区求解程序,Commun。数字。方法工程,13,239-247(1997)·Zbl 0878.73073号
[32] Schrefler,文学学士。;西蒙尼,L。;Turska,E.,热工水力力学问题中的标准交错和交错牛顿格式,计算。方法应用。机械。工程,14493-109(1997)·Zbl 0891.76073号
[33] 黄,M。;Zienkiewicz,O.C.,《土-孔-流体耦合动力学问题的新无条件稳定交错解程序》,国际数值杂志。方法。工程,431029-1052(1998)·Zbl 0945.74069号
[34] 马萨诸塞州费利帕。;Park,K.C.,耦合机械系统的交错瞬态分析程序:公式,计算。方法应用。机械。工程师,24,61-111(1980)·Zbl 0453.73091号
[35] Farhat,C。;帕克,K.C。;Dubois-Pelerin,Y.,耦合热弹性问题瞬态有限元分析的无条件稳定交错算法,计算。方法应用。机械。工程,85,349-365(1991)·Zbl 0764.73081号
[36] Armero,F。;Simo,J.C.,非线性耦合热塑性的先验稳定性估计和无条件稳定乘积公式算法,国际J塑性,9749-782(1993)·Zbl 0791.73026号
[37] Simo,J.C.,非线性热传导、塑性和粘塑性演化的时间离散变分问题的非线性稳定性,计算。方法应用。机械。工程,88,111-131(1991)·Zbl 0751.73066号
[38] Araüjo,A.,关于分裂方法的B-稳定性的注释,计算。目视检查。科学。,6, 53-57 (2004) ·Zbl 1067.65077号
[39] Hundsdorfer,W.H。;Spijer,M.L.,关于Runge-Kutta方法的B-稳定性的注释,Numer。数学。,36, 319-331 (1981) ·Zbl 0441.65065号
[40] Romero,I.,《非线性热力系统的热力学一致性时间步进算法》,国际J·数值。方法。工程,79,706-732(2009)·Zbl 1171.80305号
[41] Richtmyer,R.D。;Morton,K.W.,《初值问题的差分方法》(1967),John Wiley and Sons·Zbl 0155.47502号
[42] Brenan,K.E。;坎贝尔,S.L。;Petzold,L.R.,微分代数方程初值问题的数值解(1996),SIAM:SIAM Philadelphia·Zbl 0844.65058号
[43] Chorin,A。;休斯·T·R。;麦克拉肯,M.F。;Marsden,J.E.,《产品公式和数值算法》,Commun。纯应用程序。数学。,205-256(1978年)·Zbl 0358.65082号
[44] Lapidus,M.L.,Trotter-Lie公式的推广,积分。Equat公司。操作。理论。,4/3, 366-413 (1981) ·Zbl 0463.47024号
[45] Turska,E。;Wisniewski,K。;Schrefler,B.A.,瞬态问题交错求解程序的误差传播,计算。方法应用。机械。工程,114177-188(1994)
[46] Vijalapura,P.K。;Govindjee,S.,《高度非线性强耦合问题的自适应混合时间步长方案》,国际期刊Numer。方法。工程,64,819-848(2005)·Zbl 1122.74525号
[47] 维贾拉普拉,P.K。;应变J。;Govindjee,S.,指数-1微分代数方程的分步方法,J.Compute。物理。,203, 305-320 (2005) ·Zbl 1063.65074号
[48] Simo,J.C。;Govindjee,S.,《塑性和粘塑性的非线性B稳定性和对称保持回归映射算法》,国际期刊Numer。方法。工程师,31151-176(1991)·Zbl 0825.73959号
[49] Schrefler,B.A。;西蒙尼,L。;Turska,E.,热工水力力学问题中的标准交错和交错牛顿格式,计算。方法应用。机械。工程,14493-109(1997)·Zbl 0891.76073号
[50] 维米尔,P.A。;Verruijt,A.,《有限元固结的精度条件》,国际期刊数字。分析。方法地质力学。,5, 1-14 (1981) ·Zbl 0456.73060号
[51] ƀenišek,A.,Biot固结理论中的存在唯一性定理,Aplik。材料。,29, 194-210 (1984) ·Zbl 0557.35005号
[52] J.Wan,L.J.Durlowsky,T.J.R.Hughes,K.Aziz,耦合地质力学-油藏流动模拟的稳定有限元方法,摘自:SPE油藏模拟研讨会,德克萨斯州休斯顿(SPE 79694),2003年。;J.Wan,L.J.Durlowsky,T.J.R.Hughes,K.Aziz,耦合地质力学-油藏流动模拟的稳定有限元方法,摘自:SPE油藏模拟研讨会,德克萨斯州休斯顿(SPE 79694),2003年。
[53] Phillips,P.J。;Wheeler,M.F.,多孔弹性混合和连续Galerkin有限元方法的耦合I:连续时间情况,计算。地质科学。,11, 131-144 (2007) ·Zbl 1117.74015号
[54] 菲利普斯,P.J。;Wheeler,M.F.,多孔弹性的混合和连续Galerkin有限元方法的耦合II:离散时间情况,计算。地质科学。,11, 145-158 (2007) ·Zbl 1117.74016号
[55] Geertsma,J.,流体压力下降对多孔岩石体积变化的影响,Trans。AIME,210,331-340(1957)
[56] Borja,R.I.,《关于饱和和非饱和多孔连续介质中的机械能和有效应力》,国际固体结构杂志。,43, 6, 1764-1786 (2006) ·Zbl 1120.74446号
[57] 阿齐兹,K。;Settari,A.,《石油储层模拟》(1979年),爱思唯尔出版社:伦敦
[58] Hughes,T.J.R.,《有限元方法:线性静态和动态有限元分析》(1987),新泽西州普伦蒂斯·霍尔:普伦蒂塞·霍尔恩格尔伍德克利夫斯·Zbl 0634.73056号
[59] Kim,J。;Tchelepi,H.A。;Juanes,R.,《流动与地质力学耦合顺序方法的稳定性、准确性和效率》,Soc.Petrol Eng.J.(2011)·Zbl 1228.74106号
[60] Aavatsmark,I.,四边形网格多点通量近似介绍,计算。地质科学。,6, 405-432 (2002) ·兹比尔1094.76550
[61] 艾萨克森,E。;Keller,H.B.,《数值方法分析》(1966),John Wiley and Sons:John Willey and Sons New York·Zbl 0168.13101号
[62] Strikwerda,J.C.,《有限差分格式和偏微分方程》(2004),SIAM:SIAM Philadelphia,PA·Zbl 1071.65118号
[63] 米加,麻省理工学院。;Paulsen,K.D。;Kennedy,F.E.,Von Neumann对Biot一般二维固结理论的稳定性分析,Int.J.Numer。方法。工程,43,955-974(1998)·Zbl 0940.76035号
[64] Marsden,J.E。;Hughes,T.J.R.,《弹性的数学基础》(1983),《普伦蒂斯·霍尔:普伦蒂斯霍尔·恩格尔伍德·克利夫斯,新泽西州》(经更正重印,多佛,纽约,1994)·Zbl 0545.73031号
[65] Butcher,J.C.,隐式Runge-Kutta方法的稳定性,BIT,15,358-361(1975)·Zbl 0333.65031号
[66] Simo,J.C。;Hughes,T.J.R.,计算非弹性(1998),Springer:Springer纽约·Zbl 0934.74003号
[67] Simo,J.C。;Taylor,R.L.,速率相关弹塑性的一致切线算子,计算。方法应用。机械。工程,38,101-118(1985)·Zbl 0535.73025号
[68] Taylor,R.L。;Simo,J.C.,平面应力弹塑性的返回映射算法,国际期刊数值。方法。工程,22649-676(1986)·Zbl 0585.73059号
[69] Thomée,V.,抛物问题的Galerkin有限元方法(1997),Springer-Verlag:Springer-Verlag Berlin·Zbl 0884.65097号
[70] 图尔斯卡,E。;Schrefler,B.A.,关于固结问题分区求解程序的收敛条件,计算。方法应用。机械。工程,106,51-63(1993)·Zbl 0783.76064号
[71] Strang,G.,《线性代数及其应用》(1988),Brooks/Cole
[72] 博尔贾,R.I。;Lee,S.R.,Cam-clay塑性,第一部分:弹塑性本构关系的隐式积分,计算。方法应用。机械。工程,78,49-72(1990)·Zbl 0718.73034号
[73] Kim,J。;切莱皮,H.A。;Juanes,R.,《耦合流和地质力学序列方法的稳定性和收敛性:固定应力和固定应变分裂》,计算。方法应用。机械。工程,200,13-16,1591-1606(2011)·Zbl 1228.74101号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。