地质力学学报  2020, Vol. 26 Issue (1): 74-83
引用本文
任浩林, 刘成林, 刘文平, 杨熙雅, 李文研. 四川盆地富顺-永川地区五峰组—龙马溪组应力场模拟及裂缝发育区预测[J]. 地质力学学报, 2020, 26(1): 74-83.
REN Haolin, LIU Chenglin, LIU Wenping, YANG Xiya, LI Wenyan. Stress field simulation and fracture development prediction of the Wufeng Formation-Longmaxi Formation in the Fushun-Yongchuan Block, Sichuan Basin[J]. Journal of Geomechanics, 2020, 26(1): 74-83.
四川盆地富顺-永川地区五峰组—龙马溪组应力场模拟及裂缝发育区预测
任浩林1,2, 刘成林1,2, 刘文平3, 杨熙雅1,2, 李文研4    
1. 中国石油大学(北京)地球科学学院, 北京 102249;
2. 中国石油大学(北京)油气资源与探测国家重点实验室, 北京 102249;
3. 中国石油西南油气田分公司页岩气研究院, 四川 成都 610051;
4. 中国石油青海油田勘探事业部, 甘肃 敦煌 736202
摘要:岩石中构造裂缝主要受控于地层所处的区域构造应力场,地应力对油气的运移、成藏和分布有着重要的作用。通过应用ANSYS有限元数值模拟方法对四川盆地富顺-永川区块的五峰组-龙马溪组进行应力场模拟及分析研究、应用构造曲率法中的三点法对该层段页岩进行张裂缝发育情况预测研究、综合模拟结果与曲率计算数据对目的区域目的层段进行裂缝发育强度综合预测研究,结果表明富顺-永川区块五峰组-龙马溪组应力高值沿背斜走向分布的规律明显,研究区的东部及西南部的背斜应力值较其他区域背斜高;曲率高值沿背斜走向分布的规律明显,研究区东部及西南部背斜曲率值较其他区域背斜高;在现今应力的作用下,研究区背斜处裂缝发育程度较高,向斜处裂缝发育程度较低,研究区东部及西南部背斜裂缝发育程度较其他区域高。
关键词曲率法    有限元    裂缝预测    龙马溪组    四川盆地    
DOI10.12090/j.issn.1006-6616.2020.26.01.008     文章编号:1006-6616(2020)01-0074-10
Stress field simulation and fracture development prediction of the Wufeng Formation-Longmaxi Formation in the Fushun-Yongchuan Block, Sichuan Basin
REN Haolin1,2, LIU Chenglin1,2, LIU Wenping3, YANG Xiya1,2, LI Wenyan4    
1. College of Geosciences, China University of Petroleum, Beijing 102249, China;
2. State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum, Beijing 102249, China;
3. Shale Gas Research Institute, PetroChina Southwest Oil & Gas Field Company, Chengdu 610051, Sichuan, China;
4. Petroleum Exploration Division, Petro China Qinghai Oilfield, Dunhuang 736202, Gansu, China
Abstract: The tectonic fractures in rocks are mainly controlled by the regional tectonic stress field where the strata are located. Tectonic stress plays an important role in the migration, accumulation and distribution of oil and gas. The stress field of the Wufeng Formation-Longmaxi Formation in the Fushun-Yongchuan Block of the Sichuan Basin was simulated and analyzed by ANSYS finite element software; The three-point method in the method of structural curvature method was used to predict the development of cracks in the shale, and the comprehensive simulation results and curvature calculation data were used to predict the fracture intensity in the target area. The results proved that the distribution of the high stress value of the Wufeng Formation-Longmaxi Formation in the Fushun-Yongchuan Block is along the anticline, and the anticline stress values in the eastern and southwestern part of the study area are higher than those in other areas; the distribution of high curvature is along the anticline, and the curvature values of the anticline in the east and southwest of the study area are higher than those in other areas; the fracture in the anticline of the study area has a higher level of development, and it's lower in the syncline, at the same time it's higher in the anticline of the eastern and southwestern part of the study area than other areas under the current stress.
Key words: curvature method    finite element    fracture prediction    Longmaxi Formation    Sichuan Basin    
0 引言

