基于球缝模型优化反演的裂缝储层评价新方法
基于FEM-DFN的页岩气井复杂裂缝扩展与优化——以长宁页岩气藏X1水平井组为例
基于FEM-DFN的页岩气井复杂裂缝扩展与优化——以长宁
页岩气藏X1水平井组为例
郑健;何永生;汪勇;蔡景顺;唐煊赫;朱海燕
【期刊名称】《断块油气田》
【年(卷),期】2024(31)3
【摘要】对长宁页岩气示范区水平井裂缝扩展形态认识不清、井间干扰频繁等问题严重制约了页岩气藏的高效开采。
对此,文中基于地质-工程一体化思想,结合适用于多井压裂分析的有限元-离散裂缝网络(FEM-DFN)方法,分析了水平井组体积裂缝扩展形态及天然裂缝对水力裂缝扩展的影响,厘清了井间干扰因素。
结果表明:在天然裂缝带发育区域,天然裂缝带对水力裂缝的诱导、捕获和阻隔作用主导了水力裂缝的扩展行为;天然裂缝带导致的水力裂缝不对称扩展、非均匀起裂是造成井间干扰的重要原因;天然裂缝带中各类裂缝的发育程度对水力裂缝扩展有重大影响。
该研究成果可为长宁页岩气藏高效水力压裂参数设计提供支撑和依据。
【总页数】9页(P415-423)
【作者】郑健;何永生;汪勇;蔡景顺;唐煊赫;朱海燕
【作者单位】四川长宁天然气开发有限责任公司;成都理工大学能源学院;成都理工大学油气藏地质及开发工程国家重点实验室
【正文语种】中文
【中图分类】TE319
【相关文献】
1.页岩气藏加密井压裂时机优化——以四川盆地涪陵页岩气田X1井组为例
2.基于复杂缝网模拟的页岩气水平井立体开发效果评价新方法——以四川盆地南部地区龙马溪组页岩气为例
3.基于FEM-DFN的页岩气储层水力压裂复杂裂缝交错扩展模型
4.基于停泵压降数据的页岩气井单段裂缝参数反演——以长宁N209井区页岩气井为例
因版权原因，仅展示原文概要，查看原文内容请购买。
裂缝等效宽度的斯通利波检测实验及反演
２０２４年３月第３９卷第２期西安石油大学学报（自然科学版）ＪｏｕｒｎａｌｏｆＸｉ’ａｎＳｈｉｙｏｕＵｎｉｖｅｒｓｉｔｙ（ＮａｔｕｒａｌＳｃｉｅｎｃｅＥｄｉｔｉｏｎ）Ｍａｒ．２０２４Ｖｏｌ．３９Ｎｏ．２收稿日期：２０２３ ０７ １３基金项目：国家自然科学基金“致密储层裂缝有效性的双尺度声波测井评价方法研究”（４２１０４１２６）第一作者：齐戈为（２０００ ），男，硕士研究生，研究方向：地球物理测井。
Ｅ ｍａｉｌ：２０２１７１０２９７＠ｙａｎｇｔｚｅｕ．ｅｄｕ．ｃｎ通讯作者：唐军（１９７９ ），男，博士，副教授，研究方向：地球物理测井和岩石物理。
Ｅ ｍａｉｌ：ｔａｎｇｊｕｎ＠ｙａｎｇｔｚｅｕ．ｅｄｕ．ｃｎＤＯＩ：１０．３９６９／ｊ．ｉｓｓｎ．１６７３ ０６４Ｘ．２０２４．０２．００４中图分类号：ＴＥ１３５文章编号：１６７３ ０６４Ｘ（２０２４）０２ ００３１ ０８文献标识码：Ａ裂缝等效宽度的斯通利波检测实验及反演齐戈为１，唐军１，蔡明１，何泽１，郑辰昌１，秦迎春２（１．长江大学地球物理与石油资源学院，湖北武汉４３０１００；２．加柏利能源技术（天津）有限公司，天津３００４５３）摘要：为了破解裂缝地层识别与评价的难题，进行缩比例模型井声波测量物理实验，结合数值模拟分析不同裂缝宽度和条数下斯通利波特征，总结水平裂缝不同宽度及条数对斯通利波衰减的影响规律，提出裂缝等效宽度的概念，定量表征裂缝带对斯通利波幅度的影响，建立基于裂缝等效宽度的反演方法和流程。
研究结果表明：随着裂缝等效宽度的增大，直达斯通利波幅度呈指数规律减小；在裂缝等效宽度保持不变的情况下，裂缝条数越多，直达斯通利波幅度越大；利用裂缝等效宽度能够实现裂缝带的定量表征。
关键词：声波测井；斯通利波；裂缝等效宽度；裂缝条数ＳｔｏｎｅｌｅｙＷａｖｅＤｅｔｅｃｔｉｏｎＥｘｐｅｒｉｍｅｎｔａｎｄＩｎｖｅｒｓｉｏｎｏｆＥｑｕｉｖａｌｅｎｔＦｒａｃｔｕｒｅＷｉｄｔｈＱＩＧｅｗｅｉ１，ＴＡＮＧＪｕｎ１，ＣＡＩＭｉｎｇ１，ＨＥＺｅ１，ＺＨＥＮＧＣｈｅｎｃｈａｎｇ１，ＱＩＮＹｉｎｇｃｈｕｎ２（１．ＣｏｌｌｅｇｅｏｆＧｅｏｐｈｙｓｉｃｓａｎｄＰｅｔｒｏｌｅｕｍＲｅｓｏｕｒｃｅｓ，ＹａｎｇｔｚｅＵｎｉｖｅｒｓｉｔｙ，Ｗｕｈａｎ，Ｈｕｂｅｉ４３０１００，Ｃｈｉｎａ；２．ＪｉａｂａｉｌｉＥｎｅｒｇｙＴｅｃｈｎｏｌｏｇｙ（Ｔｉａｎｊｉｎ）Ｃｏ．，Ｌｔｄ．，Ｔｉａｎｊｉｎ３００４５３，Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ：Ｔｈｅｐｈｙｓｉｃａｌｅｘｐｅｒｉｍｅｎｔｏｆａｃｏｕｓｔｉｃｗａｖｅｍｅａｓｕｒｅｍｅｎｔｗａｓｃａｒｒｉｅｄｏｕｔｉｎａｓｃａｌｅｄｍｏｄｅｌｗｅｌｌｉｎｏｒｄｅｒｔｏｉｄｅｎｔｉｆｙａｎｄｅｖａｌｕ ａｔｅｆｒａｃｔｕｒｅｆｏｒｍａｔｉｏｎ．ＴｈｅｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆＳｔｏｎｅｌｅｙｗａｖｅｓｕｎｄｅｒｄｉｆｆｅｒｅｎｔｗｉｄｔｈａｎｄｎｕｍｂｅｒｏｆｈｏｒｉｚｏｎｔａｌｆｒａｃｔｕｒｅｓｗｅｒｅａｎａｌｙｚｅｄｔｈｒｏｕｇｈｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎ，ａｎｄｔｈｅｉｎｆｌｕｅｎｃｅｓｏｆｔｈｅｗｉｄｔｈａｎｄｎｕｍｂｅｒｏｆｈｏｒｉｚｏｎｔａｌｆｒａｃｔｕｒｅｓｏｎｔｈｅａｔｔｅｎｕａｔｉｏｎｏｆＳｔｏｎｅｌｅｙｗａｖｅｓｗａｓｓｕｍｍａｒｉｚｅｄ．Ｏｎｔｈｉｓｂａｓｉｓ，ｔｈｅｃｏｎｃｅｐｔｏｆｅｑｕｉｖａｌｅｎｔｆｒａｃｔｕｒｅｗｉｄｔｈｗａｓｐｒｏｐｏｓｅｄｆｏｒｔｈｅｆｉｒｓｔｔｉｍｅ，ｔｈｅｉｍｐａｃｔｏｆｆｒａｃｔｕｒｅｚｏｎｅｓｏｎＳｔｏｎｅｌｅｙｗａｖｅａｍｐｌｉｔｕｄｅｗａｓｑｕａｎｔｉｔａｔｉｖｅｌｙｃｈａｒａｃｔｅｒｉｚｅｄ，ａｎｄａｎｉｎｖｅｒｓｉｏｎｍｅｔｈｏｄａｎｄｐｒｏｃｅｓｓｂａｓｅｄｏｎｅｑｕｉｖａｌｅｎｔｆｒａｃｔｕｒｅｗｉｄｔｈｗｅｒｅｅｓｔａｂｌｉｓｈｅｄ．Ｔｈｅｒｅｓｅａｒｃｈｒｅｓｕｌｔｓｓｈｏｗｔｈａｔａｓｔｈｅｅｑｕｉｖａｌｅｎｔｆｒａｃｔｕｒｅｗｉｄｔｈｉｎｃｒｅａｓｅｓ，ｔｈｅａｍｐｌｉｔｕｄｅｏｆｔｈｅｄｉｒｅｃｔＳｔｏｎｅｌｅｙｗａｖｅｄｅｃｒｅａ ｓｅｓｅｘｐｏｎｅｎｔｉａｌｌｙ．Ｕｎｄｅｒｔｈｅｓａｍｅｅｑｕｉｖａｌｅｎｔｆｒａｃｔｕｒｅｗｉｄｔｈ，ｔｈｅｍｏｒｅｔｈｅｎｕｍｂｅｒｏｆｆｒａｃｔｕｒｅｓ，ｔｈｅｇｒｅａｔｅｒｔｈｅａｍｐｌｉｔｕｄｅｏｆｔｈｅｄｉｒｅｃｔＳｔｏｎｅｌｅｙｗａｖｅ．Ｔｈｅｅｑｕｉｖａｌｅｎｔｆｒａｃｔｕｒｅｗｉｄｔｈｃａｎｂｅｕｓｅｄｔｏｑｕａｎｔｉｔａｔｉｖｅｌｙｃｈａｒａｃｔｅｒｉｚｅｔｈｅｆｒａｃｔｕｒｅｚｏｎｅ．Ｋｅｙｗｏｒｄｓ：ａｃｏｕｓｔｉｃｌｏｇｇｉｎｇ；Ｓｔｏｎｅｌｅｙｗａｖｅ；ｅｑｕｉｖａｌｅｎｔｆｒａｃｔｕｒｅｗｉｄｔｈ；ｎｕｍｂｅｒｏｆｆｒａｃｔｕｒｅｓ［Ｃｉｔａｔｉｏｎ］齐戈为，唐军，蔡明，等．裂缝等效宽度的斯通利波检测实验及反演［Ｊ］．西安石油大学学报（自然科学版），２０２４，３９（２）：３１ ３８．ＱＩＧｅｗｅｉ，ＴＡＮＧＪｕｎ，ＣＡＩＭｉｎｇ，ｅｔａｌ．Ｓｔｏｎｅｌｅｙｗａｖｅｄｅｔｅｃｔｉｏｎｅｘｐｅｒｉｍｅｎｔａｎｄｉｎｖｅｒｓｉｏｎｏｆｅｑｕｉｖａｌｅｎｔｆｒａｃｔｕｒｅｗｉｄｔｈ［Ｊ］．ＪｏｕｒｎａｌｏｆＸｉ’ａｎＳｈｉｙｏｕＵｎｉｖｅｒｓｉｔｙ（ＮａｔｕｒａｌＳｃｉｅｎｃｅＥｄｉｔｉｏｎ），２０２４，３９（２）：３１ ３８．西安石油大学学报（自然科学版）引 言裂缝作为致密储层重要的渗流通道和储集空间，是勘探与开发研究的重点［１］。
基于叠前AVAZ反演的火山岩裂缝预测方法研究
（b）断裂相干属性 图4 裂缝与断层分布
将本文方法与椭圆拟合法的裂缝预测结果 进行对比，椭圆拟合法预测结果显示裂缝在研究 区中部较发育，和目的层均方根振幅属性的强振 幅规律基本一致，而在北部断裂带附近无明显异 常，说明椭圆拟合法较易受振幅影响。相比于椭 圆拟合方法，本文叠前AVAZ反演方法裂缝预测结 果更加突出了强各向异性的裂缝发育区，减小了 强振幅对预测结果的影响，对断裂带附近的裂缝
图2 不同入射角下振幅随方位角的变化
2.2 叠前 AVAZ 反演及应用效果分析
研究区C5井成像测井显示，目的层火石岭组 井周裂缝较为发育，主要包括高导缝、诱导缝和 高阻缝。裂缝走向玫瑰图见图3，三种类型裂缝的 走向都为近南-北向，近似垂直于最大水平主应力 方向。
基于校正后的OVT道集开展叠前AVAZ反演， 得到各向同性梯度、各向异性梯度（裂缝发育密 度）和方位角（裂缝发育方向）数据。图4a为沿 目的层裂缝发育矢量图，颜色代表裂缝发育密 度，矢量棒方向代表裂缝方向，图4b为沿层提取 的断裂相干属性。图4显示裂缝与断层分布相关
Ｓｔｕｄｙ ｏｎ ｖｏｌｃａｎｉｃ ｒｏｃｋ ｆｒａｃｔｕｒｅ ｐｒｅｄｉｃｔｉｏｎ ｍｅｔｈｏｄ ｂａｓｅｄ ｏｎ ｐｒｅｓｔａｃｋ ＡＶＡＺ ｉｎｖｅｒｓｉｏｎ Ｂｉ Ｃｈｅｎｃｈｅｎ，Ｙａｎｇ Ｑｉｎｌｉｎ，Ｘｉｅ Ｗｅｉ，Ｌｉａｎｇ Ｚｈｉｑｉａｎｇ，Ｃｈｅｎｇ Ｙｕｎ Sinopec Geophysical Research Institute Co.，Ltd.，Nanjing 211103
1 叠前 AVAZ 反演方法原理
根据等效介质理论，可将定向发育裂缝的地 层可看作HTI介质模型，在实际中地震波反射系数 与反射角和传播方向均有关。当来自各个方向的 地震波穿过裂缝时，其振幅、速度、旅行时差都
岩石物理驱动的储层物性参数非线性地震反演方法
岩石物理驱动的储层物性参数非线性地震反演方法潘新朋;刘志顺;高大维;王璞;郭振威;柳建新【期刊名称】《地球物理学报》【年(卷),期】2024(67)3【摘要】叠前地震反演和岩石物理反演分别是获取弹性参数和物性参数的重要手段,两者结合有助于实现储层参数预测并精细刻画储层特征.储层物性参数的反演依赖于岩石物理模型,在进行物性参数反演时可以将复杂的岩石物理模型做泰勒展开,进而得到其一阶或高阶的近似表达式,然而这会降低模型的精确性并增加反演的误差.为了提高储层物性参数反演的稳定性和准确性,本文以碎屑岩储层为例,提出了岩石物理驱动的储层物性参数非线性地震反演方法.首先,基于贝叶斯框架和高斯分布约束条件,从叠前地震数据中实现纵、横波速度及密度等弹性参数的反演.其次,通过碎屑岩岩石物理模型建立起弹性参数与物性参数之间的联系.最后,利用粒子群算法进行全局寻优获得较为准确的孔隙度、泥质含量和含水饱和度等物性参数.合成数据和实际资料测试结果验证了所提方法的可行性和准确性,反演结果与测井数据吻合较好,可有效指示含气储层区域,本文方法在储层预测和评价方面具有广泛的应用前景.【总页数】18页(P1237-1254)【作者】潘新朋;刘志顺;高大维;王璞;郭振威;柳建新【作者单位】中南大学地球科学与信息物理学院;有色资源与地质灾害探查湖南省重点实验室;有色金属成矿预测与地质环境监测教育部重点实验室【正文语种】中文【中图分类】P631【相关文献】1.岩石物理驱动的储层裂缝参数与物性参数概率地震反演方法2.储层弹性与物性参数地震叠前同步反演的确定性优化方法3.高分辨率非线性储层物性参数反演方法和应用4.地震和测井资料联合反演储层物性参数的方法5.地震岩石物理驱动的裂缝预测技术研究现状与进展(Ⅰ)——裂缝储层岩石物理理论因版权原因，仅展示原文概要，查看原文内容请购买。
储层裂缝的综合表征技术
储层裂缝的综合表征技术储层裂缝是指沉积岩层中的裂缝或缝隙，它们对油气储层的渗透性和储集性能有着重要影响。
因此，准确地综合表征储层裂缝的特征是油气勘探开发中的关键问题之一。
本文将介绍几种常用的储层裂缝综合表征技术。
一维裂缝参数表征技术是一种简单而有效的方法。
它通过测量岩心样品中裂缝的长度、宽度和密度等参数，来表征裂缝的分布和特征。
这种方法适用于储层裂缝较少且分布较均匀的情况。
通过对一系列岩心样品的测试和分析，可以获得储层裂缝的统计特征，进而对裂缝的形成机制和储层的渗透性进行定量评价。
声波测井是一种常用的储层裂缝表征技术。
它利用声波在岩石中的传播特性，通过测量声波数据来推断储层中裂缝的存在和性质。
声波测井可以提供储层裂缝的方位、长度、宽度和密度等信息，从而为储层评价和开发提供重要依据。
同时，声波测井还可以结合其他地球物理数据，如电阻率测井和密度测井等，来进行综合解释和分析，提高裂缝表征的准确性和可靠性。
地震反演技术也被广泛应用于储层裂缝的综合表征。
地震反演是一种将地震数据与地下介质模型进行匹配的方法，通过反演得到地下介质的物理参数，如波速和密度等。
在储层裂缝表征中，地震反演可以提供裂缝的几何形状、分布范围和孔隙度等参数。
这种方法适用于储层裂缝较多、分布较复杂的情况。
然而，地震反演技术需要大量的地震数据和高性能计算设备支持，所以在实际应用中具有一定的限制。
数值模拟方法是一种基于物理模型的储层裂缝表征技术。
它通过建立储层裂缝的数值模型，利用数值计算方法求解模型中的物理过程和现象，从而得到裂缝的形态和分布规律。
数值模拟方法可以模拟储层裂缝的形成、演化和影响等过程，为储层开发和优化提供科学依据。
然而，数值模拟方法需要考虑多种物理过程和复杂的计算模型，所以在实际应用中需要具备一定的专业知识和技术水平。
储层裂缝的综合表征技术包括一维裂缝参数表征、声波测井、地震反演和数值模拟等方法。
每种方法都有其适用的场景和优势，可以根据实际需求选择合适的技术手段来进行储层裂缝的综合表征。
微地震方法的裂缝监测与储层评价
微地震方法的裂缝监测与储层评价李政;常旭;姚振兴;王一博【摘要】在油田水力压裂微地震事件定位结果的基础上,结合有效微地震事件的时空分布、震级大小、地震矩、震源半径、应力降和b值等地震学参数进行综合研究,并结合研究区域的地质背景和测井资料对水力压裂诱发的裂缝网络进行几何形态分析和应力解释.本文提出的微地震综合分析解释方法可对压裂后储层物性进行综合评价,有利于对储层改造效果进行预测,对油田的水力压裂施工具有指导意义.【期刊名称】《地球物理学报》【年(卷),期】2019(062)002【总页数】13页(P707-719)【关键词】微地震;震源参数;b值;裂缝解释;储层评价【作者】李政;常旭;姚振兴;王一博【作者单位】中国科学院地质与地球物理研究所,北京 100029;中国科学院地球科学研究院 ,北京 100029;中国科学院大学,北京 100049;中国科学院地质与地球物理研究所,北京 100029;中国科学院地球科学研究院 ,北京 100029;中国科学院地质与地球物理研究所,北京 100029;中国科学院地球科学研究院 ,北京 100029;中国科学院地质与地球物理研究所,北京 100029;中国科学院地球科学研究院 ,北京 100029【正文语种】中文【中图分类】P6310 引言页岩油气资源开发的核心技术是水平井以及水压致裂技术，其中采用微地震监测手段识别水力压裂裂缝的分布和发育尤为重要.然而，在实际生产中，微地震技术的应用并不完善，利用微地震所能获取的各种信息对压裂裂缝的识别以及对压裂后的储层评价也一直是一个需要研究的问题.20世纪80年代，Thill(1972)研究了脆性岩石微破裂诱发的地震信号，Schuster(1978)将水力压裂微地震监测技术引入了石油领域，随即被迅速应用于非常规油气的压裂效果监测中(Maxwell and Urbancic, 2001; Shapiro et al., 2002).近年来，微地震监测的理论方法研究得到了快速发展，在众多的方法研究中，与压裂裂缝解释密切相关的是震源定位和震源机制反演.在震源定位方法的研究中，地震波的干涉原理和衰减特性得到了利用(王璐琛等, 2016; 常旭等, 2018)，另一方面，通过微地震和地面观测数据的联合反演也提高了震源定位的精度(刁瑞等, 2017)，震源定位精度的提高直接影响储层压裂裂缝的解释.目前在微地震反演的实际应用研究中，震源位置数据得到普遍的应用(Eisner et al., 2010;Maxwell,2010)，而对震源机制以及与震源机制相关的地震学数据的应用尚不充分.事实上，微地震的震源机制与储层介质的特性、裂缝发育的产状、地层应力的方向有直接关系(翟鸿宇等, 2016;李晗和姚振兴, 2018， Chen et al.,2018)，利用微地震反演获得的地震学参数对压裂效果进行评价和解释至关重要.传统的反射地震资料解释主要根据地震波的传播规律和地质特征，把各种地震波信息转变为构造、地层岩性等信息，将地震剖面转变为地质剖面进行解释.对于非常规油气的勘探开发，需要实施储层压裂改造后的裂缝解释与评价，压裂施工中的微地震监测资料提供了不同于反射地震资料的储层评价信息.微地震数据利用地震学的研究方法得到震源的时间空间位置、震源机制等信息，利用这些信息同样可以对储层的物理性质进行分析和解释.因此，本文研究了微地震资料的储层解释方法，以某油田实际资料为基础，提出了微地震数据的分析方法和流程，在压裂区微地震事件定位的基础上，对有效微地震事件的时空分布、震级、地震矩、震源半径、应力降和地震b值等地震学参数进行综合研究，并结合研究区域的地质背景、测井等资料，对压裂区的水力压裂裂缝进行了分析，对研究区的储层性质与水力压裂效果进行了评价.1 研究区域的地质背景及微地震监测数据本文的研究区域位于鄂尔多斯盆地南部，按地层的分布形态划分为：北部伊盟隆起、南部渭北隆起、西部西缘断褶带和天环坳陷(天环向斜)、东部晋西挠褶带以及中部的陕北斜坡(西倾单斜构造)等五个一级构造单元.盆地内致密油层的主要特征为：(1)致密油藏源储互层共生，平面上主要分布在湖盆中部有利砂体中.主力层段受多级坡折带影响，重力流发育，指示了储层的非均一性特征；(2)储层孔喉结构复杂，实验室可见孔隙含量低，以微(纳米)孔隙类型为主，指示了储层的致密特征；(3)储层致密砂岩天然微裂缝发育(每10 m发育天然裂缝约2～3条)，水平两向应力差4～7MPa，有利于在水力压裂作用下形成复杂裂缝网络(杨华等, 2015).研究区内压裂井与监测井井位如图1所示.该井有效储层纵、横向变化较大，非均质性强，埋藏深度为1640 m，平均油层厚度为11.9 m，平均孔隙度为9.7%，填隙物含量为14.7%，是典型的致密储层.通过对岩心裂缝的测量统计和测井资料的构造裂缝发育的研究可知，水平压裂井所在地区裂缝密度较大(牛小兵等, 2014).图1 压裂井和监测井井位(a)俯视图； (b) 侧视图.Fig.1 Horizontal well and monitoring well location(a) Plan view; (b) Side view.施工单位在本文研究区域内展开了“水力喷砂分段多簇混合水力压裂”储层改造实验.实验中，压裂井垂深1588 m，水平段长度为1205 m，改造段数为14段，每段各有两个射孔.施工排量为6 m3·min-1(油套同注)，平均套管迫压为20.7 MPa，单级入地液量为800 m3，施工入地总液量为11200 m3，支撑剂总砂量为600 m3，单级砂量为43 m3，其中喷射用20/40目石英砂为3 m3，压裂用100目粉陶为1.5 m3，40/70目覆膜砂为30 m3，20/40目覆膜砂为12 m3，最高砂比28%，平均砂比7.2%.微地震监测井共设置12个三分量检波器，对水力压裂过程中形成的微地震事件进行记录.技术人员对采集到的14段压裂中的微地震事件进行了识别和预处理，提供了5832个具有企业标准的微地震事件的信息.图2 采用的微震事件数量与震级Fig.2 The number of Micro-seismic events and its magnitude2 微地震数据对压裂裂缝的几何解释2.1 微震事件的时空分布与裂缝的关系对这些微地震数据的各项参数进行了再分析，对水力压裂储层改造的效果进行了评价.图2为微地震事件按照震级大小统计的全压裂井段微地震事件柱状图.由图2可知，绝大多数的微地震事件震级集中在-2.8级至-2级之间.图3是用不同颜色表示的14个压裂段微震事件的空间位置在水平面的投影，图中由蓝色到红色的每一种颜色代表一个压裂段.本文对14段中已完成定位的5832个微地震事件进行了综合分析，根据微地震事件的时空信息、震级大小、震源半径、地震矩、应力降、b值等地震学参数进行了研究，并结合研究区的地质资料对研究区压裂裂缝的分布和发育进行了半定量分析与定性解释.微地震事件的空间分布可以提供裂缝的空间位置.由图3可知，各压裂段监测到的微地震事件沿水平井两侧分布，且分布半径约400 m，预示压裂形成的裂缝分布在以水平井为轴线，半径约400 m的椭圆柱形空间范围内.根据色标可知，深蓝色为第1段压裂段，深红色为第14段压裂段.由图3可知，每一段压裂的微震事件可能出现在相邻段内，这表明在对储层进行压裂改造的过程中可能激活了相邻未压裂段的原生裂缝，使得在相邻压裂段中产生了微地震事件，因此，根据微地震事件的时空分布可以定性分析压裂裂缝的变化.图3表现出的微地震事件时空分布说明该井的压裂对储层产生复杂联通的裂缝网络有利.在全部压裂段中，第六压裂段微地震事件分布表现出明显的异常.本文看到除了主裂缝(图3中蓝线所示)以外，还产生了一条与主裂缝相交的次生裂缝(图3中红线所示).在对压裂事件生长过程的研究中，本文发现第五段压裂过程中，尚未实施压裂的第六段区域在沿着不同于主裂缝的方向上产生了数个震级较大的微地震事件(例如图3子图中红色箭头所示).在第五段压裂结束后，第六段压裂首先沿着之前产生过大微震事件的方向(图4)开始破裂，然后再沿主裂缝方向上产生大量微地震事件.图4是第5段压裂微地震事件按照时间顺序排列的分布图，色标右侧为压裂时间，色标左侧为产生的微地震事件个数.由图4可以看出，第六段压裂共产生微地震事件694个，其中有近300个微地震事件在22 min左右产生且分布位置基本沿着次生裂缝的方向，剩下近400个微地震事件在约3 h之内产生且基本沿着主裂缝方向分布.第六段微震事件的时空分布特征说明了第六段区域存在着与主裂缝走向不一致的原生裂缝，此裂缝极有可能在第五段的压裂过程中被提前激活，从而产生了交叉的裂缝网络.2.2 微震事件震级大小与裂缝的关系为了分析微地震事件与水力压裂缝之间的关系，本文将所有微地震事件按照震级的大小进行排列.图5是全压裂段微地震事件震级大小及空间位置分布，从蓝色到红色分别表示了微地震震级从小到大的趋势.根据图5可知，微地震事件的震级从-3.2至-1.2之间不等，大部分微震事件的震级集中在-2.8至-2之间.从微地震事件的整体分布上可以看出，压裂区域中部震级相对较小且分布非常集中，两侧压裂段的微震事件震级相对中部较大且分布相对分散.根据研究区储层具有较强非均匀性的特点，可定性解释为压裂区储层中部岩性泥质含量较高，脆性较差，不易破裂，而两侧储层岩性泥质含量较低，脆性较好，较易破裂.图3 各压裂段微震事件时间空间分布Fig.3 Spatial and temporal distribution of micro-seismic events in fracturing segments图4 第六段压裂事件时间空间分布图(按微震事件发生顺序排列)Fig.4 Spatial and temporal distribution belonging to the sixth section micro-seismic events (in order of micro-seismic events)2.3 裂缝走向及破裂面产状根据各压裂段微地震事件的空间位置，对每个压裂段的微地震事件进行了三维空间直线拟合与裂缝走向分析，同时根据主裂缝的长度和方位角拟合出每个压裂段的破裂面(图6).表1为各段裂缝的长度、方位、角度、破裂面的面积以及方位、角度.由此，本文可以得到，除了第六段，其余压裂段的破裂面近似为铅锤面，而第六段的破裂面近似为水平面.根据破裂面产状的差异以及该区域储层裂缝发育的特点可以推测，第六段的破裂面与其他各段破裂面的形成机制不同，其他各段的破裂主要由压裂形成，而第六段的破裂可能与原生构造相关，如果第六段所在的储层存在原生裂缝，则可能在实施第五段压裂时被激活.这一推测与图3的微地震发震时间的分布吻合，即尚未实施第六段压裂时，第六段上的储层岩石已经发生破裂并产生了大震级的微地震事件.图5 全压裂段微震事件的震级与空间位置分布Fig.5 Magnitude and spatial distribution of micro-seismic events in full fracturing图6 各压裂段裂缝主破裂面Fig.6 Main fracture surface of fracturing sections 所以图6可以解释为第五段的压裂过程提前激活了第六段的原生裂缝，进而导致两个相邻压裂段的裂缝形成了新的交错网络，两段裂缝主破裂面连通后形成了一个近乎水平的破裂面.2.4 压裂体积估算所谓压裂体积又称改造体积，是指在水力压裂过程中，脆性岩石产生剪切滑移诱使天然裂缝不断扩张，继而形成裂缝网络，增加改造体积，提高产量和最终采收率.对每个压裂段的微地震事件，将其点集所构成的三维凸包作为该压裂段的改造体，将三维凸包的体积作为压裂体积；将三维点集用最小二乘法拟合得到的平面作为主破裂面，主破裂面与三维凸包的截面作为主破裂面面积；将三维点集用最小二乘拟合得到的直线作为主破裂线，主破裂线被主破裂面截断的长度作为主裂缝的破裂长度.对每段压裂的有效微地震三维数据体构建了多面体包络面，从而得到了压裂井每段的压裂体积(图7)，具体数值见表1.表1 各压裂段裂缝参数裂统计Table 1 Statistics of fracturing results in 14 fracturing sections压裂段改造体积(106m3)破裂面积(105m2)破裂面走向(°)破裂面倾角(°)主裂缝长度(m)主裂缝走向(°)主裂缝倾角(°)18.20051.976283.072472.7558753.8282.91790.648026.57661.556981.399 381.528813.6581.71993.152936.37421.934179.631984.4782844.3579.67771 .161445.47271.506878.997577.5754788.0779.16300.512655.86641.599778. 23285.5345927.6478.38720.37863.49591.212884.19246.6331939.4978.5228 0.75926-10.36800.289784.19246.6331385.85103.20460.715473.63061.305878.69818 1.6197879.7678.37020.321484.11121.206477.810975.6734738.02101.46962 .867993.59691.114679.416770.7342718.2678.95071.3916104.96021.566874 .710578.5051765.33105.12011.0948112.3920.876073.660988.7193689.1073 .83711.3046123.92791.213174.485973.8945759.5775.26922.5781132.02310 .878777.492088.1805659.3777.46200.3995143.96901.228076.660282.49647 73.7276.83360.6338图7 各压裂段压裂体积分布Fig.7 Fracturing volume distribution of 14fracturing sections3 微地震裂缝的应力解释微地震裂缝的应力解释指在完成研究区域水平压裂段微地震事件定位工作的基础上，利用微地震事件的地震矩、应力降、震源半径、b值等地震学参数进行综合研究，对储层的力学性质进行解释.3.1 地震矩、应力降与震源半径及其应力解释地震矩(M0)是地震释放能量大小的直接度量(Kanamori,1977)，其由震源位错理论给出定义：其中μ为剪切模量，s为断裂面积，为平均位错量.由虎克定律给出应力降(Δσ)与应变(Δe)之间的关系：其次，利用布龙(Brune)圆盘形位错模式，震源等效圆位错半径(r)和拐角频率(f0)的关系为β是S波波速(Brune, 1970; Svanes, 1971).将研究区的全部微地震事件的地震矩，即三维空间分布的能量点集，进行了三次样条插值处理，然后沿水平井所在的水平面做三维数据体的截面，得到该区地震矩的分布图(图8).假设在本次压裂实验中，地下储层空间得到了充分的水力压裂，则图8所示的结果即表示该研究区可能释放能量的整体分布情况.由此可以发现储层空间中部释放能量较少，而其周围释放能量较多.与微地震震级分布(图5)、压裂缝长度和破裂面面积(图6)的分析结果比较，地震矩的三维空间分析(图8)能更好地评价储层空间的三维连续性变化.与表示地震释放能量的参数地震矩相比，地震应力降评价的则是地震断层机制和地下块体释放能量的行为.微地震事件中的应力降应特指岩石受到应力急剧变化时破裂所造成的应力降低的行为，水力压裂过程中发生的微地震事件的应力降则多在2000～10000 Pa之间.在常温高压下，岩石破裂应力降越大，岩石破裂强度也越高(臧绍先, 1984) .本文对研究区内14个压裂段的微地震应力降数据同样进行了三维空间的数据插值处理，得到了水平井所在水平面上的三维数据体横截面(图9).结果显示，研究区中部应力降较小，而其周围的应力降反而较大.这与图8中显示的地震矩分布图情况类似，然而应力降中心的低值区域更为集中.其是否能对破裂区的范围有更好的指示意义还需进一步研究.本文同时也对震源半径进行了和地震矩、应力降两个参数相同的插值处理，图10是震源半径所得到的结果.与图8和图9相同，图10中间除了极个别大震级的微地震事件是高值以外，研究区域整体中间为低值，两侧为高值区.3.2 b值及其应力解释1941年Gutenberg和Richter通过对大量地震资料研究发现：地震震级M与大于等于震级M的地震数目N之间存在如下公式所示的幂律分布关系(Gutenberg and Richter, 1942)，即G-R关系：logN=a-bM式中，a，b为常数，是描述地震带内地震震级频度分布特征的重要参数，可以根据地震资料通过公式计算获得.其中，a反映平均地震活动水平；b反映大小地震的比例关系.在地震预报领域中，研究发现：大震前震源及附近区域经常会出现某些震级档内的地震增多或减小，导致出现大小地震比例失调，b值减小的异常现象，此外，区域应力积累水平升高是大地震发生的必要条件，因此认为，b值反映了地应力状态，二者呈反比关系，且b值的高低与岩石介质的特性有关，比如岩石的脆性、弹性、塑性、破裂程度等(韩骏和姚令侃, 2015).岩石学实验中，声发射活动与地震活动的机制最为接近，在统计参数上与地震活动性的可对比性也最强.Scholz通过岩石破裂实验发现，岩石随着应力的增加，b值出现明显下降(Scholz, 1968).通过对5种岩石测量记录曲线和数据进行分析处理，得到与Scholz实验类似的结果，在岩石达到破裂应力之前b值下降较快.Wyss等(2000)研究矿山岩石的破裂行为发现：岩体内构造应力的大小与b值图8 研究区地震矩Fig.8 Seismic moment in study area图9 研究区应力降Fig.9 Stress drop in study area图10 研究区微地震事件震源半径Fig.10 Source radius of the micro-seismic events in study area成反比，低b值区往往具有更高的应力积累.本文对研究区域的14个水平压裂段进行了b值的计算.首先将每段的水平压裂的微地震事件个数按照震级大小进行分段统计.根据最小二乘法，分别算出每段的b 值.图11为14个压裂段分别统计的b值.其中，位于整个压裂段中部的第7、8、9段的b值分别为2.59、2.36、2.56，明显大于两侧压裂段的b值(两侧b值平均值为1.95).根据岩石构造与b值的关系，本文推测，整个压裂段中部b值较高的区域对应的应力较小，这一推测与研究区中心部位微震事件的震级小于周边区域的现象可相互印证.图12为全压裂段微地震事件的b值，为2.44，相对两侧压裂段b 值较大.图11 各压裂段b值(1-8段)Fig.11 b-value of 1-8 fracturing sections图11 (续) 各压裂段b值(9-14段)Fig.11 (continued) b-value of 9-14 fracturing sections图12 全部微地震事件b值Fig.12 b-value from all micro-seismic events研究区划分为20 m×20 m的网格点，对每个网格点四周50 m×70 m的区域内的微震事件进行统计并计算b值，计算结果作为该格点处的地震b值，然后对所有格点的b值进行插值处理.为保证统计计算中有足够的样本量使计算结果稳定且可靠，每个单元格内的微地震数量最少为15个.由此，本文得到了研究区域的b 值分布图(图13).根据b值分布图以及图11各段b值分析，本文认为b值小于2.3的区域应力降较大，属于潜在破裂区.由图13可知，通过b值扫描，本文将分散的微震事件的分布转换成了储层连续受力情况的分布，这对储层裂缝的区域连通性有非常好的指示意义并对压裂施工有着重要的指导意义.图13 研究区b值分布Fig.13 b-value in study area图14 研究区水力压裂与测井资料联合解释(a) 测井解释，油田提供； (b) 微地震监测.Fig.14 Joint interpretation of hydraulic fracturing and logging data in the research area(a) Logging data interpretation, provided by oil field; (b) Micro-seismic monitoring interpretation.4 微地震与测井资料的综合解释将微地震裂缝的定性与定量分析图件和测井解释图件结合(图14)，本文可以对储层的岩性与裂缝分布有更好的约束和指示.压裂区域中部微地震事件震级相对较小且分布紧凑，两侧压裂段震级偏大且分布相对分散.根据声波(AC)以及自然伽马(GR)测井曲线可以得出压裂井段中部储层岩石泥质含量高，脆性较差，孔隙度较小，通过含烃曲线可以看出该部分含油气少，储层较差，而压裂井段两侧岩石泥质含量低，脆性较好，孔隙度较大，油气含量大，储层性质较好.5 结论本文提出了利用微地震事件的地震学信息进行储层评价与解释的方法：(1)微地震事件的时间-空间分布可以定性分析压裂裂缝的时间-空间分布，还可以定量计算压裂裂缝的长度、方位角、破裂面的面积和储层的改造的体积.对事件时空分布的异常点进行重点研究，可以进一步推测储层中原生裂缝与压裂裂缝的关系；(2)微地震事件的震级、震级大小的分布可以指示储层岩性的非均一性，对储层岩性和脆性作出定性分析；(3)微地震的地震矩、应力降、震源半径的研究可以分析储层压裂后岩石破裂程度以及能量传播的方向；(4)微地震b值可以对地应力状态进行研究，b值扫描可以对储层潜在的破裂区进行预测，可以将分散的微震事件的分布转换成储层连续受力情况的分布，这对储层裂缝的连通性以及水力压裂方案的设计有重要的指导意义.本文提出的研究方法用于水力压裂微地震监测实际资料的解释，得出的结果不仅与测井资料给出的岩性信息吻合，而且给出了储层裂缝发育和应力分布的新的认识，获得了储层裂缝的几何形态、原生裂缝与压裂裂缝的关系、储层改造体积、应力状态的分布以及潜在的破裂区域.致谢感谢长庆油田对本文的研究提供的帮助和支持.References【相关文献】Brune J N. 1970. Tectonic stress and the spectra of seismic shear waves from earthquakes. Journal of Geophysical Research: Atmospheres, 75(26): 4997-5009.Chang X, Li Z, Wang P, et al. 2018. Micro-seismic location based on frequency attenuation compensation. Chinese Journal of Geophysics (in Chinese), 61(1): 250-257, doi:10.6038/cjg2018L0396.Chen H C,Meng X B, Niu F L, et al. 2018. Microseismic monitoring of stimulating shale gas reservoir in SW China: 2. spatial clustering controlled by the preexisting faults and fractures. Journal of Geophysical Research: Solid Earth, 123(2): 1659-1672.Diao R, Wu G C, Shang X M, et al. 2017. Joint correction method based on 3D seismic and surface microseismic data. Chinese Journal of Geophysics (in Chinese), 60(1): 283-292, doi: 10.6038/cjg20170123.Eisner L, Williams-Stroud S, Hill A, et al. 2010. Beyond the dots in the box: Microseismicity-constrained fracture models for reservoir simulation. The Leading Edge,29(3): 326-333.Gutenberg B, Richter C F. 1942. Earthquake magnitude, intensity, energy, and acceleration. Bulletin of the Seismological Society of America, 32(3): 163-191.Han J, Yao L K. 2015. Crustal stress evaluation method based on seismic activity parameter b-value. Railway Standard Design (in Chinese), 59(7): 36-39, 127.Kanamori H. 1977. The energy release in great earthquakes. Journal of Geophysical Research, 82(20): 2981-2987.Li H, Yao Z X. 2018. Microseismic focal mechanism inversion in frequency domain based on general dislocation point model. Chinese Journal of Geophysics (in Chinese), 61(3): 905-916, doi: 10.6038/cjg2018L0237.Maxwell S C,Urbancic T I. 2001. The role of passive microseismic monitoring in the instrumented oil field. The Leading Edge, 20(6): 636-639.Maxwell S C, Rutledge J, Jones R, et al. 2010. Petroleum reservoir characterization using downhole microseismic monitoring. Geophysics, 75(5): 75A129-75A137.Niu X B, Hou G T, Zhang J Z, et al. 2014. Assessment criteria for member 6 and 7 tight sand in Ordos basin and its applications. Geotectonica et Metallogenia (in Chinese), 38(3): 571-579.Scholz C H. 1968. The frequency-magnitude relation of microfracturing in rock and its relation to earthquakes. Bulletin of the Seismological Society of America, 58(1): 399-415. Schuster C L. 1978. Detection within the wellbore of seismic signals created by hydraulic fracturing. ∥48th Ann. Internat Mtg., Soc. Expi. Geophys.. Expanded Abstracts.Shapiro S A, Rothert E, Rath V, et al. 2002. Characterization of fluid transport properties of reservoirs using induced microseismicity. Geophysics, 67(1): 212.Svanes K. 1971. Correction [to “Tectonic stress and the spectra, of seismic shear waves from earthquakes”]. Journal of Geophysical Research, 76(20): 5002-5002.T hill R E. 1972. Acoustic methods for monitoring failure in rock. ∥Proceedings of the 14 U.S. Symposium on Rock Mechanics. New York: American Rock Mechanics Association, 1972.Wang L C, Chang X, Wang Y B. 2016. Locating micro-seismic events based on interferometric traveltime inversion. Chinese Journal of Geophysics (in Chinese), 59(8): 3037-3045, doi: 10.6038/cjg20160826.Wyss M, Schorlemmer D, Wiemer S. 2000. Mapping asperities by minima of local recurrence time: San Jacinto-Elsinore fault zones. Journal of Geophysical Research: Solid Earth, 105(B4): 7829-7844.Yang H, Liu X S, Yan X X. 2015. The relationship between tectonic-sedimentary evolution and tight sandstone gas reservoir since the late Paleozoic in Ordos Basin. Earth Science Frontiers (in Chinese), 22(3): 174-183.Zang S X. 1984. Earthquake stress drop and the stress drops of rock fracture. Acta Seismologica Sinica (in Chinese), 6(2): 182-193.Zhai H Y, Chang X, Wang Y B, et al. 2016. Inversion for microseismic focal mechanisms inattenuated strata and its resolution. Chinese Journal of Geophysics (in Chinese), 59(8): 3025-3036, doi: 10.6038/cjg20160825.附中文参考文献常旭, 李政, 王鹏等. 2018. 基于频率衰减补偿的微地震定位方法.地球物理学报, 61(1): 250-257, doi: 10.6038/cjg2018L0396.刁瑞, 吴国忱, 尚新民等. 2017. 三维地震与地面微地震联合校正方法. 地球物理学报, 60(1): 283-292, doi: 10.6038/cjg20170123.韩骏, 姚令侃. 2015. 基于地震活动性参数b值的地应力评估方法研究. 铁道标准设计, 59(7): 36-39, 127.李晗, 姚振兴. 2018. 基于“剪切+张裂”一般位错模型频率域求解微震震源机制. 地球物理学报, 61(3): 905-916, doi: 10.6038/cjg2018L0237.牛小兵, 侯贵廷, 张居增等. 2014. 鄂尔多斯盆地长6-长7段致密砂岩岩心裂缝评价标准及应用. 大地构造与成矿学, 38(3): 571-579.王璐琛,常旭, 王一博. 2016. 干涉走时微地震震源定位方法. 地球物理学报, 59(8): 3037-3045, doi: 10.6038/cjg20160826.杨华, 刘新社, 闫小雄. 2015. 鄂尔多斯盆地晚古生代以来构造-沉积演化与致密砂岩气成藏. 地学前缘, 22(3): 174-183.臧绍先. 1984. 地震应力降与岩石破裂应力降. 地震学报, 6(2): 182-193.翟鸿宇,常旭, 王一博等. 2016. 含衰减地层微地震震源机制反演及其反演分辨率. 地球物理学报,59(8): 3025-3036, doi: 10.6038/cjg20160825.。
基于天然裂缝破坏行为的页岩储层压裂微地震事件预测
基于天然裂缝破坏行为的页岩储层压裂微地震事件预测李强;尹成;王俊力;罗浩然【摘要】微地震事件的分布范围和密度与页岩气井的产量有较好的正相关关系,经济、高效地进行微地震事件的压前预测和压后评估,对压裂增产优化设计具有重要的指导意义.针对微地震监测技术成本高和施工受限的问题,基于渗流力学和岩石力学理论,结合天然裂缝张性和剪切破坏准则,建立了水力裂缝动态扩展过程中微地震事件模拟预测的流固耦合数学模型,给出了数值求解方法.采用涪陵页岩气示范区HF-X水平井的储层地质参数和施工参数,对页岩储层压裂过程中剪切和张性微地震事件的动态演化和展布进行数值模拟,模拟结果与实时监测的微地震数据吻合度高于80％.结果表明,该理论方法能够经济、高效和可靠地预测微地震事件的展布.【期刊名称】《石油物探》【年(卷),期】2018(057)006【总页数】7页(P878-883,891)【关键词】页岩;压裂;微地震;破坏准则;动态演化;数值模拟;压前预测;压后评估【作者】李强;尹成;王俊力;罗浩然【作者单位】西南石油大学地球科学与技术学院 ,四川成都 610500;西南石油大学地球科学与技术学院 ,四川成都 610500;中国石油天然气股份有限公司西南油气田分公司蜀南气矿 ,四川泸州 646000;西南石油大学地球科学与技术学院 ,四川成都610500【正文语种】中文【中图分类】P631页岩储层渗透率极低,但通常发育有大量天然裂缝,水平井多段分级压裂使储层产生复杂裂缝网络,为气体流动提供“高速通道”,是实现页岩气工业开采的有效技术[1-2]。
伴随着大量低粘液体在大排量的条件下注入,多条水力裂缝同时延伸,压裂液向地层滤失,扰动增产区域地层压力场和应力场[3-4],以剪切或张性破坏的形式激活周围处于闭合状态的天然裂缝[5-10],形成具有较高表观渗透率的复杂裂缝网络,极大地提高了页岩气井的产量。
天然裂缝或岩石在破坏的瞬间发生应力松弛,部分应变能以弹性波的形式释放,即产生微地震事件[11-13]。
国内非常规油气地球物理勘探技术现状及进展
国内非常规油气地球物理勘探技术现状及进展 水
赵万金 李海亮 杨午 阳
（ 中国石油勘探 开发研 究院西北分院 ，甘肃省兰州市 ７ ０ ２ ） ３ ０ ０
随 着 中国经济 持续 快速 的发 展 ，油气供需 矛 盾 已 经成 为制约 中国经济和社会发展 的主要瓶颈之一 。 ２ １ 年我 国原油 产量突破 ２×１ 。， ００ ０ｔ 但对 外依赖 度达到
１ 中国非常规油气资源分布 及开发现 状
非 常规 油 气指成藏机 理 、 赋存状态 、 布规律 及勘 分
第一作者 简介 ：赵万金 ，男 ，工程师 ，硕士 ， ０６ ２ ０ 年毕业于 中国石油大学 （ 华东 ） ，目前从事地 球物理勘探研 究工作。
收 稿 日期 ：２ １—１—０ ；修 改 日期 ：２ １—０ — ５ ０１ ２ ３ ０ ２ ５ ２
，
一 中国石油勘探 ２ １年第４ ３ ｒ ０２ 期 ６
开 发利 用起 来 ，对 于保 障 国家 能源安 全具 有重要 的战 探 明储量 ２ ８ Ｘ １ ｍ ． ０ 。我 国页 岩 气资 源量 达 １ Ｘ 略 意义乜 １ 。
１ ，其 中可 采储量 ２ ６Ｘ１ ，与美 国相 当 ；目 ０ｍ ． ０ ｍ。
文献 调研 发现 ，非 常规 油 气勘探 开发理 论 聚焦于 前 ，我 国页岩 气资 源仍处 于未 开发状 态 ，开发潜 力 巨 成 藏地 质理论 研究 和钻 井开发 技术研 究 ，对地球 物理 大 ； 国致密 砂岩 天然 气资源 量约 为 １ ４ ０ ， 我 ． ８Ｘ １ １ ３ ｍ
游离气）渗透率等参数 ； 、 ②烃源岩的生烃潜力评价 ， 主 ＡＶＡ裂 缝检测 技术 以及 多波 多分 量预 测技 术 。 要 包括 干酪 根的识 别 与类型 划分 、 有机 质含量 、 热成 熟 （） ４ 烃类检测技术： 主要体现在吸收和衰减两个方 度等一 系列指标 的定性或 定量解释 ； ③岩石 力学参数和 面 ， 并从 不 同角度有针 对性 地发展 了多项检测 技术 ， 例 裂 缝发 育指标 的评价 ［ 。 １ ４ 如 基于 流体 因子 检测技 术 、基于 神 经 网络 的 气层识 别
呼图壁地区超深致密砂岩多尺度裂缝地震预测技术
呼图壁地区超深致密砂岩多尺度裂缝地震预测技术作者：徐亚楠王晓涛习宇霄王贤叶迪苏艳丽来源：《新疆地质》2024年第01期摘要：针对呼图壁地区超深致密砂岩不同尺度裂缝预测问题，创新形成基于叠前反演和叠后多属性融合的多尺度裂缝预测方法。
据裂缝发育规模，将裂缝划分为大尺度裂缝和小尺度裂缝，其中大尺度裂缝采用叠后地震属性识别，小尺度裂缝成像测井约束反演识别；利用叠后地震数据提取进行机器深度学习、凌乱性检测与曲率增强进行大尺度裂缝检测，将三者结果进行融合以实现大尺度裂缝精细刻画；应用基于OVT域地震数据的叠前方位各向异性属性与叠前反演两种方法进行小尺度裂缝预测；综合大、小尺度裂缝预测结果分析目标区裂缝发育的规律，预测不同尺度裂缝发育有利区。
研究结果表明，在呼图壁地区超深储层应用多尺度裂缝地震预测方法获得的预测结果与实测钻孔结果吻合度较高。
关键词：致密砂岩储层；多尺度裂缝；呼图壁背斜；OVT域地震数据；裂缝预测裂缝是深层致密砂岩储层的重要储集空间和渗流通道，对准南呼图壁地区油气勘探开发具有重要的实际意义。
裂缝具明显的多尺度特征，不同尺度天然裂缝对致密低渗透储层的作用和对油气的影响不同，大尺度裂缝通常影响油气的保存；中小尺度裂缝控制了致密储层的渗流系统和储集作用，是影响致密砂岩储层高产及稳产的关键。
由于不同尺度裂缝的作用不同，需要分尺度研究天然裂缝的发育规律，才能更加科学有效地指导致密低渗透油气的勘探开发[1]。
呼图壁地区目的层埋深普遍超7 km，厚10～20 m，存在目的层埋深大、地震成像信噪比低、分辨率不足等问题，导致以往储层裂缝预测不能满足部署高效评价井的生产需要。
因此，提出对呼图壁地区砂岩储层不同尺度裂缝多技术手段预测方法，对多尺度裂缝的不同地震响应特征，分别采用叠后地震属性、叠前反演等技术手段进行预测，最后将多技术手段结果进行融合并优选[2-6]。
主要技术步骤可以概括为（图1）：①在叠后地震资料上先利用各种属性开展大尺度裂缝研究，如深度学习、凌乱性检测、曲率增强等属性均对大尺度裂缝响应较明显；②据地质露头、构造分析和钻录井等信息获得的构造裂缝的类型、产状、组系、密度等特征的认识，对叠后属性进行约束，将各类属性进行融合，优选得到大尺度裂缝预测结果；③用成像测井识别的裂缝来标定地震方位频率梯度属性及叠前方位各向异性反演，得到小尺度裂缝预测结果；④将以上两种尺度裂缝进行融合，最终得到储层裂缝体系的综合预测成果。
基于测井反演的CT扫描煤储层特征评价——以大宁-吉县区块为例
210中国深层煤层气资源丰富，位于鄂尔多斯盆地东南缘的大宁-吉县区块开展的深层煤层气勘探开发工作，证实深层煤层气勘探开发潜力巨大。
目前，国内外研究煤层储层特征的实验方法很多，常见方法有压汞、低温液氮吸附[1-2] 、CT扫描[3-4] 、X 射线计算机层析扫描[5]及核磁共振[6] 等。
本文主要采用测井与CT扫描技术相结合的方法，对大宁-吉县区块的煤岩储层特征分析，未压裂选层提供技术支撑。
1 煤层测井响应特征在测井曲线上，煤层的电性特征和放射性特征特别明显，阵列感应电阻率曲线（M2RX，M2R1，M2R2）、声波时差曲线（DT24）和补偿中子曲线（CNL）呈明显高幅值反映，密度曲线（DEN）和自然伽玛曲线（GR）呈明显的低幅值反映，自然电位曲线（SP）反映不明显。
灰岩顶板岩性致密，电阻率越高，密度越大。
灰岩的电阻率和密度呈高幅值，声波时差和自然伽玛曲线呈低幅值。
补偿中子视石灰岩孔隙度近似为零，自然电位幅值相对为零。
泥岩底板，泥岩层黏土能吸附水分子和放射性元素等微细小物质，所以泥岩的电阻率呈低幅值，自然伽玛曲线呈高幅值。
自然电位幅值相对为零（图1）。
图1 深层煤岩测井曲线特征图基于测井反演的CT扫描煤储层特征评价——以大宁-吉县区块为例武晓磊 李秀芳 周婷婷 陈明 孙俊义 余莉珠中石油煤层气有限责任公司临汾采气管理区 山西 临汾 041000 摘要：针对大宁-吉县区块煤岩发育特征，通过CT扫描与常规测井手段相结合，采用反演煤岩虚拟灰度剖面和岩心CT扫描煤层结构精细划分煤层技术，建立3条煤岩组分曲线：镜质组含量、壳质组+惰质组含量、矿物含量，划分夹矸和宏观煤岩类型，计算煤岩段有机质丰度和确定煤层结构类型。
煤岩灰度剖面定量出煤层段具体的宏观煤岩组分含量，有效实现单井煤层有机质丰度评价，结合薄夹矸和有机质丰度分布，可有效为压裂选段提供优化指导。
关键词：深8煤 CT扫描 灰度剖面 大宁-吉县区块Ｒｅｓｅｒｖｏｉｒ ｆｅａｔｕｒｅ ｅｖａｌｕａｔｉｏｎ ｂａｓｅｄ ｏｎ ｃｏｎｖｅｎｔｉｏｎａｌ ｍｅｔｈｏｄｓ ａｎｄ ｈｏｌｏｇｒａｐｈｉｃ ｄｉｇｉｔａｌ ｃｏｒｅ ｄｅｓｃｒｉｐｔｉｏｎ Ｔａｋｉｎｇ ｔｈｅ Ｄａｎｉｎｇ Ｊｉｘｉａｎ ｂｌｏｃｋ ａｓ ａｎ ｅｘａｍｐｌｅＷｕ Ｘｉａｏｌｅｉ，Ｌｉ Ｘｉｕｆａｎｇ， Ｚｈｏｕ Ｔｉｎｇｔｉｎｇ，Ｌｉ Ｈａｏｚｈｅｎｇ，Ｃｈｅｎ Ｍｉｎｇ，Ｓｕｎ Ｊｕｎｙｉ，Ｙｕ ＬｉｚｈｕPetroChina Coalbed Methane Co. ，Ltd. ，Linfen 041000Abstract ：Based on the development characteristics of coal and rock in the Daning Jixian block ，a combination of CT scanning and conventional logging techniques was used to invert the virtual gray scale profile of coal and rock and fine classify the coal seam structure using core CT scanning. Three coal and rock composition curves were established: vitrinite content ，crust+inertinite content ，mineral content ，classification of gangue and macroscopic coal and rock types ，calculation of organic matter abundance in coal and rock sections ，and determination of coal seam structure types. Quantitative determination of the specific macroscopic coal rock component content in coal seam sections using the gray scale profile of coal rock effectively achieves the evaluation of organic matter abundance in single well coal seams. Combined with the distribution of thin gangue and organic matter abundance ，it can effectively provide optimization guidance for fracturing section selection.Keywords ：Shen8 coal ；CT scanning ；Grayscale profile ；Daning Jixian Block2112 反演煤岩虚拟灰度剖面2.1 反演技术方法（1）宏观煤岩组分标定，全直径毫米扫描灰度数据与宏观煤岩组分拟合：煤岩组分、工业分析数据与全直径灰度值（G）标定，建立灰度值与煤岩组分定量关系：镜煤+亮煤，Mj=a*G+b 矿物，Mk=c*G+d （2）扫描煤岩井段测井数据与灰度（G）数据拟合：G=a*DEN+b*AC+c*GR+d，拟合后计算全井煤岩段虚拟灰度值（Gx）。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
基于球缝模型优化反演的裂缝储层评价新方法
刘诗琼;石强;汤继周;刘堂晏;赵文君;陈鹏
【期刊名称】《地球物理学报》
【年(卷),期】2022(65)4
【摘要】裂缝型储层评价是复杂油气田评价的难题之一,而且,针对裂缝型储层的勘探需求也在日益增加.本文将裂缝型储层中的裂缝归纳为裂缝与孔隙两个基本类型,然后,叠加裂缝与孔隙,在不同尺寸的配置关系下形成球缝模型.根据球缝模型提出孔隙结构参数的寻优方法,结合现有的反演算法,形成针对核磁测井回波的优化反演算法,基于优化反演结果可获取裂缝均值宽度、孔隙半径、分选系数、峰态和其他孔隙结构参数.对比分析孔隙结构参数与储层品质因子及油气产量之间的关系,可以发现孔隙结构参数与储层品质因子具有很好的相关性,与储层测试产量也有良好的相关性.对比结果说明,基于球缝模型反演的储层孔隙结构参数,对于裂缝型储层的物性具有很好的刻画效果.因此,本文的研究成果对于提高裂缝型储层的评价质量具有潜在的应用前景.
【总页数】10页(P1451-1460)
【作者】刘诗琼;石强;汤继周;刘堂晏;赵文君;陈鹏
【作者单位】西南石油大学油气藏地质及开发工程国家重点实验室;同济大学海洋地质国家重点实验室;中国石油勘探开发研究院
【正文语种】中文
【中图分类】P631
【相关文献】
1.页岩储层多缝应力干扰下裂缝起裂角计算新模型
2.一种基于多点地质统计的储层反演新方法
3.基于粒子群优化算法的小波神经网络缝洞型储层识别模型
4.碳酸盐岩缝洞型储层有效性评价新方法
因版权原因，仅展示原文概要，查看原文内容请购买。