CN111814364B 一种岩溶储层演化数值模拟方法 （中国石油化工股份有限公司）

揭示碳酸盐岩缝洞型储层的形成机制提供理论2步骤3.3、采用有限体积和有限元的方法对多类型储层流-固-化耦合数学模型进行离3sf为质量源/汇项。up和RC分别为所述力学平衡方程、所述流体质量守恒方程和所述带有化学反度。7.根据权利要求6所述的岩溶储层演化数值模拟方4KarstAquiferEvolutioninFracturedRocks的论文，建立了溶洞形成过程的管流模型。2009年KaufmannG.又发表了题为ModellingKarstGeomorphologyonDifferent[0005]2005年PangaM.等在AIChEJournal上发表题为Two-scalecontinuummodel发表题为ModelingandAnalysisofWormholeFormationinReactiveDissolution5[0017]在一个实施方式中，所述步骤1.1[0024]D·(De-bp⃞)+png=0(3)6ρb分别为流体密度和岩石骨架与孔隙内流体体积平均密度，为达西f和Cs分别为孔隙内的钙离子浓度和固-液交界面上的钙离子浓度，De为弥散系数张[0031]步骤3.3、采用有限体积方法和有限元的[0036]在一个实施方式中，其中步骤3.3中对所述数值离散得到的式(19)采用Newton-7[0046]图1是本发明实施例中的基于嵌入式离散裂缝的流-固-化耦合岩溶演化数值模拟[0051]图6是本发明实施例中的岩溶演化构造裂缝模型模拟得到了100年后储层孔隙度[0052]图7是本发明实施例中的岩溶演化构造裂缝模型模拟得到了200年后储层孔隙度[0053]图8是本发明实施例中的岩溶演化构造裂缝模型模拟得到了300年后储层孔隙度[0056]如图1所式为本发明是实例中的岩溶演化数值模拟方法流程示意图，下面结合对本发明中的岩溶演化数值模拟方法做一个具体[0057]本发明提出的是一种基于嵌入式离散裂缝的流-固-化全耦合的岩溶演化数值模式离散裂缝方法(EDFM)建立裂缝与基质耦合地质模型的优势可概括为：1)裂缝建模方便，8行网格剖分(其中顶点可不在模拟区域范围内，非平面天然裂缝可以通过多个平面裂缝拼[0066]在一个优选的实施例中，步骤1.1[0070]下面对系统的控制方程进行具体的介绍，所述流-固-化耦合系统的控制方程包9f和Cs分别为孔隙内的钙离子浓度和固-液交界面上的钙离子浓度，De为弥散隙结构，可通过建立孔隙结构与物理化学参数之间的关系(Carman-Kozeny方程)来考虑物[0088]在步骤2中得到多类型储层流-固-化耦合数学模型以后，需要对该数学模型进行[0091]步骤3.3、采用有限体积和有限元的方式(11)所示的六面体单元有限元插值函[0131]根据岩溶演化数学模型特点，设计相应的C++程序设计框架，按照模块化设计理展性强的设计目标。如图4所示，整个岩溶演化模拟系统KarstEvolSys按功能分为三个部变形和弹塑性变形的本构模型模块ConstituModel、根据网格自动生成两点或多点流动格[0141]使用本发明提出的方法模拟含天然构造裂缝的储层演化过程，主要的模拟流程流体压力和温度分布分别按照静水压力和地温梯度(按照25℃/km)分别计算得到，二氧化