勘探现状表明,在有利的保存单元里页岩气的甜点受裂缝的发育程度严格控制,因此对裂缝系统的研究是油气勘探开发过程中重点关注的方面(Finkbeiner et al., 1997宋岩等,2013邹才能等,2015刘树根等,2016颜彩娜等,2016董大忠等,2018)。裂缝类型多而复杂,裂缝的发育主要受外因(构造应力场)和内因(岩相岩性、层厚等)的联合控制(Du Rouchet,1981张守仁等,2004Nelson,2011黄生旺等,2017张学敏等,2018)。因此,目前对于裂缝的发育情况以及发育特征大多是从力学成因角度,通过对构造应力场使用不同的方法进行模拟来进行评价:部分研究者在研究过程中选择以建立地质模型为渠道(丁文龙等,2016),通过对岩相的划分、不同岩相岩石力学性质参数进行分析,将地质模型转化为岩石力学模型来对应力场进行模拟;部分研究者选择通过数学公式推导(张守仁等,2004),采用对岩石应力-应变曲线进行数学处理、划分微小单元、在单元内采用有限元法逼近应变曲线的方法,对目标地层的弹性问题进行计算,从而采集地层的应力场数据并作出模拟;部分研究者采用数字化的岩石分析手段对裂缝的发育程度进行准确计算(Ramandi et al., 2017);更多的研究者选择使用有限元软件对岩层进行受力分析,计算地层中地应力的分布(Al-Ruwaili and Chardac, 2003蒋有录等,2005魏春光等,2006刘广峰等,2009Petricca et al. 2013吴林强等,2014Mohtarami et al., 2017)。有限元方法相比其他方法的优势在于将传统的地质学问题和地质建模问题转化为数学与计算方面的求解问题,即将地层视为材料整体,将各类岩石视为不同种类的材料,将大区域、大面积的地质应力分析转化为对材料的受力分析,一定程度上简化求取和模拟应力场的步骤。

文章以四川盆地富顺-永川地区前期工程资料为基础,对岩石力学参数数据进行了分析,采用古应力场有限元数值模拟方法预测了剪裂缝的分布特征,同时采用构造曲率法预测了张裂缝的发育程度,最后通过归一化处理,从而对研究区块五峰组—龙马溪组进行裂缝发育综合强度的预测,为后期气藏开发选区提供数据支撑。

1 区域地质背景

富顺-永川区块位于四川省自贡市富顺县至重庆市永川区之间(图 1),面积约4300 km2,区域构造位置位于川南地区永川-宜宾双行雁列褶皱束上。乐光禹等(1996)认为,川南地区构造边界由三条主要断裂带交错衔接(东界:长寿-遵义断裂带;南界:古蔺-盐津断裂带;北西界:华蓥山-青山岭断裂带)并围限成直角三角形,三角形三边分别受力,造就了川南地区独特的构造形态,永川-宜宾双行雁列褶皱束的构造形态主要受直角三角形块体的斜边控制,其构成主要构造形迹的最大主应力方向为北西—南东向。永川-宜宾双行雁列褶皱束的构造形态以一系列大型的长轴紧密背斜夹宽缓向斜的褶皱构造特征为主,研究区位于该褶皱束的中上部,由于研究区面积仅占川南构造带面积一小部分,且远离构造带边界,因此在构造单元划分上处于最低的三级构造单元。对其构造形态的分析也较为简单,主要分析研究区的背斜向斜特征及组合关系。研究区内的几个背斜走向相差较小,背斜走向大多为北东—南西向,紧闭的背斜与宽缓的向斜形成侏罗山式的箱型褶皱(图 2),背斜的枢纽长度比较长,大部分斜向穿过整个研究区。

图 1 富顺-永川区块地理位置图 Fig. 1 The location of the Fushun-Yongchuan Block

图 2 富顺-永川区五峰组底面三维构造图 Fig. 2 The 3-D structural picture of the bottom of the Wufeng Formation in the Fushun-Yongchuan Block

富顺-永川区块目的层五峰组—龙马溪组在研究区内广泛发育,该区五峰组—龙马溪组沉积相为水体宁静、还原程度较高的深水陆棚亚相,主要沉积黑色碳质页岩、硅质页岩或含粉砂页岩、夹薄板状硅质岩和碳质泥页岩(郭英海等,2004牟传龙等,2014王同等,2015),该亚相黑色页岩分布稳定,沉积厚度大,有机质丰度高,是四川盆地页岩气勘探开发最主要的目的层段(图 3)(刘文平等,2017)。受沉积环境和古水深的控制,不同地区和不同沉积期的矿物组成特征存在较大差异,并与有机质的富集关系密切。

