本章主要采用静水有限元模拟方法,结合成像测井分析结论和本区构造特点,对现代地应力场进行分析模拟。有限元方法最核心的思想是将储层岩石分割为一个个离散的单元,然后在各个单元内进行分析,构建平衡方程,最后将各个单元拼接在一起,得到模拟结果。目前有限元模拟方法主要基于力学变分原理,采用加权方法对基本力学方程进行离散处理。根据上述静水有限元模拟基本原理,需要构建8点有限元网格的型函数,本次模拟令......
2023-08-20
根据裂缝的成像测井识别方法,收集了研究区及附近区域的30口井的成像测井资料(图13),进行了裂缝识别分析,共识别出裂缝800余条天然裂缝(表2)。
图13 裂缝识别
表2 研究区及附近区域的裂缝产状表(部分)
续表
续表
续表
根据本区成像测井解释结果(表2),研究层段裂缝共可以分为天然裂缝和诱导裂缝两种类型。诱导裂缝包括井壁崩落(BO)和水力压裂缝(DIT)两种。BO缝的走向多为南北向,DIT的走向多为近东西向。在本文第3章中将对BO缝和DIT缝进行详细论述。本区天然裂缝可以分为两组,一组近南北向裂缝,另一组为近东西向裂缝,且以近南北向裂缝为主。根据裂缝分组的结果,本研究区一共存在2组裂缝(图14):一组为北北东方向,另一组为北东东方向。在每一个裂缝组内,需要确定该裂缝组的裂缝密度和间距,从而为本文第5章的裂缝多相渗流模拟做好数据准备。
图14 裂缝分组图
在确认了裂缝分组结果之后,裂缝的走向和倾角属性已经在分组结果中获得体现。另一个重要的裂缝属性是裂缝密度属性。
裂缝密度特征是后期裂缝建模和多相渗流耦合的重要指标。一般裂缝密度指单位长度内与某一方向垂直的裂缝条数:
其中M为裂缝条数,L为测量长度。
裂缝密度是裂缝建模的重要依据。单井裂缝密度可以通过成像测井资料解释和岩心观察获得,但是对于井间的裂缝,则没有直接的地质资料提供支持,所以井间裂缝的产状分析是裂缝建模的一个难题。单井裂缝密度仍然是开发裂缝建模技术的重要基础资料。除了裂缝密度,裂缝密度的倒数,即裂缝间距D也是影响裂缝储层中多相渗流的重要参数:
由于裂缝密度特指某一方向上的单位长度的裂缝数量,所以需要在裂缝的两个分组内分别计算裂缝密度(表3)。在建立各个分组的裂缝密度之后,还需要确认裂缝孔隙度fφ:
其中h为裂缝开度,D为根据裂缝密度计算的裂缝间距。
由于裂缝开度h一般要远远小于裂缝间距D,所以裂缝孔隙度计算公式可以简写为:
表3 裂缝密度表(部分)
续表
(组I 为近东西向裂缝,组II为近南北向裂缝。有效裂缝均为成像测井可探测的宏观裂缝)
由本次研究所识别出的裂缝产状分布数据可知,研究区块第一组东西向的裂缝密度为0.04 m-1,第二组南北向的裂缝密度为0.18 m-1。本次研究受成像测井分辨率限制,所识别裂缝的平均开度为800μm左右。根据裂缝孔隙度计算公式,第一组裂缝的平均孔隙约为0.003%,第二组裂缝的平均孔隙度约为0.016%。根据裂缝各个组系的裂缝密度和裂缝孔隙度以及平均开度,就可以开展计算各个组系的裂缝渗透率。裂缝渗透率的计算将在第5章多相渗流分析中给予详细论述。
有关储层裂缝建模与多相流耦合理论研究的文章
本章主要采用静水有限元模拟方法,结合成像测井分析结论和本区构造特点,对现代地应力场进行分析模拟。有限元方法最核心的思想是将储层岩石分割为一个个离散的单元,然后在各个单元内进行分析,构建平衡方程,最后将各个单元拼接在一起,得到模拟结果。目前有限元模拟方法主要基于力学变分原理,采用加权方法对基本力学方程进行离散处理。根据上述静水有限元模拟基本原理,需要构建8点有限元网格的型函数,本次模拟令......
2023-08-20
根据本文第4章的研究成果,SEFE裂缝建模方法最终所输出的裂缝建模结果包括一系列应力场分布模型、一个裂缝密度分布模型和一个裂缝离散片元模型。本次研究所建立的裂缝属性模型,全部基于SEFE裂缝建模方法的输出结果完成。为了对比SEFE裂缝建模的效果,将断层控制的DFN裂缝建模结果和SEFE裂缝建模结果进行了对比,如表10。可见SEFE裂缝建模结果在断层间断区域的裂缝建模效果较DFN更好,与井点裂缝分布数据更为吻合。......
2023-08-20
恢复后的层面模型,即被认定为一个可接受的变形前模型,而构造正演只需对恢复过程进行回放即可获得。目前,有些研究者开始认为裂缝的分形维特征的本质在于介质本身的属性特征,如材料刚度。在裂缝建模的基础上,研究人员提出了多尺度裂缝建模与多相流耦合分析模型。这种方法在进行裂缝性介质内的多相流模拟时,将裂缝与孔隙视为统一的均匀介质进行处理。......
2023-08-20
由以上方程可知,整个双孔双渗模型基本控制方程组共包含S wm,S om,S wf,S of,Pw m,Pom,Pw f,Pof 共8个未知量,方程组的数目也是8个。Pcowm,Pcowf,Krwm,Krom,则由两相流实验得到。......
2023-08-20
无网格扩展有限元方法的工作原理是允许网格在不进行重新劈分的情况下模拟网格单元中的不连续面。无网格扩展有限元方法就是基于这一点,在传统有限元方法的基础上引入非连续性函数表示裂缝处的非连续位移。在非各向同性的储层岩石介质中,裂缝尖端奇性场依赖于裂缝尖端的具体位置,因此要准确地模拟裂缝尖端奇性场,需要不断追踪裂缝的扩展轨迹,降低模拟速度。一种替代方案是基于裂缝扩展沿断裂径迹扩展行为的扩展有限元模拟方法。......
2023-08-20
由以上方程知,双孔单渗模型基本控制方程组有S wm,S om,S wf,S of,Pw m,Pom ,Pw f,Pof共8个未知量,方程组的数目也是8个。Pcowm,Pcowf,Krwm,Krom,K rwf和K rof则由两相流实验得到。......
2023-08-20
图69 SEFE裂缝模拟结果图70 裂缝开度模拟结果图71 裂缝滑移程度模拟结果利用多场克里金模拟系统,综合成像测井裂缝解释资料和裂缝扩展正反演系统的模拟结果,最终建立了大情字井地区上白垩统青山口组的裂缝模型。模拟结果表明,单井点裂缝属性预测误差MAEcv=6.4%。......
2023-08-20
对于第二种裂缝系统反演方法,目前较为常用的方法有梯度法、牛顿迭代法、蒙特卡洛模拟法等。裂缝反演系统,就转化为求解目标函数的极值问题。因此搜索步距应当和目标函数在最新搜索点处的梯度成反比。确定了搜索方向,但还面临着裂缝反演问题中的一个重要障碍。这里,本次研究在共轭向量搜索的基础上,加入模拟退火......
2023-08-20
相关推荐