2. 山东省油藏地质重点实验室, 山东青岛 266580;
3. 中国石油大学(华东)深层油气重点实验室, 山东青岛 266580;
4. 胜利油田东胜精攻石油开发集团股份有限公司, 山东东营 257000
2. Reservoir Geology Key Laboratory of Shandong Province, Qingdao 266580, China;
3. Key Laboratory of Deep Oil and Gas in China University of Petroleum (East China), Qingdao 266580, China;
4. Dongsheng Petroleum Development Group Company Limited of Shengli Oilfield, Dongying 257000, China
天然裂缝的成因受多种因素的综合影响, 其发育状态具有很强的不确定性, 分布规律具有很强的随机性[1-4]。如何更加准确地刻画储层中裂缝的空间展布, 仅从构造或者岩相单一手段分析裂缝的发育程度存在明显缺陷[5-6]。从裂缝发育的岩石来看, 不同部位岩石脆性组分的含量、岩石结构、岩石粒度等的差异导致对裂缝控制作用的不同, 即沉积相控制下的岩性为主要影响因素[7-8]; 从裂缝形成机制来看, 采用有限元数值模拟构造应力场分析方法, 认为应力分布情况是主要影响因素[9-14]; 从裂缝与周围构造体间关系来看[15-16], 裂缝发育位置距断层远近为主要影响因素; 从地震资料入手, 利用蚂蚁追踪、相干体属性等分析方法对识别裂缝具有很好的效果[17]。从建模手段来看, 逐渐从连续裂缝模型发展到离散裂缝网络模型[18-21]。针对致密砂岩储层, 笔者采用基于示性点过程的离散裂缝网络建模方法, 在裂缝发育地质模式的控制下, 综合裂缝成因力学机制, 利用裂缝密度体建立天然裂缝模型。
1 裂缝建模研究思路传统裂缝建模方法对成像测井、岩心的依赖程度较大, 存在成本高, 井间误差大等问题。裂缝密度体的提出可以有效解决缺乏特殊测井和取心井, 以及井少、井距大的情况下裂缝预测与建模的问题。同时分尺度充分发挥不同信息的准确度, 尤其对于中小尺度裂缝, 在建立断层相关裂缝发育概念模式的约束下, 依据岩石力学试验探究裂缝成因, 从机制角度上对断层相关裂缝发育模式进行解释; 同时利用构造应力场有限元数值模拟手段, 对研究区造缝期应力场分布进行模拟, 基于岩石破裂准则定量化表征裂缝空间分布。综合岩相、构造和应力场模拟结果, 建立三维裂缝密度体, 采用离散裂缝网络建模方法, 建立天然裂缝空间模型(图 1)。
大芦湖油田樊162地区处于东营凹陷博兴洼陷东部(图 2)。博兴洼陷西靠高青凸起, 北邻平南构造, 南接鲁西隆起, 东部以石村断裂带分界。洼陷内部继承性地发育金家-樊家鼻状构造, 将洼陷划分为东西两个次级洼陷, 樊162地区位于鼻状构造带的东翼, 高青断层下降盘。沙三段区域构造面貌较为简单, 总体呈南高北低, 倾角3°~4°的单斜。在控洼断层高青-平南断层、博兴断层的共同作用下, 区域地质应力以近南北向伸展为主兼有右旋张扭应力。区内发育一系列北东东向及近东西向的盆倾正断层, 落差为10~400 m。
樊162储层中天然裂缝发育主要方向为北东东向及近南北向, 少量裂缝发育北东向。该地区受高青-平南断层、博兴断层的控制, 区域地质应力以近南北向伸展为主。
根据取心井岩心观察统计结果, 取心井均发育天然裂缝, 发育裂缝程度100%。裂缝类型以构造缝为主, 其中又以同时受剪应力和张应力作用下的张剪性裂缝为主, 发育位置多集中于粉细砂岩的朵叶体微相中。主要发育高角度裂缝, 裂缝倾角多集中在65°~90°, 以近乎垂直的倾角为主。发育少量的斜交裂缝, 未见低角度裂缝发育。裂缝开度范围在40 μm~2 mm, 主要集中在40~80 μm, 整体上开度较小, 闭合缝占比较高。开度较大的裂缝绝大多数被白云石充填, 占比15%。裂缝面平直, 有擦痕、阶步等特征(图 3)。通过岩心对裂缝密度进行分析得出, 研究区平均裂缝密度3条/m, 裂缝密度超过3条/m的取心井数占到了75%, 属天然裂缝较发育地带。
岩性差异性是影响裂缝发育的内因。沉积作用导致岩石组分、岩性组合、结构、单层厚度等的差异, 进而导致裂缝发育程度的差异性。建立研究区岩相空间展布, 为裂缝预测提供约束条件。采用沉积微相约束, 序贯指示模拟的方法, 建立起不同沉积微相内裂缝发育程度的差异性与相关性。
研究区为洪水型湖底扇沉积, 主要为水道砂体, 物源来自东南方向。发育主水道、水下天然堤、辫状水道、辫状水道间、水道前缘砂坝和朵叶体等6种微相类型。在沉积微相平面展布上, 不同时期水道主要呈北西-南东向, 树枝状分布, 物源逐渐向东偏移; 辫状河道间砂体分布广泛; 水道前缘砂坝分布在辫状水道末梢端(图 4)。由于物源来自东南方向, 使得砂岩含量及砂岩厚度呈现出由东南向西北逐渐减小、减薄的现象,导致水道末梢端的水道前缘砂坝砂岩厚度较小。
分别统计裂缝在不同沉积微相下岩性的分布以及发育裂缝密度之间的关系(图 5)。在水道前缘砂坝、朵叶体中以含有脆性矿物较多的粉砂岩、细砂岩为主, 因朵叶体中发育砂泥岩互层,裂缝在朵叶体中发育一般, 在水道前缘砂坝较为发育; 辫状水道间裂缝发育程度次之, 岩性主要以中细砂岩为主, 粒度变粗; 水道中裂缝发育程度最差, 岩性以砂砾岩、含砾中粗砂岩为主。
在沉积微相的约束下, 单井上岩性解释作为硬数据, 基于序贯指示模拟算法, 建立岩相模型, 其概率分布一致性较高(图 6)。
分析距断层远近与裂缝密度间关系, 随着距主断层距离的逐渐增大, 裂缝线密度呈现逐渐递减趋势(图 7); 同时裂缝开度呈现逐渐递减的趋势。从图 7中可以看出, 随着距主断层距离的逐渐增大, 裂缝密度比裂缝开度更具敏感, 其相关性高于裂缝开度。这表明在断层形成后, 断裂表面的有效应力与差应力同时减小, 使得裂缝开度变化上所需能量明显低于裂缝密度变化时的能量。
从岩性、层厚、所处构造位置等因素的综合分析, 基于天然裂缝发育特征建立三维空间裂缝发育概念模式(图 8, 其中Shmin为水平最小主应力, Shmax为水平最大主应力, Sv为垂向主应力)。将断层内部结构划分为断层核、诱导裂缝带; 按照受力状态划分为剪切带、破碎带、弱变形带。
选取研究区取心井岩心样品共计44块, 进行单轴压缩、巴西劈裂以及常规三轴实验, 其中单轴压缩与常规三轴实验样品为直径25 mm, 高度50 mm的圆柱; 巴西劈裂实验采用直径25 mm, 高度16.6 mm的圆柱, 表面粗糙度大于5级(Ra < 2.5 μm, GB1031-83), 垂向偏差角度小于0.11°, 水平偏差小于5 μm。
单轴压缩实验共计7组, 每组两个样品取平均(表 1); 巴西劈裂实验共计6组, 每组3个样品取平均(表 2); 常规三轴压缩实验共计4组, 每组设置在5、10、15 MPa下各一个样品(表 3)。
取自靠近断层核部密度最大的断裂角砾岩, 其轴向抗压强度最大, 弹性模量最大; 在平面上, 随着距主断层距离的增加, 岩性粒度逐渐变细, 密度逐渐变小, 应力应变曲线呈现出轴向抗压强度逐渐降低, 应变值逐渐变大的趋势; 在经过断层部位的垂向上, 岩性粒度变化由上到下从细粒到粗粒再到细粒, 密度从小到大再到小, 轴向抗压强度从小到大再到小, 弹性模量从小到大再到小, 应变值从大到小再到大的趋势(图 9(a))。
宏观裂缝形成之前是大量微裂缝聚集的过程。结合单轴压缩下岩样的应力应变曲线, 将微裂缝演化至宏观裂缝划分为5个阶段:第1阶段为初始压实阶段, 对应应力应变曲线的OA段, 该阶段岩石内部裂隙逐渐被压缩闭合, 但卸载后可完全恢复; 第2阶段为线弹性变形阶段, 曲线接近直线, 对应AB段; 第3阶段曲线发生偏离, 出现塑性变形, 岩石内部出现平行最大主应力方向的微裂隙, 应力越大, 裂隙越多, 对应曲线的BC段; 第4阶段岩石内部裂隙迅速增多, 裂隙密度增大, 即将到达岩石最大承载力, 对应曲线CD段; 第5阶段岩石内部微裂隙被贯通, 裂隙密度降低, 裂缝开度增大, 应力释放, 形成宏观裂缝, 对应曲线DE段(图 9(b))。
2.4.3 构造应力场数值模拟研究区在三叠纪由于华北、华南板块碰撞导致NWW-NW向压性构造发育; 侏罗纪晚期至白垩纪初期, 太平洋板块向欧亚板块以NNW斜向低角度(10°)俯冲, 引起郯庐断裂左行走滑, 产生右旋张扭构造应力场控制下的NE和NW向两组断陷的发育; 古近纪时期, 太平洋板块向欧亚板块俯冲, 引起大规模的伸展拉张, 产生NE、NEE、NW、NWW向断陷的发育。综合构造运动演化, 最终确定在喜山期研究区整体受到拉张应力, 兼具右旋剪切力。通过施加不同的边界载荷, 反复调试, 最终确定在模型的左上和右下边界施加68 MPa最大主张应力, 左下和右上边界施加32 MPa最小主压应力。为了模拟研究区喜山期所受到的右旋剪切应力, 对模型的左上角与右下角施加一组垂直于最大主应力方向的剪切力6 MPa。同时对模型底部Z方向, 模型左角与右角施加X、Y方向的位移约束(图 10)。
模拟结果显示:研究区水平最小主应力为44~87 MPa(图 11(a)); 水平中间主应力为22~58 MPa(图 11(b)); 水平最大主应力为-5~40 MPa(图 11(c))。在断层核部位置及周围出现明显应力异常; 断层核部因岩石已破坏, 应力得到释放, 出现明显应力低值; 断层核部周围环带出现明显应力高值, 并且呈现出环带状递减趋势, 这与岩石破裂过程应力-应变曲线一致, 进一步验证了裂缝发育分带性特征。总体规律呈现断层核部处于应力低值区, 诱导裂缝带上盘部位更易出现应力高值。3个水平应力的分布基本相似, 呈现出受断层控制的状态。
在拉张应力和右旋剪切力的作用环境下, 采用格里菲斯准则和库仑-摩尔准则对是否破裂进行定性判别, 引入格里菲斯二维破裂准则的计算公式。
当σ1+3σ3≥0时, 破裂准则为
$ {\left( {{\sigma _1} - {\sigma _3}} \right)^2} - 8\left( {{\sigma _1} + {\sigma _3}} \right){\sigma _{\rm{T}}} = 0. $ | (1) |
当σ1+3σ3 < 0时, 破裂准则为
$ {\sigma _3} + {\sigma _{\rm{T}}} = 0. $ | (2) |
式中, σ1、σ3分别为最大和最小主应力, MPa; σT为张应力, MPa。
库仑-摩尔定律是指当材料沿破坏产生的剪切面破裂, 其剪应力τ既要克服材料本身的内聚力C, 也要克服作用在剪切面上的摩擦力σtan φ。引入库仑-摩尔破裂准则的函数:
$ \left| \tau \right| = C + \sigma {\rm{tan}}\varphi . $ | (3) |
式中, |τ|为材料抗剪强度, MPa; C为内聚力, 材料本身固有的剪切强度, MPa; σ为作用在物体上的法向力, 即正应力, MPa; φ为内摩擦角, (°); tan φ为内摩擦系数。
结合研究区主要发育高角度张剪性裂缝的特征, 基于构造应力场模拟结果, 采用格里菲斯和摩尔-库仑准则相结合的方式半定量判断裂缝发育程度。引入张破裂系数和剪破裂系数, 其计算公式为
$ {I_{\rm{t}}} = {\sigma _{\rm{T}}}/{\sigma _{\rm{t}}}, $ | (4) |
$ {I_{\rm{n}}} = \tau /\left| \tau \right|. $ | (5) |
式中, It为张破裂系数; σt为材料抗张强度, MPa; In为剪破裂系数; τ为剪应力, MPa。
实际情况下裂缝破裂绝大多数并未达到抗张强度或抗剪强度即开启, 且考虑实际研究区裂缝特征为张剪性, 介于剪破裂与张破裂之间, 因此为了更加准确地表征致密砂岩中构造裂缝的发育程度, 综合张破裂系数和剪破裂系数, 构建综合破裂指数R:
$ R = \left( {a{I_{\rm{t}}} + b{I_{\rm{n}}}} \right)/2. $ | (6) |
依据综合破裂指数, 以单井上裂缝密度作为硬数据, 以综合破裂指数作为空间约束, 采用序贯指示方法, 分层随机模拟出裂缝密度, 圈定划分出储层中裂缝发育等级区域, 分发育区、发育过渡区和不发育区3个等级(表 4)。
在研究区西北部和中部广泛分布有Ⅰ~Ⅱ级裂缝发育区, 综合破裂指数大于1.5, 裂缝密度约为4, 其中在构造隆起的高部位, 因地层弯曲程度比其他部位要大, 该地区综合破裂指数最大, 属裂缝发育Ⅰ区。在研究区西南部主要分布Ⅱ~Ⅲ级裂缝发育区, 综合破裂指数分布在1.0~1.5, 裂缝密度约1.5~4条/ m。由于研究区南部很少受构造活动的影响, 地势相对平坦, 砂体厚度较北部厚, 该地区综合裂缝发育指数最小, 裂缝密度约为1.5条/m, 属裂缝发育过渡Ⅱ区和不发育Ⅲ区(图 12)。
基于前期裂缝发育概念模式的建立, 引入距断层远近裂缝发育指数对概念模式进行定量化表征, 其计算公式为
$ z = \sigma \left[ {1 - {e^{\lambda (1 - x)}}} \right]. $ | (7) |
式中, z为裂缝发育程度, z值处在0~1之间, 1代表裂缝发育程度高, 0代表裂缝不发育; σ为比例因子; λ为裂缝发育程度变化率; x为距断层距离, m。
依据距断层远近裂缝发育指数, 划分裂缝发育Ⅰ区、发育过渡Ⅱ区、不发育Ⅲ区3个裂缝发育区域, 划分标准如表 5所示。
以断层中轴线作为计算起始点(图 13(a)), 建立距断层远近的三维属性体(图 13(b))。依据距断层远近的裂缝发育指数, 在属性体的基础上划分为3个区域, 其中红色区域为发育Ⅰ区, 裂缝发育; 绿色区域为发育过渡Ⅱ区, 裂缝一般发育; 紫色区域为不发育Ⅲ区, 裂缝不发育(图 13(c))。采用离散裂缝网络建模方法构建裂缝三维模型(图 13(d)), 较好展示裂缝与断层间关系。
在岩相模型的控制下, 将综合破裂指数控制下的裂缝密度模型, 与距断层远近指数控制下的裂缝密度模型, 采用趋势建模手段, 进行融合(图 14)。
离散裂缝网络建模是目前建立天然裂缝模型较为常用的方法, 通过依次定义裂缝空间分布、几何形态、开度和方位来建立接近真实地下裂缝空间分布模型。对于裂缝空间分布, 采用趋势建模建立的裂缝密度体作为输入; 裂缝几何形态定义为面片状四边形, 伸展率为2, 裂缝长度服从指数分布; 开度采用岩心实测校正后数据; 方位采用应力场模拟结果进行约束。建立的离散裂缝网络模型如图 15所示。
选取裂缝密度分布跨度较大的5口井, 以每口井轨迹中线为中心点, 周围25 m内4个网格内的裂缝密度作为统计数据源, 得到5口井裂缝密度依次为0.7、1.2、2.9、3.2、4.0。与动态生产数据中这5口井的含水率做相关性分析结果也可看出, 井点附近裂缝密度与含水率具有较好的线性相关性, 相关系数为0.817 2(图 16)。实际生产动态资料与建立的裂缝模型具有很好的相关性, 验证了所建裂缝模型的有效性。
(1) 以粉、细砂岩为主的水道前缘砂坝微相, 裂缝较为发育, 其微相岩性具有粒度细、孔隙度低、脆性矿物含量高的特征。中扇亚相以中细砂岩为主的辫状河道间微相, 裂缝发育一般。内扇中以砂砾岩为主的河道微相, 裂缝不易发育。
(2) 中国东部正断层相关裂缝发育具有如下特征:裂缝密度与距离断层远近呈反比趋势; 断层上盘相较下盘更易裂缝发育; 砂岩中裂缝发育程度大于泥岩中裂缝发育程度; 薄层砂体裂缝发育程度大于厚层砂体。应力聚集阶段对应原状地层, 应力集中阶段对应诱导裂缝带, 应力释放阶段对应断层核。
(3) 在研究区西北部和中部, 构造隆起高部, 断层上盘位区域, 裂缝较为发育。在研究区东南部, 地势较为平坦, 砂体较厚部位, 裂缝发育较少。
(4) 在裂缝发育部位具有高脆性、近断层、占高点、主上盘、薄砂层和沿长轴的特征。
[1] |
NELSON R A. Geologic analysis of naturally fractured reservoirs[M]. Houston: Gulf Professional Company, 2001: 12-15.
|
[2] |
BOURBIAUX B, BASQUET R, DANIEL J M, et al. Fractured reservoirs modelling: a review of the challenges and some recent solutions[J]. First Break, 2005, 23(1087): 1395-1410. |
[3] |
马存飞, 董春梅, 栾国强, 等. 泥页岩自然流体压力缝类型、特征及其作用:以中国东部古近系为例[J]. 石油勘探与开发, 2016, 43(4): 580-589. MA Cunfei, DONG Chunmei, LUAN Guoqiang, et al. Types, characteristics and effects of natural fluid pressure fractures in shale: a case study of the Paleogene strata in Eastern China[J]. Petroleum Exploration and Development, 2016, 43(4): 580-589. |
[4] |
MA C, ELSWORTH D, DONG C, et al. Controls of hydrocarbon generation on the development of expulsion fractures in organic-rich shale: based on the Paleogene Shahejie Formation in the Jiyang Depression, Bohai Bay Basin, East China[J]. Marine and Petroleum Geology, 2017, 86(7): 1406-1416. |
[5] |
李阳, 侯加根, 李永强. 碳酸盐岩缝洞型储集体特征及分类分级地质建模[J]. 石油勘探与开发, 2016, 43(4): 600-606. LI Yang, HOU Jiagen, LI Yongqiang. Features and hierarchical modeling of carbonate fracture-cavity reservoirs[J]. Petroleum Exploration and Development, 2016, 43(4): 600-606. |
[6] |
樊建明, 屈雪峰, 王冲, 等. 鄂尔多斯盆地致密储集层天然裂缝分布特征及有效裂缝预测新方法[J]. 石油勘探与开发, 2016, 43(5): 740-748. FAN Jianming, QU Xuefeng, WANG Chong, et al. Natural fracture distribution and a new method predicting effective fractures in tight oil reservoirs of Ordos Basin, NW China[J]. Petroleum Exploration and Development, 2016, 43(5): 740-748. |
[7] |
林承焰, 梁昌国, 任丽华, 等. 碎屑岩潜山沉积相及其与裂缝发育的关系[J]. 中国石油大学学报(自然科学版), 2008, 32(4): 1-6. LIN Chengyan, LIANG Changguo, REN Lihua, et al. Sedimentary facies and relation to degree of fracture development of clastic buried hills[J]. Journal of China University of Petroleum (Edition of Natural Science), 2008, 32(4): 1-6. DOI:10.3321/j.issn:1673-5005.2008.04.001 |
[8] |
周永胜, 张流. 裂缝性储集层岩芯裂缝统计分析[J]. 世界地质, 2000, 19(2): 117-124. ZHOU Yongsheng, ZHANG Liu. Statistics and analysis of drilling core fractures of fractured reservoirs[J]. World Geology, 2000, 19(2): 117-124. DOI:10.3969/j.issn.1004-5589.2000.02.003 |
[9] |
KAI J, DING W L, HUANG W H, et al. Fractures of lacustrine shale reservoirs, the Zhanhua Depression in the Bohai Bay Basin, eastern China[J]. Marine and Petroleum Geology, 2013, 48(5): 113-123. |
[10] |
DING W, FAN T, YU B, et al. Ordovician carbonate reservoir fracture characteristics and fracture distribution forecasting in the Tazhong area of Tarim Basin, Northwest China[J]. Journal of Petroleum Science & Engineering, 2012, 87(3): 62-70. |
[11] |
王珂, 张惠良, 张荣虎, 等. 超深层致密砂岩储层构造裂缝特征及影响因素:以塔里木盆地克深2气田为例[J]. 石油学报, 2016, 37(6): 715-727. WANG Ke, ZHANG Huiliang, ZHANG Ronghu, et al. Characteristics and influencing factors of ultra-deep tight sandstone reservoir structural fracture: a case study of Keshen-2 gas field, Tarim Basin[J]. Acta Petrolei Sinica, 2016, 37(6): 715-727. |
[12] |
徐珂, 戴俊生, 商琳, 等. 高尚堡油田深层油藏南区现今地应力场预测及应用[J]. 中国石油大学学报(自然科学版), 2018, 42(6): 19-29. XU Ke, DAI Junsheng, SHANG Lin, et al. Prediction of current in-situ stress field and its application of southern area of deep reservoir in Gaoshangpu Oilfield[J]. Journal of China University of Petroleum (Edition of Natural Science), 2018, 42(6): 19-29. DOI:10.3969/j.issn.1673-5005.2018.06.003 |
[13] |
戴俊生, 汪必峰, 马占荣. 脆性低渗透砂岩破裂准则研究[J]. 新疆石油地质, 2007, 28(4): 393-395. DAI Junsheng, WANG Bifeng, MA Zhanrong. Research on cracking principles of brittle-permeability sands[J]. Xinjiang Petroleum Geology, 2007, 28(4): 393-395. DOI:10.3969/j.issn.1001-3873.2007.04.001 |
[14] |
杨少春, 牛海瑞, 宋明水, 等. 车排子地区石炭系挤压逆冲构造区断层共生裂缝发育程度定量表征[J]. 中国石油大学学报(自然科学版), 2017, 41(5): 1-8. YANG Shaochun, NIU Hairui, SONG Mingshui, et al. Quantitative characterization of development degree of fault-related fracture in compressional thrust structure block of the Carboniferous in Chepaizi area[J]. Journal of China University of Petroleum (Edition of Natural Science), 2017, 41(5): 1-8. DOI:10.3969/j.issn.1673-5005.2017.05.001 |
[15] |
GHOSH K, MITRA S. Structural controls of fracture orientations, intensity, and connectivity, Teton anticline, Sawtooth Range, Montana[J]. AAPG Bulletin, 2009, 93(8): 995-1014. DOI:10.1306/04020908115 |
[16] |
FANG J, ZHOU F, TANG Z. Discrete fracture network modelling in a naturally fractured carbonate reservoir in the Jingbei Oilfield, China[J]. Energies, 2017, 10(2): 183-198. DOI:10.3390/en10020183 |
[17] |
吴键, 孙圆辉, 王彬, 等. 准噶尔DX18区块裂缝性火山岩储集体三维地质建模[J]. 石油勘探与开发, 2012, 39(1): 92-99. WU Jian, SUN Yuanhui, WANG Bin, et al. 3D geological modeling of fractured volcanic reservoir bodies in Block DX18 in Jungger Basin, NW China[J]. Petroleum Exploration and Development, 2012, 39(1): 92-99. |
[18] |
姚军, 王子胜, 张允, 等. 天然裂缝性油藏的离散裂缝网络数值模拟方法[J]. 石油学报, 2010, 31(2): 284-288. YAO Jun, WANG Zisheng, ZHANG Yun, et al. Numerical simulation method of discrete fracture network for naturally fractured reservoirs[J]. Acta Petrolei Sinica, 2010, 31(2): 284-288. |
[19] |
彭仕宓, 索重辉, 王晓杰, 等. 整合多尺度信息的裂缝性储层建模方法探讨[J]. 西安石油大学学报(自然科学版), 2011, 26(4): 1-7. PENG Shimi, SUO Chonghui, WANG Xiaojie, et al. A modeling method for fractured reservoirs using multi-scale information[J]. Journal of Xi' an Shiyou University (Natural Science Edition), 2011, 26(4): 1-7. DOI:10.3969/j.issn.1673-064X.2011.04.001 |
[20] |
郎晓玲, 郭召杰. 基于DFN离散裂缝网络模型的裂缝性储层建模方法[J]. 北京大学学报(自然科学版), 2013, 49(6): 964-972. LANG Xiaoling, GUO Zhaojie. Fractured reservoir modeling method based on discrete fracture network model[J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 2013, 49(6): 964-972. |
[21] |
赵力彬, 石石, 肖香姣, 等. 库车坳陷克深2气藏裂缝-孔隙型砂岩储层地质建模方法[J]. 天然气工业, 2012, 32(10): 10-13. ZHAO Libin, SHI Shi, XIAO Xiangjiao, et al. Geologic modeling of fractured and porous sandstone reservoirs in the Keshen-2 gas pool of the Kuqa Depression, Tarim Basin[J]. Natural Gas Industry, 2012, 32(10): 10-13. DOI:10.3787/j.issn.1000-0976.2012.10.003 |