图 3 四川盆地古生界(部分)地层综合柱状图(刘文平等,2017) Fig. 3 Comprehensive histogram of the Paleozoic strata (partial) in the Sichuan Basin (Liu et al., 2017)

该区断层发育较为复杂,狭窄的背斜脊部地层倾角大,复杂的地质构造导致裂缝分布情况不明确,急需开展应力场模拟及裂缝预测的研究,为后续工作及其他研究区的勘探开发提供技术支撑。

2 有限元数值模拟法预测剪裂缝发育程度

研究选择以现有资料为基础,以有限元数值模拟方法为手段,通过应力场的数值模拟对裂缝主要形成时期的古应力场进行了预测,结合破裂准则评价了研究区的剪裂缝发育情况。

2.1 地质模型及力学模型的建立

(1) 以矿权界线为边界,模型底面通过五峰组底面的数字化构造等高线图联立Surfer软件共同建立(刘涛和杨凤鹏,2002Moaveni et al., 2003崔芳鹏等,2008赵维涛和陈孝珍,2009),顶面及四个侧面通过布尔操作生成;由于三维地质建模的复杂性,且目的层厚度(350~500 m)远小于研究区的长宽(长97.5 km,宽40.8 km),因此将目的层作等厚、均质处理,厚度(即模型高度)取各井目的层厚度平均值457 m;(2)模型的力学参数对模拟结果有重要的影响(陈运平等,2008),因此模型力学参数的设定需具有代表性。力学参数(杨氏模量、泊松比)的数据来源为大量测井数据解释而来的动态数据(图 4),泊松比值大部分落在0.12~0.24区间内、平均值0.17,杨氏模量值大部分落在2×104~4×104 MPa区间内、平均值3.0×104 MPa,力学参数与深度变化间关系不明显,且其分布范围较为集中。为弱化特殊的高值或低值点在参数设定方面带来的影响,杨氏模量设定为中位数值2.84×104 MPa,泊松比设定为中位数值0.17。

a—泊松比—深度关系图;b—杨氏模量—深度关系图 图 4 富顺-永川区块五峰组—龙马溪组力学参数分布 Fig. 4 Distribution of rock's mechanical parameters in the Wufeng Formation—Longmaxi Formation, Fushun-Yongchuan Block

(3) 断层的处理方式沿用已有研究中的处理方法:假定为断裂带,采用断裂带方式处理,适当降低断层两侧适当距离岩性的杨氏模量(降低至60%),同时增加泊松比(增加0.02)(周灿灿,2003刘广峰等,2009宋伟等,2017)。

(4) 已建好的模型采用Solid20节点实体,通过软件网格化功能自动划分单元,共建立51084个三维立体单元,283239个节点,平均每平方千米65个节点,满足了常规要求的每平方千米3~5个输出点(穆龙新,2009)。

2.2 模型受力设定

模型所受应力主要包括垂向的上覆岩层重力与水平的应力。目的层垂向受力大小与深度的拟合关系良好(图 5a),因此采用该拟合关系所得函数来计算模型中节点的垂向力大小,软件中自动计算并代入相关运算;使用前期勘探过程中测定的研究区及周边的燕山期主幕水平应力(σ)大小(60~90 MPa)共12个数据(表 1)进行直线拟合,通过拟合直线计算目的层所受水平应力大小(图 5b),经计算目的层不同深度处受水平主应力大小的范围为80~110 MPa,模型载荷加载处(对应埋深3660 m处地层)受力约70 MPa。

a—垂向上覆岩层压力-深度关系;b—水平应力-深度关系 图 5 富顺-永川区块五峰组—龙马溪组模型所受应力值拟合图 Fig. 5 Stress value fitting of the Wufeng Formation—Longmaxi Formation model in the Fushun-Yongchuan Block

表 1 四川盆地燕山期主幕水平主应力大小计算表 Table 1 Computation chart of horizontal stress of Yenshanian main curtain in the Sichuan Basin

乐光禹等(1996)通过对川南地区构造形态模拟实验的分析认为,川南地区直角三角形构造在形成时期所受地应力为垂直斜边的挤压力以及二直角边处的反作用力(支撑力),而研究区所处的永川-宜宾褶皱束的构造形态受上述的挤压力控制。结合模型自身形状特点,该挤压力在研究区处被分解为自北向南、自西向东两个分力。

