CN113919196B 储层三维应力场模拟方法、模拟系统、终端、存储介质 （中国石油大学(华东)）

本发明属于油气田开发地质学和石油工程在地震相约束下，利用增强曲率属性跟踪裂缝续面导致的局部强度弱化和应力集中对整个岩程序实现地质模型和地应力有限元模型的网格2在地震相约束下，利用增强曲率属性跟踪裂缝线，实现在基于建立的裂缝确定性模型，利用相关网格转换程序进行地质模型步骤五，通过从地表至目的层的全时窗叠前三维地震反演得到步骤六，把断层作为不连续体引入到三维岩石力学参数不同层位的地应力实验结果对测井计算地应力步骤十一，在确定地应力反演工区以及参演参数后，32.根据权利要求1所述的基于确定性多尺度裂缝建模和竞争粒子群算法的储层三维应通过收集研究区野外露头、单井岩心观察、录井岩屑及前期地质研究采用分组方式进行弥散裂缝的确定性建模对小尺度弥散裂缝采用蚁群跟踪算法进行拓扑参数，建立相关天然裂缝网络进行大尺度离散裂缝分析，i3.根据权利要求1所述的基于确定性多尺度裂缝建模和竞争粒子群算法的储层三维应s4.根据权利要求1所述的基于确定性多尺度裂缝建模和竞争粒子群算法的储层三维应4nsMPa裂缝的正应力与裂缝的刚度成正比，裂缝的尺寸则与裂缝的刚度成反比s5.根据权利要求1所述的基于确定性多尺度裂缝建模和竞争粒子群算法的储层三维应6.根据权利要求1所述的基于确定性多尺度裂缝建模和竞争粒子群算法的储层三维应所述步骤十具体包括：有限元网格采用随机搜索方法对扭曲、拉长的7.一种根据权利要求1～6任意一项所述基于确定性多尺度裂缝建模和竞争粒子群算法的储层三维应力场模拟方法的基于确定性多尺度裂缝建模和竞争粒子群算法的储层三三维地质建模模块，用于利用解释全层系地震波数据得到相关地质5测井曲线校正和质量分析模块，用于利用岩石力学和地学参数和地应力解释，并利用岩心力学参数开展单轴/三轴压缩试验和纵横声波速度测试三维属性模型模块，用于利用叠前地震反演建立研究区域的纵横波和天然裂缝的三维岩石力学场模型校正模块，用于根据天单井地层孔隙压力模型模块，用于分析研究区沉积规律和异常压局部地应力场方向获取模块，用于收集工区不同反演网格边界模块，用于对目标储层设置相应的载荷和约束竞争粒子群算法求解模块，用于运用竞争粒子群算行权利要求1～6任意一项所述基于确定性多尺度裂缝建模和竞争粒子群算法的储层三维得所述处理器执行权利要求1～6任意一项所述基于确定性多尺度裂缝建模和竞争粒子群67[0003]目前在三维地应力场的设置和模拟方面普遍存在如下几类问题:①较少的考虑断大优势，基于有限元网格的有限元地应力模拟可以较准确地获取目标区域的应力分布状赖于近井筒裂缝密度曲线，并利用蒙特卡洛等统计试验法实现天然裂缝随机分布处理假8[0012](1)现有技术没有充分利用地震属性信息对井间预测结果进行约束，特别是在裂[0013](3)现有技术多采用强制转换，或者剔除和牺牲局部节点数据的方法实现地质角自然邻点插值算法进行节点赋值，实现将地质模型属性结果传导至有限元的应力模型中，[0014](3)现有技术没有引入竞争粒子群优化算法进行三维地应力场求解计算，使得获9在地应力建模过程中开发的网格交互方案将有助于后续实现地质-油藏-压裂等多个模拟[0027]步骤五，通过从地表至目的层的全时窗叠前三维地震反[0035]在本发明一实施例中，所述步骤三利用三维地震进行区块天然裂缝的刻画和描MPa[0056]在本发明一实施例中，所述步骤七利用现场钻井或者测井实石力学参数和地应力解释，并利用岩心力学参数开展单轴/三轴压缩试验和纵横声波速度[0066]天然裂缝的三维岩石力学场模型校正模块，用于根据天[0068]网格转换算法应用模块，采用全局和局部双定位和跟踪方法寻找节点和单元信了竞争粒子群优化算法实现三维地应力场求解计算，较现有线性叠加地应力反演技术对[0080]图1是本发明实施例提供的基于确定性多尺度裂缝建模和竞争粒子群算法的储层[0083]图4是本发明实施例提供的考虑天然裂缝弱化的岩石力学数值模拟应用数据图。[0084]图5是本发明实施例提供的三维地应力场建模应用结果图。图5(a)最小水平主应[0089]如图1所示，本发明提供一种基于确定性多尺度裂缝建模和竞争粒子群算法的储力学参数和地应力解释，并利用岩心力学参数开展单轴/三轴压缩试验和纵横声波速度测层的形变控制在模型空间内，另外网格单元的位移和应力集中随着远离研究区而逐渐变[0100]本发明还提供一种基于确定性多尺度裂缝建模和竞争粒子群算法的储层三维应石力学参数和地应力解释，并利用岩心力学参数开展单轴/三轴压缩试验和纵横声波速度[0105]天然裂缝的三维岩石力学场模型校正模块，用于根据天[0107]网格转换算法应用模块，采用全局和局部双定位和跟踪方法寻找节点和单元信MPa[0137]步骤7：利用单井声波-电阻率和密度测井数据开展孔隙压力预测，优选伊顿法、层压力，如果平衡时间不足，得到的静压力可以作为地层压力的下限(特别是在低渗透地性则需要结合现场中途测试或者岩石物理实验得到其规新粒子群位置和速度调整更新粒子的适应度值获取全局最优的钻井诱导缝和井眼崩落裂缝走向FMI成像测井结果形成的各层位成像玫瑰图获取的地应者适应性变化遵循本公开的一般性原理并包括本公开未公开的本技术领域中的公知常识