2. 中国石化石油勘探开发研究院 北京 100083;
3. 中海油深圳分公司研究院 广州 510240
勘探开发实践表明,裂缝性油藏是我国含油气盆地中重要的油气藏类型,合理高效地开发此类油藏对我国石油工业的持续发展及能源战略安全具有重要意义(袁士义等,2004; 李新景等,2007; 曾联波等,2010; 王珂等,2015)。由于油藏中普遍发育复杂的裂缝系统,裂缝不但影响着油气的储集和运移,而且对不同阶段的生产开发也产生着重要影响(梁官忠等,2001; 杨永林等,2004; 曾联波等,2007; 陈淑利等,2008; 李跃纲等,2012; 赵向原等,2015),科学评价裂缝发育特征,明确其形成机理、 控制因素及分布规律,可为有效开发裂缝性油气藏提供地质依据。渤海海域Z油田沙一下生物灰岩油藏为典型的裂缝性油藏,油藏地质条件复杂,注采矛盾突出。前人并没有对该区天然裂缝展开系统研究,对影响开发效果的裂缝类型及其特征认识不清。为此,本文综合利用岩心、 成像测井及分析测试等资料,对裂缝的发育特征进行了研究,并在研究区构造演化特征的基础上,对裂缝的成因机理进行了分析。
1 地质背景渤海湾地区是我国华北含油气盆地的重要组成部分,位于华北克拉通东北缘,属于新生代裂谷盆地(杨明慧,2009; 滕长宇等,2014; 陈飞等,2015; 李理等,2015)。该地区囊括渤中坳陷、 黄骅坳陷、 济阳坳陷和下辽河坳陷海域部分等多个渤海湾盆地次级构造单元带以及多个二级构造单元。这些不同级别构造单元在盆地结构上存在较大差异,根据这种差异性又可将该区划分为辽东湾区、 渤南区、 渤西区及渤中—渤东区4个区带(杨中轩,1984; 李涛等,2009; Huang et al., 2012)。渤海Z油田位于渤海湾渤西区,地处歧口凹陷,是渤海湾地区有利的石油富集区之一(图 1)。该区新生界地层从下至上包括孔店组、 沙河街组、 东营组、 馆陶组、 明化镇组及平原组,其中古近系沙河街组是Z油田主力含油气层系,该组地层从上至下划分为3段,分别为沙一段、 沙二段、 沙三段,沙一段又分为上、 下两部分,其中沙一下底部为一套湖相碳酸盐沉积,厚度约为35 m,富含油气,是主力产油气层之一,主要由生物石灰岩、 鲕粒石灰岩、 灰泥石灰岩和泥粉晶白云岩等组成(肖丽华等,2005; 朱筱敏等,2007; 杨有星等,2012; 李燕等,2014)。受构造演化影响,该地区断层交错,构造复杂(王光奇等,2003; 陈斌等,2005; 许立青等,2015),储层中裂缝发育。
渤海Z油田沙一下生物灰岩地层由于地貌原因无法进行露头观察,研究天然裂缝特征主要依靠岩心资料进行。根据31口井岩心裂缝描述及薄片观察,按照裂缝成因类型,该区生物灰岩储层发育构造裂缝和成岩裂缝两种类型,其中以构造裂缝为主。按照裂缝力学性质,可将沙一下储层构造裂缝进一步划分为剪切裂缝和扩张裂缝两类。剪切裂缝形成时,岩石破裂方向平行于裂缝面方向,裂缝面相对平直光滑,规模较大,产状稳定,缝面上可见擦痕、 阶步及羽饰构造等现象(图 2a、 图 2b、 图 2c); 不同于剪切裂缝,扩张裂缝形成时岩石破裂方向垂直于裂缝面方向,缝面粗糙不平,产状不稳定,延伸较短,发育程度有限,在岩心上比较少见,薄片上可见扩张裂缝具有明显的上述特征,并表现为明显的绕颗粒的分布特点(图 2d、 图 2e、 图 2f)。该区成岩裂缝主要为沉积物在成岩过程中由于压实和压溶等地质作用形成的规模较小的裂缝系统,如在薄片上可见层理缝一般为毫米至厘米级,多顺微层理面发育,呈锯齿状延伸,具有弯曲、 断续、 尖灭的特点(图 3a); 有些溶蚀裂缝呈不规则带状或串珠状展布(图 3b)。此外,还可见缝合线构造(图 3c),多与微层面平行展布,呈锯齿状展布,大多数被泥质或有机质充填。
上述构造裂缝和成岩裂缝由于成因机制和发育特征不同,对储层中流体的储集及渗流能力的影响存在较大差别,因此在注水开发过程中应区别加以对待。研究区构造裂缝规模较大,延伸范围广(平面上可达十几米至几十米),各向异性特征明显,对注水开发的影响最大,若注入水沿构造裂缝流动,极易引起方向性水窜甚至造成油井水淹,降低开发效果; 成岩裂缝规模较小,但在某些部位发育数量较多,对改善储层的储集性能作用明显,而注水过程中主要影响着局部渗流,对开发的影响相对构造裂缝来说较弱。根据这种差别,为了提高注水开发效果,须对构造裂缝发育特征进行深入研究,以为研究区制定有效的开发方案提供地质依据。
根据岩心及成像测井资料,Z油田沙一下储层主要发育北东—南西向、 北西—南东向及近东西向3组构造裂缝,以前两组为主(图 4),绝大多数构造裂缝为高角度缝,其中倾角大于60°的裂缝所占比例在80%以上,低角度构造缝基本不发育(图 5)。裂缝规模上(纵向高度、 平面长度及地下开度),由于岩心及成像测井均不能观测到裂缝的全貌且缺乏露头,若要获取相关参数需要采取其他方法: 1)裂缝纵向高度,实际工作中裂缝高度参数的获取可在划分岩石力学层之后直接在岩心上进行统计(曾联波, 2004,2008),根据这种方法得到该区构造裂缝高度一般小于170 cm,主要分布在10~30 cm以内,反映出裂缝主要在岩层内发育的特点; 2)裂缝平面长度,由于随机分布的裂缝平面长度具有分形特征(Pan et al., 2010; Parashar and Reeves, 2012),采取数学分形理论可求取裂缝长度,通过这种方法预测沙一下储层裂缝平面长度主要分布在35 m以内,某些规模较大裂缝的长度可达68 m; 3)裂缝地下开度,岩心或薄片观察可测得某些裂缝的地表开度,但这并不能代表地下环境的裂缝真实开度,需要对其进行校正,据高温高压三轴试验可知,裂缝开度与其所受的静封闭压力有关,随着裂缝面所受到的静封闭压力增大,裂缝开度呈负指数函数递减(曾联波等,1998),用这种关系可以间接获取裂缝地下开度; 通过统计沙一下储层岩心及薄片裂缝开度并进行校正,可知该区微裂缝开度在100 μm以内,主要分布在50 μm以内,部分规模较大的宏观裂缝开度可达几百微米甚至毫米级,这可能与灰岩裂缝的规模及溶蚀作用有关。
研究表明,Z油田沙一下储层裂缝的发育特征受岩性、 岩层厚度、 断层等地质因素控制作用明显。1)岩性的差异决定了岩石力学性质的差异,从而影响了裂缝的发育程度(曾联波等,1998)。根据岩心裂缝观察可知,该区白云岩类岩层裂缝发育程度最高,其次为白云质灰岩和泥灰岩类岩石,生物灰岩裂缝发育程度相对较差(图 6),这是由于随着岩石中方解石、 白云石等脆性矿物的含量越高,岩石脆性程度越大,裂缝的发育程度越高,而生物灰岩由于其物性最好,在相同构造应力场作用下容易产生塑性变形,不易形成构造裂缝。2)岩层厚度控制了天然裂缝规模,这是由于裂缝的发育受岩石力学层控制(Huang and Angelier, 1989),一般情况下岩石力学层与单一岩性层具有较好的对应关系,在一定的岩层厚度范围内,随着岩层厚度不断增大,裂缝规模越大,但裂缝线密度减小(图 7)。3)裂缝发育程度与断层的分布密切相关。通过对断层不同部位取心井裂缝统计表明,裂缝在断层上盘比下盘发育,在断层端部及转换部位比在断层中部更发育; 在距断层一定的距离范围内,裂缝发育程度较高,随着与断层面距离增大,裂缝发育程度变弱直至趋于稳定,这表明断层带附近不同部位由于受断层附近应力扰动的影响使得裂缝发育程度存在明显的差异。
由于研究区构造裂缝对注水开发的影响较大,因此本文主要从地质角度分析构造裂缝的成因机理,为研究裂缝的分布规律提供依据。根据岩心观察可知,该区沙一下储层构造裂缝并不是只在一期构造应力场作用下形成,主要有以下两方面的证据:
(1)从岩心构造裂缝倾角及缝面上的擦痕特征来看存在两种情况,即缝面两侧岩石的运动方向要么沿着裂缝倾向,要么沿着裂缝走向,说明裂缝的形成是在两种不同的构造应力环境下形成,反映裂缝至少存在两个形成时期。 图 2a显示裂缝面两侧的运动方向平行于裂缝面倾向,表明该组裂缝形成时的应力状态为最大主压应力垂直,中间主应力和最小主应力水平,在该应力状态下形成了与最大主压应力方向呈锐夹角的高角度剪切裂缝,为正断层型裂缝; 图 2b显示裂缝面两侧的运动方向沿裂缝走向,表明裂缝形成时最大主压应力为水平方向,与 图 2a中裂缝形成时的应力环境完全不同。
(2)从岩心及薄片观察来看,构造裂缝之间存在明显的切割、 限制及充填差异等关系,反映裂缝在不同时期形成的特点。如 图 8a所示A10井岩心发育3组构造裂缝相互切割,其中有两组裂缝明显近垂直相交特征,反映出至少在两期构造应力环境下形成; 图 8b所示薄片观察早期构造裂缝形成被矿物充填后,在后期的构造应力作用下再次形成晚期裂缝并切割早期裂缝,晚期裂缝无充填(被蓝色铸体充填),充分证明了裂缝是在两期不同的构造应力背景下形成。
那么上述构造裂缝分别是在哪几期构造应力场作用形成的呢?这需要根据裂缝特征并结合该区的构造演化背景进行分析。歧口凹陷位于黄骅坳陷中部,而黄骅坳陷是渤海湾新生代盆地的一级坳陷,该区在新生代始新世开始接受沉积,构造演化上经历了古近纪断陷和新近纪—第四纪拗陷两大阶段(吴元燕等,2000; 马寅生等,2002; 李燕等,2014)。其中始新世孔店期开始,在北西西—南东东向伸展(北北东—南南西向挤压)应力作用下,在盆地进入断陷初期,发育一系列北北东向同沉积断裂; 至渐新世早期(沙三段),构造应力逐渐右偏为北西—南东向伸展(北东—南西向挤压)作用转变,此时北东东向断裂活动加剧,裂陷进入扩张期,翘倾断块开始形成,发育大型的断鼻、 背斜及斜坡构造; 至沙二段—沙一段时期,进入断陷发育稳定期,此时同沉积断裂的方向以北东向为主; 至渐新世晚期(东营组),同沉积断裂活动逐渐减弱,断陷发育进入了衰退期; 至古近纪末期,即东营组沉积期及以后,断陷活动结束,此期间构造应力继续右偏,转变为北北西—南南东向伸展(北东东—南西西向挤压)作用,但伸展作用减弱,挤压作用增强,盆地发生反转,早期断裂活动由引张转变为走滑,并发育走滑断裂系统(马寅生等,2002; 周立宏等,2011)。新近纪时期开始,受水平方向拉伸和走滑应力控制的断陷及走滑作用被重力控制的热沉降作用所取代,盆地进入了裂后拗陷阶段。
根据上述分析可知,沙一段储层沉积以后,经历了始新世—渐新世裂陷中后期北西—南东向引张应力场及东营末期北东东—南西西向反转挤压构造应力场两期构造应力作用,综合沙一下储层构造裂缝特征分析判断,该区构造裂缝主要在以上两期构造应力场下形成。其中,成因类型为 图 2a的裂缝为早期裂缝,主要在始新世—渐新世裂陷中后期北西—南东向引张应力场作用下形成,与断层同生,其方位一般为北东—南西向,尤其在断层附近发育,裂缝走向与断层走向基本一致(图 9a); 成因类型为 图 2b的裂缝为晚期裂缝,主要在东营末期北东东—南西西向反转挤压构造应力场下形成,方位包括北东—南西向、 北西—南东向,为剪切裂缝,除这两组裂缝外,该期构造应力场还可以形成以近东西向为主的扩张裂缝(图 9b)。由于现今北东—南西向裂缝的形成经过了早、 晚两期裂缝的叠加,因此表现为该区优势裂缝方位(图 4)。
(1)渤海湾地区Z油田沙一下生物灰岩油藏主要发育构造裂缝及成岩裂缝两种类型,其中构造裂缝又分为剪切裂缝及扩张裂缝两类。大多数构造裂缝为高角度缝,主要方位为北东—南西向、 北西—南东向及近东西向; 裂缝纵向高度一般小于170 cm,裂缝平面长度主要分布在35 m以内,裂缝开度主要在100 μm以内; 裂缝发育特征受岩性、 岩层厚度、 断层等地质因素控制。
(2)根据裂缝的发育特征及构造演化特征综合分析可知,研究区构造裂缝主要在始新世—渐新世裂陷中后期北西—南东向引张应力场及东营末期北东东—南西西向反转挤压应力场两期构造应力作用下形成,其中早期裂缝主要方位为北东—南西向,晚期裂缝为北东—南西向、 北西—南东向及近东西向。
[1] | 陈斌,廖细明,杨宏伟等.2005. 周清庄油田精细储层预测.西南石油学院学报,27(3): 11-13. |
[2] | Chen Bin,Liao Ximing,Yang Hongwei et al.2005.The meticulous reservoir prediction of the Zhouqingzhuang oilfield.Journal of Southwest Petroleum Institute,27(3): 11-13. |
[3] | 陈飞,胡光义,范廷恩等.2015. 渤海海域W油田新近系明化镇组河流相砂体结构特征.地学前缘,22(2): 207-213. |
[4] | Chen Fei,Hu Guangyi,Fan Tingen et al.2015.Sandbody architecture and sequence stratigraphy of fluvial facies,Neogene Minghuazhen Formation,W oilfield,Bohai Bay.Earth Science Frontiers,22(2): 207-213. |
[5] | 陈淑利,孙庆和,宋正江.2008. 特低渗透裂缝型储层注水开发中后期地应力场变化及开发对策.现代地质,22(4): 64-72. |
[6] | Chen Shuli,Sun Qinghe and Song Zhengjiang.2008.Changes of ground stress field and development policy in later half period of water injection for the extremely low permeability reservoirs with fractures.Geoscience,22(4): 64-72. |
[7] | 李理,赵利,刘海剑等.2014. 渤海湾盆地晚中生代-新生代伸展和走滑构造及深部背景.地质科学,50(2): 446-472. |
[8] | Li Li,Zhao Li,Liu Haijian et al.2014.Late Mesozoic to Cenozoic extension and strike slip structures and deep background of Bohai Bay Basin.Chinese Journal of Geology,50(2): 446-472. |
[9] | 李涛,王宗秀.2009. 华北上地壳构造单元划分.地学前缘,16(4): 216-224. |
[10] | Li Tao and Wang Zongxiu.2009.Tectonic division of the upper crust of the North China block.Earth Science Frontiers,16(4): 216-224. |
[11] | 李燕,金振奎,杨有星等.2014. 歧口凹陷沙一下亚段湖相碳酸盐岩储层质量差异机理.天然气地球科学,25(2): 201-211. |
[12] | Li Yan,Jin Zhenkui,Yang Youxing et al.2014.Mechanism of reservoir quality difference of lacustrine carbonate in Lower Member 1 Shahejie Formation in Qikou sag.Natural Gas Geoscience,25(2): 201-211. |
[13] | 李新景,胡素云,程克明.2007. 北美裂缝性页岩气勘探开发的启示.石油勘探与开发,34(4): 647-654. |
[14] | Li Xinjing,Hu Suyun and Cheng Keming.2007.Suggestions from the development of fractured shale gas in North America.Petroleum Exploration and Development,34(4): 392-400. |
[15] | 李跃纲,巩磊,曾联波等.2012. 四川盆地九龙山构造致密砾岩储层裂缝特征及其贡献.天然气工业,32(1): 22-26. |
[16] | Li Yuegang,Gong Lei,Zeng Lianbo et al.2012.Characteristics of fractures and their contribution to the deliverability of tight conglomerate reservoirs in the Jiulongshan structure,Sichuan Basin.Nature Gas Industry,32(1): 22-26. |
[17] | 梁官忠,王群超,夏欣等.2001. 阿北安山岩油藏裂缝发育特征对开发的影响.石油勘探与开发,28(2): 64-72. |
[18] | Liang Guanzhong,Wang Qunchao,Xia Xin et al.2001.The effect of growth fracture of Abei andesite reservoir on oilfield development.Petroleum Exploration and Development,28(2): 64-72. |
[19] | 马寅生,曾庆利,张兴.2002. 黄骅坳陷新生代构造应力场演化的光弹模拟与石油地质条件分析.地质力学学报,8(3): 219-228. |
[20] | Ma Yingsheng,Zeng Qingli and Zhang Xing.2002.Photoelastic experimental modelling of evolution of Cenozoic tectonic stress field and analysis of petroleum geological condition in Huanghua depression.Journal of Geomechanics,8(3): 219-228. |
[21] | 滕长宇,邹华耀,郝芳.2014. 渤海湾盆地构造差异演化与油气差异富集.中国科学:地球科学,44(4): 579-590. |
[22] | Teng Changyu,Zou Huayao and Hao Fang.2014. Control of differential tectonic evolution on petroleum occurrence in Bohai Bay Basin.Science China Earth Sciences,57(5): 1117-1128. |
[23] | 王珂,张荣虎,戴俊生等.2015. 低渗透储层裂缝研究进展.地球科学与环境学报,37(2): 44-58. |
[24] | Wang Ke,Zhang Ronghu,Dai Junsheng et al.2015.Review on low-permeability reservoir facture.Journal of Earth Sciences and Environment,37(2): 44-58. |
[25] | 王光奇,漆家福,岳云福.2003. 歧口凹陷及周缘新生代构造的成因和演化.地质科学,38(2): 230-240. |
[26] | Wang Guangqi,Qi Jiafu and Yue Yunfu.2003.Formation and evolution of the Cenozoic tectonics within and surrounding the Qikou sag.Chinese Journal of Geology,38(2): 230-240. |
[27] | 吴元燕,付建林,周建生等.2000. 歧口凹陷含油气系统及其评价.石油学报,21(6): 18-22. |
[28] | Wu Yuanyan,Fu Jianlin,Zhou Jiansheng et al.2000.Evaluation of hydrocarbon system in Qikou sag.Acta Petrolei Sinica,21(6): 18-22. |
[29] | 肖丽华,孟元林,牛嘉玉等.2005. 歧口凹陷沙河街组成岩史分析和成岩阶段预测.地质科学,40(3): 346-362. |
[30] | Xiao Lihua,Meng Yuanlin,Niu Jiayu et al.2005.Diagenetic history and diagenetic stages prediction of Shahejie Formation in the Qikou sag.Chinese Journal of Geology,40(3): 346-362. |
[31] | 许立青,李三忠,索艳慧等.2015. 渤海湾盆地大岐口凹陷断裂系统与陆内拉分断陷.地质科学,50(2): 489-502. |
[32] | Xu Liqing,Li Sanzhong,Suo Yanhui et al.2015.Fault system and basin prototype of the Great Qikou sag,Bohai Bay Basin.Chinese Journal of Geology,50(2): 489-502. |
[33] | 杨明慧.2009. 渤海湾盆地变换构造特征及其成藏意义.石油学报,30(6): 817-823. |
[34] | Yang Minghui.2009.Transfer structure and its relation to hydrocarbon exploration in Bohai Bay Basin.Acta Petrolei Sinica,30(6): 817-823. |
[35] | 杨永林,张莉,杨亚娟等.2004. 鄯善油田裂缝体系及对注水开发效果的影响.西北大学学报(自然科学版),34(1): 97-100. |
[36] | Yang Yonglin,Zhang Li,Yang Yajuan et al.2004.The fracture system and its effect on waterflooding in Shanshan oilfield.Journal of Northwest University(Natural Science Edition),34(1): 97-100. |
[37] | 杨有星,金振奎,高白水.2012. 歧口凹陷缓坡区湖泊碳酸盐岩优质储层特征及成因.石油学报,34(1): 97-100. |
[38] | Yang Youxing,Jin Zhenkui and Gao Baishui.2012.Characteristics and origin of lacustrine carbonate high-quality reservoirs in the slope area of Qikou sag.Acta Petrolei Sinica,34(1): 97-100. |
[39] | 杨中轩.1984. 渤海湾盆地的构造单元及构造成因.石油勘探与开发,(2): 27-31. |
[40] | Yang Zhongxuan.1984.The tectonic elements and the origin of the structures in Bohai Bay Basin.Petroleum Exploration and Development,(2): 27-31. |
[41] | 袁士义,宋新民,冉启泉.2004. 裂缝性油藏开发技术.北京:石油工业出版社.1-335. |
[42] | Yuan Shiyi,Song Xinmin and Ran Qiquan.2004.Development Technology of Fractured Reservoir.Beijing: Petroleum Industry Press.1-335. |
[43] | 曾联波,任德生.1998. 大民屯凹陷静52断块油藏裂缝参数研究.断块油气田,5(4): 29-32. |
[44] | Zeng Lianbo and Ren Desheng.1998.Study of facture parameter of Jing 52 faulted reservoir in Damingtun depression.Fault-Block Oil and Gas Field,5(4): 29-32. |
[45] | 曾联波.2004. 低渗透砂岩油气储层裂缝及其渗流特征.地质科学,39(1): 11-17. |
[46] | Zeng Lianbo.2004.Fissure and its seepage characteristics in low-permeable sandstone reservoir.Chinese Journal of Geology,39(1): 11-17. |
[47] | 曾联波,李忠兴,史成恩等.2007. 鄂尔多斯盆地上三叠统延长组特低渗透砂岩储层裂缝特征及成因.地质学报,81(2): 175-180. |
[48] | Zeng Lianbo,Li Zhongxing,Shi Cheng'en et al.2007.Characteristics and origin of fractures in the extra low-permeability sandstone reservoirs of the Upper Triassic Yanchang Formation in the Ordos Basin.Acta Geologica Sinica,81(2): 175-180. |
[49] | 曾联波.2008. 低渗透砂岩储层裂缝的形成与分布.北京:科学出版社.1-169. |
[50] | Zeng Lianbo.2008.Formation and Distribution of Factures in Low-Permeability Sandstone Reservoir.Beijing: Science Press.1-169. |
[51] | 曾联波,柯式镇,刘洋等.2010. 低渗透油气储层裂缝研究方法.北京:石油工业出版社.1-187. |
[52] | Zeng Lianbo,Ke Shizheng,Liu Yang et al.2010.Study Methods of Facture in Low-Permeability Oil and Gas Reservoir.Beijing: Petroleum Industry Press.1-187. |
[53] | 赵向原,曾联波,王晓东等.2015. 鄂尔多斯盆地宁县-合水地区长6、长7、长8储层裂缝差异性及开发意义.地质科学,50(1): 274-285. |
[54] | Zhao Xiangyuan,Zeng Lianbo,Wang Xiaodong et al.2015.Differences of natural fracture characteristics and their development significance in Chang 6,Chang 7 and Chang 8 reservoir,Ningxian-Heshui area,Ordos Basin.Chinese Journal of Geology,50(1): 274-285. |
[55] | 周立宏,卢异,肖敦清.2011. 渤海湾盆地歧口凹陷盆地结构构造及演化.天然气地球科学,22(3): 373-382. |
[56] | Zhou Lihong,Lu Yi and Xiao Dunqing.2011.Basinal texture structure of Qikou sag in Bohai Bay Basin and its evolution.Natural Gas Geoscience,22(3): 373-382. |
[57] | 朱筱敏,董艳蕾,郭长敏等.2007. 歧口凹陷沙河街组一段层序格架和储层质量分析.沉积学报,25(6): 935-941. |
[58] | Zhu Xiaomin,Dong Yanlei,Guo Changmin et al.2007.Sequence framework and reservoir quality of Sha 1 Member in Shahejie Formation,Qikou sag.Act Sedmentologica Sinica,25(6): 935-941. |
[59] | Huang L,Liu C Y,Zhou X H et al.2012. The important turning points during evolution of Cenozoic Basin offshore the Bohai Sea: Evidence and regional dynamics analysis.Science China Earth Sciences,55(3): 476-487. |
[60] | Huang Q and Angelier J.1989. Fracture spacing and its relation to bed thickness.Geological Magazine,126(4): 355-362. |
[61] | Pan J B,Lee C C,Lee C H et al.2010. Application of fracture network model with crack permeability tensor on flow and transport in fractured rock.Engineering Geology,116(1-2): 166-177. |
[62] | Parashar R and Reeves D M.2012. On iterative techniques for computing flow in large two-dimensional discrete fracture networks.Journal of Computational and Applied Mathematics,236(18): 4712-4724. |
2. Exploration & Production Research Institute, SINOPEC, Beijing 100083;
3. Research Institute of Shenzhen Branch Company, CNOOC, Guangzhou 510240