综上,模型受力情况最终确定为(图 6):自上而下的重力挤压(载荷)与下伏岩层提供的支撑(约束);模型北侧面、西侧面受垂直模型面的挤压(载荷,由水平主应力分解而来)与南侧面、东侧面的支撑(约束)。图 6中的模型为研究构建的研究区目的层地质模型,其顶底面形状、厚度、材料参数等固有属性按上文所述进行设置,模型中的深色有限元为上文所述采用断裂带方式处理的含断层的龙马溪组页岩。

图 6 富顺-永川区块五峰组—龙马溪组模型载荷与约束示意图 Fig. 6 Load and constraint of the Wufeng Formation—Longmaxi Formation in the Fushun-Yongchuan Block
2.3 模拟结果

进行有限元模拟后,获得研究区目的层燕山期差应力分布图(图 7)及最大主应力分布图(图 8)。图中压力数值的正值表示挤压力,负值表示拉张力。由于构造应力值的高、低值区一般分别对应裂缝的发育、不发育区(沈国华,2008),因此可通过某区域地应力的大小间接判断裂缝的发育程度,并可通过地应力的性质(挤压力或拉张力)对可能存在的裂缝的性质进行预测。

图 7 富顺-永川区块五峰组—龙马溪组差应力分布图 Fig. 7 Distribution of stress intensity in the Wufeng Formation—Longmaxi Formation, Fushun-Yongchuan Block

图 8 富顺-永川区块五峰组—龙马溪组最大主应力分布图 Fig. 8 Distribution of maximum principal stress in the Wufeng Formation—Longmaxi Formation, Fushun-Yongchuan Block

图 7分析得出,研究区目的层差应力分布范围为13~271 MPa,其范围较大;差应力受构造的控制作用明显,紧闭背斜的差应力(120~271 MPa)均比宽缓向斜(40~70 MPa)大,高值区沿断层和背斜呈条带状分布的现象明显;东部背斜的差应力(150~271 MPa)较西部(150 MPa左右)更大,东部背斜差应力达到最大值(271 MPa);目前大部分井位于差应力低值区,而井03、井04、井14等井位于差应力相对高值区。

图 8分析得出,研究区目的层最大主应力在背斜处多为挤压力,而在宽缓向斜处则均为拉张力;东部及西部背斜的最大主应力(41~93 MPa)较中部(15~41 MPa)大,高值区沿断层和背斜呈条带状分布的现象明显,东部背斜最大主应力达到最大值(170 MPa);目前大部分井位于最大主应力低值区,而井03、井04、井14、井16等井位于最大主应力相对高值区。由该图分析,研究区目的层的背斜应发育挤压成因缝,向斜应发育拉张成因缝,裂缝发育区应沿断层呈条带状分布,且东部及西部的背斜较研究区其他背斜有更发育的裂缝。

3 构造曲率法预测张裂缝发育程度 3.1 方法实现

构造曲率法是通过计算构造面主曲率,综合区域不同位置的曲率大小进行裂缝发育区预测的一种方法。由于地层的曲率值在一定程度上反映岩层弯曲的程度,因此可以利用地层某处的曲率值来评价该处因构造应力产生的张裂缝的发育情况(郭科等,1998)。

曲率的计算及裂缝发育区预测选用曲率法中较为简洁的三点法进行(曹润荣和刘宗彦,2008):在网格化等高线数据中,垂直构造走向的方向上依次提取相邻三点的海拔高度,该三点坐标及海拔高度即可确定一个圆,根据坐标及海拔高度即可计算圆心坐标,进而求出三点的平均曲率值K

$ p = a + b + c $ (1)
$ K = \frac{{\sqrt {p \times (p - 2a) \times (p - 2b) \times (p - 2c)} }}{{abc}} $ (2)

式中abc为三角形三边(即每组每两点点距),K为经过该组三点的平均曲率。据此,求得全区共437122个网格点的曲率值,并将该值以等值线形式绘制成图。

3.2 裂缝预测结果

经计算及绘制等值线图后,叠合区域断层分布图可得到研究区的曲率数值分布图(图 9)。

图 9 富顺-永川区块五峰组—龙马溪组曲率分布-断层位置叠合图 Fig. 9 Overlay of curvature distribution and fault position in the Wufeng Formation—Longmaxi Formation, Fushun-Yongchuan Block

图 8分析得出,富顺-永川地区五峰组—龙马溪组曲率值分布范围为0~1.6×10-3,紧闭背斜的曲率值(0.4×10-3~1.6×10-3)较宽缓向斜(0~0.4×10-3)高,曲率高值沿断层和背斜展布的特征明显;东部及西部背斜的曲率(0.8×10-3~1.6×10-3)较中部(0.4×10-3~0.8×10-3)大,东部背斜曲率达到最大值(1.6×10-3);目前大部分井位于曲率低值区,而井03、井4、井14等井位于曲率相对高值区。曲率值的分布特征与应力场模拟结果中应力的分布特征十分相似:二者高值均沿断层展布、研究区东部及西部背斜的二者数值较中部大、二者均在东部背斜处达到最大值。

4 裂缝发育程度的综合预测

裂缝有不同的力学成因类型,其中以挤压力产生的剪切缝和拉张力产生的张性缝为主要的裂缝成因。已有调查表明,裂缝的发育程度是不同成因类型的裂缝综合叠加的结果。因此,在实际评价过程中,往往采用归一化处理方法,综合不同成因裂缝的发育程度来对整体的裂缝发育程度进行综合预测。在上述应力场模拟及曲率计算的基础上,选取区域内应力值(曲率值)最大的点作为标尺,计算各点应力值(曲率值)与应力值(曲率值)最大点的相对大小(0~1),最终将两种值的相对大小暂按权重1:1相加(该权重有待与后期工程测试等提供的现场资料结合并进行修改),获得各点裂缝强度(即研究区内各点裂缝发育的程度)的数值,从而获得初步的裂缝强度的分布图(图 10)。

图 10 富顺-永川区块五峰组—龙马溪组裂缝强度分布图 Fig. 10 Distribution of crack strength in the Wufeng Formation—Longmaxi Formation, Fushun-Yongchuan Block

图 10分析得出,富顺-永川地区五峰组—龙马溪组裂缝强度分布范围为0~0.5,紧闭背斜的裂缝强度(0.05~0.5)较宽缓向斜(0~0.05)高,裂缝强度高值沿断层和背斜展布的特征明显;东部及西部背斜的裂缝强度(0.15~0.5)较中部(0.05~0.15)大,东部背斜裂缝强度达到最大值(0.5);目前大部分井位于裂缝强度低值区,而井03、井04、井12、井14等井位于裂缝强度相对高值区。结合图 2中的区域形态构造特征和图 10中的裂缝强度发育特征进行分析可得出,裂缝强度的高值区与褶皱发育区重合程度较大,裂缝强度的高值点均位于褶皱的枢纽上,与此同时断层的走向也与褶皱的走向一致、且平行于条带状的裂缝强度高值区,说明了裂缝发育程度受构造部位的控制十分明显。同时,研究区褶皱发育、曲率高值与剪切力高值集中于背斜处的特征明显,研究区背斜处张剪裂缝均有可能大量发育。

5 结论

(1) 富顺-永川区块五峰组—龙马溪组的背斜具有较高的差应力值、挤压应力值,其中位于富顺-永川区块东部及西南部的背斜在横向对比中可发现具有全区最高的应力值,根据工程中获得的经验,差应力的高值区在钻井时可能出现井壁坍塌和卡钻现象,完井后可能出现套管损坏,而低差应力的区域在进行压裂操作时难以控制水力裂缝的走向,因此后期施工可以依据文中模拟结果对钻井和开发方案进行调整。

(2) 富顺-永川区块五峰组—龙马溪组的背斜具有较高的曲率值,曲率大小的分布规律与应力场模拟中各应力值的分布规律一致,位于富顺-永川区块东部及西南部的五峰组—龙马溪组背斜在横向对比中可发现具有更发育的张裂缝,由于龙马溪组为海相地层,大多数张裂缝出现不同程度的方解石充填现象,其渗透率受充填情况的影响较大,后期开发需考虑方解石充填带来的储层非均质变化和敏感程度变化。

(3) 裂缝发育程度受构造部位的控制十分明显,研究区曲率高值与剪切力高值集中于背斜枢纽处的特征明显,研究区背斜枢纽处张剪裂缝均有可能大量发育,综合来看地层褶皱枢纽是裂缝发育程度最高的区域,该区域的渗透率可能较高,连通性较好,但过于发育的裂缝网络可能会对钻井、完井和后期的压裂及开发带来一定影响,下一步的设计方案需有针对性的考虑。

参考文献/References
AL-RUWAILI S B, CHARDAC O, 2003. 3D model for rock strength & in-situ stresses in the Khuff formation of Ghawar field, methodologies & applications[C]//Middle east oil show. Bahrain: Society of Petroleum Engineers: 206-219.
CAO R R, LIU Z Y. Application of curvature of face-trend surface fitting method in fracture prediction[J]. Computer Applications of Petroleum, 2008, 16(3): 12-14. (in Chinese with English abstract)
CHEN Y P, ZHAO C B, LIN G. Mechanical properties of deep earth rocks and their roles in the investigation of continental deformation processes[J]. Geotectonica et Metallogenia, 2008, 32(3): 276-284. (in Chinese with English abstract)
CUI F P, HU R L, LIU Z L, et al. Surfer software platform based complex three-dimensional geological digital models for pre-processing of FLAC3D[J]. Journal of Engineering Geology, 2008, 16(5): 699-702. (in Chinese with English abstract)
DING W L, ZENG W T, WANG R Y, et al. Method and application of tectonic stress field simulation and fracture distribution prediction in shale reservoir[J]. Earth Science Frontiers, 2016, 23(2): 63-74. (in Chinese with English abstract)
DONG D Z, SHI Z S, SUN S S, et al. Factors controlling microfractures in black shale:a case study of Ordovician Wufeng Formation-Silurian Longmaxi Formation in Shuanghe Profile, Changning area, Sichuan Basin, SW China[J]. Petroleum Exploration and Development, 2018, 45(5): 763-774. (in Chinese with English abstract)
DU ROUCHET J. Stress fields, a key to oil migration[J]. AAPG Bulletin, 1981, 65(1): 74-85.
FINKBEINER T, BARTON C A, ZOBACK M D. Relationships among in-situ stress, fractures and faults, and fluid flow:Monterey formation, Santa Maria basin, California[J]. AAPG Bulletin, 1997, 81(12): 1975-1999.
GUO K, XU Z Y, NI G S, et al. Research on the main curvature method and its application to cracky oil-gas deposits[J]. Computing Techniques for Geophysical and Geochemical Exploration, 1998, 20(4): 335-337. (in Chinese with English abstract)
GUO Y H, LI Z F, LI D H, et al. Lithofacies palaeogeography of the Early Silurian in Sichuan area[J]. Journal of Palaeogeography, 2004, 6(1): 20-29. (in Chinese with English abstract)
HUANG S W, JIANG M L. Analysis on structural fracture characteristics in Ansai Yanhewan Area[J]. Liaoning Chemical Industry, 2017, 46(10): 973-975. (in Chinese with English abstract)
JIANG Y L, ZHANG L, LU X S, et al. Application of the tectonic stress simulation based on ANSYS in Kelasu region of Kuche depression[J]. Natural Gas Industry, 2005, 25(4): 42-44. (in Chinese with English abstract)
LIU G F, LU H J, HE S L. Application of finite element analysis in reservoir in-situ stress research[J]. Science Technology and Engineering, 2009, 9(24): 7430-7435. (in Chinese with English abstract)
LIU S G, DENG B, ZHONG Y, et al. Unique geological features of burial and superimposition of the lower Paleozoic shale gas across the Sichuan basin and its periphery[J]. Earth Science Frontiers, 2016, 23(1): 11-28. (in Chinese with English abstract)
LIU T, YANG F P. Proficient in ANSYS[M]. Beijing: Tsinghua University Press, 2002. (in Chinese)
LIU W P, ZHANG C L, GAO G D, et al. Controlling factors and evolution laws of shale porosity in Longmaxi Formation, Sichuan Basin[J]. Acta Petrolei Sinica, 2017, 38(2): 175-184. (in Chinese with English abstract)
MOAVENI S, 2003. Finite element analysis: theory and application with ANSYS[M]. 2nd ed. Publishing House of Electronics Industry.
MOHTARAMI E, BAGHBANAN A, HASHEMOLHOSSEINI H. Prediction of fracture trajectory in anisotropic rocks using modified maximum tangential stress criterion[J]. Computers and Geotechnics, 2017, 92: 108-120.
MOU C L, GE X Y, XU X S, et al. Lithofacies palaeogeography of the Late Ordovician and its petroleum geological significance in Middle-Upper Yangtze Region[J]. Journal of Palaeogeography, 2014, 16(4): 427-440. (in Chinese with English abstract)
MU L X. Prediction of reservoir fractures[M]. 2009. (in Chinese)
NELSON P H. Pore-throat sizes in sandstones, siltstones, and shales:reply[J]. AAPG Bulletin, 2011, 95(8): 1448-1453.
PETRICCA P, CARAFA M M C, BARBA S, et al. Local, regional, and plate scale sources for the stress field in the Adriatic and Periadriatic region[J]. Marine and Petroleum Geology, 2013, 42: 160-181.
RAMANDI H L, MOSTAGHIMI P, ARMSTRONG R T. Digital rock analysis for accurate prediction of fractured media permeability[J]. Journal of Hydrology, 2017, 554: 817-826.
SHEN G H. Application of the finite element numerical simulation method in fracture prediction[J]. Petroleum Geology and Recovery Efficiency, 2008, 15(4): 24-26, 29. (in Chinese with English abstract)
SONG W, MA X J, LV F L, et al. The two-dimensional plane model study of tectonic ground fissure stress field in Hejian County of Cangzhou city, Hebei Province[J]. Shanghai Land & Resources, 2017, 38(3): 83-89. (in Chinese with English abstract)
SONG Y, JIANG L, MA X Z. Formation and distribution characteristics of unconventional oil and gas reservoirs[J]. Journal of Palaeogeography, 2013, 15(5): 605-614. (in Chinese with English abstract)
WANG T, YANG K M, XIONG L, et al. Shale sequence stratigraphy of Wufeng-Longmaxi Formation in southern Sichuan and their control on reservoirs[J]. Acta Petrolei Sinica, 2015, 36(8): 915-925. (in Chinese with English abstract)
WEI C G, LEI M S, WAN T F, et al. Numerical simulation of palaeotectonic stress field of Yingcheng Fm in Gulong-Xujiaweizi area:prediction and comparative study of tectoclase development area[J]. Oil & Gas Geology, 2006, 27(1): 78-84, 105. (in Chinese with English abstract)
WU L Q, LIU C L, LI B, et al. Numerical simulation of tectonic stress field and prediction of oil-favored areas:a case study of the third member of Qingshankou Formation in Guizijing region of Qian'an area, Songliao Basin[J]. Journal of Geomechanics, 2014, 20(4): 339-351. (in Chinese with English abstract)
YAN C N, JIN Z J, ZHAO J H, et al. Comparison of Marcellus shale in United States and Longmaxi formation shale in southern China[J]. Geological Science and Technology Information, 2016, 35(6): 122-130. (in Chinese with English abstract)
YUE G Y, DU S Q, HUANG J J, et al. Principle of structural compounding-combine[M]. Chengdu: The Chengdu University of Science and Technology Press, 1996. (in Chinese)
ZHANG S R, WAN T F, CHEN J P. Tectonic stress field modeling and fracture prediction in strata in Xiaoquan-Xinchang area, western Sichuan depression[J]. Oil & Gas Geology, 2004, 25(1): 70-74, 80. (in Chinese with English abstract)
ZHANG X M, YIN S, SHI C L. Developmental characteristics and controlling factors of fractures in tight sandstone of Shanxi Formation, southern Qinshui Basin[J]. Marine Origin Petroleum Geology, 2018, 23(3): 43-52. (in Chinese with English abstract)
ZHAO W T, CHEN X Z. Foundation of finite element method[M]. Beijing: Science Press, 2009. (in Chinese)
ZHOU C C, 2003. Studies on the structure mode of Baigezhuang region and the identification and prediction of structure fracture of reservoirs[D]. Guangzhou: Guangzhou Institute of Geochemistry, Chinese Academy of Sciences3. (in Chinese with English abstract)
ZOU C N, DONG D Z, WANG Y M, et al. Shale gas in China:characteristics, challenges and prospects (I)[J]. Petroleum Exploration and Development, 2015, 42(6): 689-701. (in Chinese with English abstract)
曹润荣, 刘宗彦. 面曲率法-趋势面拟合法在裂缝预测中的应用[J]. 石油工业计算机应用, 2008, 16(3): 12-14.
陈运平, 赵崇斌, 林舸. 深部岩石力学性质及其在大陆构造变形过程研究中的作用[J]. 大地构造与成矿学, 2008, 32(3): 276-284.
崔芳鹏, 胡瑞林, 刘照连, 等. 基于Surfer平台的FLAC3D复杂三维地质建模研究[J]. 工程地质学报, 2008, 16(5): 699-702.
丁文龙, 曾维特, 王濡岳, 等. 页岩储层构造应力场模拟与裂缝分布预测方法及应用[J]. 地学前缘, 2016, 23(2): 63-74.
董大忠, 施振生, 孙莎莎, 等. 黑色页岩微裂缝发育控制因素:以长宁双河剖面五峰组-龙马溪组为例[J]. 石油勘探与开发, 2018, 45(5): 763-774.
郭科, 胥泽银, 倪根生. 用主曲率法研究裂缝性油气藏[J]. 物探化探计算技术, 1998, 20(4): 335-337.
郭英海, 李壮福, 李大华, 等. 四川地区早志留世岩相古地理[J]. 古地理学报, 2004, 6(1): 20-29.
黄生旺, 姜萌蕾. 安塞沿河湾地区构造裂缝特征分析[J]. 辽宁化工, 2017, 46(10): 973-975.
蒋有录, 张乐, 鲁雪松, 等. 基于ANSYS的应力场模拟在库车坳陷克拉苏地区的初步应用[J]. 天然气工业, 2005, 25(4): 42-44.
刘广峰, 陆红军, 何顺利. 有限元法开展油气储层地应力研究综述[J]. 科学技术与工程, 2009, 9(24): 7430-7435.
刘树根, 邓宾, 钟勇, 等. 四川盆地及周缘下古生界页岩气深埋藏-强改造独特地质作用[J]. 地学前缘, 2016, 23(1): 11-28.
刘涛, 杨凤鹏. 精通ANSYS[M]. 北京: 清华大学出版社, 2002.
刘文平, 张成林, 高贵冬, 等. 四川盆地龙马溪组页岩孔隙度控制因素及演化规律[J]. 石油学报, 2017, 38(2): 175-184.
牟传龙, 葛祥英, 许效松, 等. 中上扬子地区晚奥陶世岩相古地理及其油气地质意义[J]. 古地理学报, 2014, 16(4): 427-440.
穆龙新. 储层裂缝预测研究[M]. 北京: 石油工业出版社, 2009.
沈国华. 有限元数值模拟方法在构造裂缝预测中的应用[J]. 油气地质与采收率, 2008, 15(4): 24-26, 29.
宋伟, 马学军, 吕凤兰, 等. 沧州河间构造地裂缝应力场二维平面模型研究[J]. 上海国土资源, 2017, 38(3): 83-89.
宋岩, 姜林, 马行陟. 非常规油气藏的形成及其分布特征[J]. 古地理学报, 2013, 15(5): 605-614.
王同, 杨克明, 熊亮, 等. 川南地区五峰组-龙马溪组页岩层序地层及其对储层的控制[J]. 石油学报, 2015, 36(8): 915-925.
魏春光, 雷茂盛, 万天丰, 等. 古龙-徐家围子地区营城组古构造应力场数值模拟-构造裂缝发育区带预测及对比研究[J]. 石油与天然气地质, 2006, 27(1): 78-84, 105.
吴林强, 刘成林, 李冰, 等. 应力场数值模拟与油藏有利区预测:以松辽盆地乾安地区归字井青三段为例[J]. 地质力学学报, 2014, 20(4): 339-351.
颜彩娜, 金之钧, 赵建华, 等. 美国Marcellus页岩与中国南方龙马溪组页岩地质特征对比及启示[J]. 地质科技情报, 2016, 35(6): 122-130.
乐光禹, 杜思清, 黄继钧, 等. 构造复合联合原理:川黔构造组合叠加分析[M]. 成都: 成都科技大学出版社, 1996.
张守仁, 万天丰, 陈建平. 川西坳陷孝泉-新场地区须家河组二-四段构造应力场模拟及裂缝发育区带预测[J]. 石油与天然气地质, 2004, 25(1): 70-74, 80.
张学敏, 尹帅, 史长林. 沁水盆地南部山西组致密砂岩裂缝发育特征及控制因素[J]. 海相油气地质, 2018, 23(3): 43-52.
赵维涛, 陈孝珍. 有限元法基础[M]. 北京: 科学出版社, 2009.
周灿灿, 2003.柏各庄地区构造样式及储层构造裂缝识别与预测[D].广州: 中国科学院研究生院(广州地球化学研究所). http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=Y540848
邹才能, 董大忠, 王玉满, 等. 中国页岩气特征、挑战及前景(一)[J]. 石油勘探与开发, 2015, 42(6): 689-701.