复杂地层钻井液漏失诊断技术系统构建.pdf
第 3 2卷 第 6期 2 0 1 5 年1 1月 钻井液与 完井液 DRI LL I NG F LUI D COMP LET I ON F LUI D Vb I . 3 2 NO. 6 No v .201 5 d o i 1 0 . 3 6 9 6 . i s s n . 1 0 0 1 5 6 2 0 . 2 0 1 5 。 0 6 . 0 2 3 复杂地层钻井液漏失诊断技术系统构建 李松 , 康毅力 , 李大奇 , 张杜杰 , 周长林 1 . 中国石油西南油气田公司工程技术研究院,四川广汉 ; 2 . 油气藏地质及开发工程国家重点实验室 西南石油大学,成都 ; 3 . 中国石化石油工程技术研究院,北京 李松等 . 复杂地层钻井液漏失诊断技术系统构建 [ J ] . 钻井液与完井液,2 0 1 5 ,3 2 6 8 9 . 9 5 . 摘要井漏不仅是最严重的储层损害方式,而且是钻井工程中长期悬而未决的重大理论和技术难题。诊断并有效 控制井漏,要从根本上认识并准确描述井漏三要素 位置、类型及强度。漏失机理及类型诊断是制定科学合理的漏失 控制技术的前提。综合利用钻前、随钻、钻后信息资料,描述和表征漏失层性质及参数,进行潜在漏失层预测。建立 了钻井液漏失诊断技术系统框架,提 出了漏失诊断具体方法,综合室内实验和数值模拟开展裂缝、孔洞的应力敏感性 和裂缝传播机制研究,预测漏失通道变形程度和漏失强度 ,为优选堵漏材料提供理论依据 ; 基于钻时、岩屑、钻井液 等工程参数对井漏的异常特征响应来识别井漏,利用实时录井参数来监测并描述井漏状态,提前预测井漏发展趋势 ; 建立了基于漏失发生机理的漏失压力模型,从漏失压力的角度诊断了漏失类型。针对井漏演化过程的认识,初步建立 了井漏诊断技术框架,为漏失控制技术提供理论支持。 关键词 漏失诊断 ; 应力敏感 ; 漏失压力 ;录井 ; 潜在漏失层 中图分类 号 T E 2 5 8 文献标识 码 A 文章编号 1 0 0 1 . 5 6 2 0 2 0 1 5 0 6 . 0 0 8 9 0 7 0 引言 国内外研究大都关注堵漏材料及堵漏技术等能够 直接解决现场井漏复杂技术难题,尚未在定量认识井 漏及漏失诊断理论方面给予同等重视,因而漏失控制 技术依赖着工程经验,而缺少足够明确漏失性质的理 论依据。S h a r 、D y k e 、B e d a 等对比分析现场漏失曲 线 的特征来 区分井漏类型的可能 [ 1 - 3 ] 。L i e t a r d通过实 时钻井液漏失记录研究了天然裂缝性储层中裂缝宽度 的随钻测量, 对于解决实际漏失问题是值得借鉴的 [4 】 。 Ma j i d i 等建立了H . B钻井流体在裂缝 中漏失的流动 模 型 ,得 出钻井液流变特性 动切力 、流性指数等 影 响裂缝性储层漏失 [ 5 】 。金衍等 2 0 0 8建立 了漏失 层位模糊识别钻前风险预测模型 ,并在预测井漏层位 风险程度中取得了很好的效果 [6 ] 。李松等提出了 “ 先 兆漏失”的概念 ,认为漏失通道变化 、漏失动态实时 监测是漏失诊断的关键 】 。沈海超等 2 0 1 3 提出了 基于井漏发生力学机理的诊断及处理技术路线 ,明确 划分潜在漏失层的纵向分布剖面 隋 】 。L i S o n g 等 [9 建 立了H . B钻井流体在一维径向可变形裂缝中的流动 模型,探讨了钻井流体的流变特性、裂缝产状、裂缝 变形及其性质等因素对井漏参数的影响,揭示了变形 裂缝中的钻井液漏失流动行为。 漏失诊断不仅要求在钻前分析尽可能多的地质和 邻井钻井资料,通过地质、力学、钻井、录井等多学 科融合 , 提取地质等特征参数, 从井漏发生因素切入, 预测潜在漏失层 , 更要在漏失过程中实时分析异常工 程参数及现象,定量分析漏失过程,诊断识别漏失类 型 、预测漏失发展趋势 ,科学地优化漏失控制技术。 1 潜在漏失层 预测技术 分析邻井地质、完钻井及测井资料,异常的测井 数据可明确反映出漏失层特征 ,预测 同区块潜在漏失 层 ,从而指导新建井提前预防或避免井漏。 基金项目 国家重点基础研究发展规划 “ 9 7 3 ” 计划 项 目课题 “ 深井复杂地层漏失与井壁失稳机理及预测” 2 0 1 0 C B 2 2 6 7 0 5 和国家科技重大专项 “ 深井复杂地层漏失与井壁失稳机理及预测” 2 0 1 1 Z X 0 5 0 0 5 . 0 0 6 . 0 0 8 H Z 联合资助。 第一作者简介孝松, 博士, 1 9 8 3年生, 2 0 1 4年获西南石油大学油气井工程专业博士学位, 现从事井壁稳定和储层改造方面的研究工作。 地址 四 川省广汉市广东路东一段2号工程技术研究院 ;邮政编码 6 1 8 3 0 0;电话 0 8 3 8 5 1 5 1 4 3 4; E ma i l l i s o n g 0 3 0 1 1 6 4 0 1 6 3 . c o m。 通讯作者 康毅力,电话 0 2 88 3 0 3 2 1 1 8; E - ma i l c wc t k y l v i p . s i n a .t o m。 学兔兔 w w w .x u e t u t u .c o m 9 0 钴井液与 完井液 2 0 1 5年 1 1月 1 . 1 漏失层测井响应 1 . 1 . 1 常规测井响应 识别漏失层 的常规测井技术包括 电阻率测井 、 声波测井和放射性测井技术。裂缝 一 孔洞型储层的常 规测井响应特征为 伽马值较低,当大缝、较大孔洞 中填充泥质时伽马值较高 ; 双侧向测井值较低,且有 正差异 ; 密度较低、中子、声波时差大 ; 小溶孔、溶 洞自然伽马值低一中等,深浅双侧向值差异不明显 ; 中子、密度声波测井曲线变化较大。对于洞径较大的 洞穴型储层,其常规测井响应表现为 深浅双侧向测 井有较大幅度正差异 ; 声波幅度产生严重衰减,声波 时差明显增大 ; 自 然伽马值高,但视溶洞是否被充填 有差异 ; 中子密度异常增大 ; 密度值较低。 1 . 1 . 2 成像测井响应特征 井壁成像测井对复杂孔隙结构如裂缝、溶孔、溶 洞、层理、井壁坍塌等进行直观而精细地描述。成像 测井包含的地层信息是通过成像图上的颜色、形态等 反映出来。其颜色可分为浅色、暗色及杂色 3 类,凭 借成像图上颜色的深浅可对岩性、裂缝等进行识别。 1 成像测井识别裂缝。裂缝在成像测井资料上 表现为缝洞处声脉冲回波幅度小、传播时间长 ; 幅度 图像 中异常面积表征缝、洞的尺寸 ,正弦线参数表征 裂缝的产状; 高导缝在F MI 图像上表现为深色 黑色 的正弦曲线。 2 成像测井识别孔洞。溶洞与溶洞周围岩块反 射情况的差异。溶洞造成传播时间的增加,因而在图 像上很易识别。溶洞在 F MI 图像上颜色发黑,形状 像小砾石,多为分散的星点状或串珠状。 1 . 2 潜在漏失层概率预测 运用工程模糊数学理论 ,建立漏失概率评价数学 模型 ,结合影响漏失通道的岩性、地质 、邻井资料及 对漏失通道较敏感的地震屙陛, 确定隶属度函数,进 而用九标度法和均方根法进行权重计算及分配,对地 层漏失概率进行综合评判 [ 10 ] o该方法集合其它漏失 预测方法的优势综合分析并定量求解,为潜在钻井液 漏失识别提供了新思路。 1 . 3 潜在漏失层的裂缝宽度变化预测 基于地层岩石裂缝多尺度性质 ,通过室内开展岩 心应力敏感性实验 ,分析地层微裂缝在工程扰动下 的 宽度动态变化。微裂缝不能代表岩体中大型裂缝的宽 度变化 ,因此应用已获取的地质力学参数 ,数值模拟 大型裂缝的宽度变化,能预测大型裂缝的宽度变化。 图 1 为室内实验模拟裂缝岩样在不同压差下裂 缝宽度与沿裂缝长度变化情况 ,改变裂缝岩样来模拟 各种复杂地层条件, 分析实验结果可获知地层微裂缝 宽度的动态变化范围。数值模拟地层岩体中大型裂缝 的宽度动态变化,裂缝宽度变化随井筒压差、裂缝长 度 、溶洞发育等 因素的影 响而大 幅增加 图 2 。给 定工程参数,便可模拟预测出裂缝宽度变化范围,从 而有针对性地随钻添加与裂缝动态宽度相匹配的堵漏 材料,及时封堵因应力敏感而扩展的动态裂缝,阻止 扩展性漏失的发生。该技术方法应用于现场试验,明 显减少了井漏事故的发生,同时也降低了井壁坍塌程 度 ,取得 了良好 的预期效果 [1 1 -1 9 ] 逞 褂 憾 有效应力/ MP a 图 1 变有效应力下不同裂缝岩样的缝宽变化 l _3 1 . 2 g 1 . 1 1 . 0 O . 5 0 .4 0 1 2 5 2 5 0 3 7 5 5 0 o 6 2 5 7 5 0 裂缝长度/ ram 图 2 不同洞径缝洞岩体裂缝宽度变化模拟 2 井漏特征诊 断技 术 2 . 1 钻时录井对井漏的响应 采用综合录井仪自 动记录钻时,根据仪器显示的 井深突变和钻时异常变化 ,可对地下的裂缝 、溶洞做 出准确判断。大的裂缝或溶洞甚至可能造成井深突变 几米 ,发生放空现象 ,而相对应的钻时则相对较小或 接近于零 。由于所钻地层存在较大尺寸的裂缝 、孔洞 网络系统,导致钻时突然降低,在钻时曲线上可明显 识别出漏失放空井段。 学兔兔 w w w .x u e t u t u .c o m 学兔兔 w w w .x u e t u t u .c o m 学兔兔 w w w .x u e t u t u .c o m 学兔兔 w w w .x u e t u t u .c o m 钻井液与 完 井液 2 0 1 5 年 1 1 月 图 1 1 a 为裂缝线性和指数变形下的漏失速率变 化情况,分析可知线性变形对应的漏失速率略小于同 等条件下指数变形对应的漏失速率, 但差别不明显, 漏失初始及结束阶段大体一致。图 1 1 b 为裂缝线性 和指数变形下的累积漏失量变化情况 , 从图中可看出, 指数变形对应的累积漏失量明显大于同等条件下的线 性变形对应 的累积漏失量。分析认为 ,裂缝线性变形 使裂缝动态宽度小于指数变形的裂缝动态宽度,从而 线性变形的漏失速率及累积漏失量均小于指数变形对 应的漏失速率及累积漏失量。 喜 0 . 6 0 . 5 嘲l 0 .4 o 3 0 . 2 躐O .1 0 . O 图 l l 不同裂缝变形行为下一维径向裂缝的钻井液漏失行为 5 井漏诊断技术综合框架 漏失诊断技术需充分利用地质、工程资料, 在易 漏失层段中钻进时进行详细的数据监测,针对各阶段 的漏失情况进行诊断分析,初步建立了漏失诊断技术 框架 图 1 2 。 观 地质 因 观工 程因 I 监 测技 术l l 1 . 。 . 南 七 I井 底 压 力 过 裂 缝 宽 度 预 测 分 析 l l 随 钻 录 井 监 测I 层 缝 层 J 上 .. .L . . . L . . . . L I . - . L 岩 及 压 哥 _ L 告 应 数 钻 岩 泥 地 性 孔 力 环 力 值 时 屑 浆 面 及 洞 预 屑 空 动 破 反 测 密 压 压 敏 模 录 返 习 压 碎 育 及 度 耗 力 感 拟 井 出 井 力 程 程 分 增 增 过 实 预 参 数 参 参 大 大 大 验 测 数 据 数 数 [ 厂 广 ’ 一 T _ _ _ _ 一 1 - - ●● -_ _ 一 1 一 一 }井 漏 潜 伏 期} _ _ 井 漏 诱 发 期 卜 噼漏 发 展 期 l } . { 随 钻 漏 失 实 时 监 测I l 判 别 并漏 发生 类型 l 采取 漏失 控 制 技术 图 1 2 井漏诊断技术综合框架 通过对钻井工程、钻井液及气体等参数等实时监 测 ,预报井漏异常情况 ,及时发现因地面或地下异常 而导致相关参数发生的变化,是诊断井漏趋势现有的 一 种技术手段。 6 结论 及 建议 1 . 钻井液漏失诊断技术是一项系统工程 ,应明确 漏失发生机理,识别潜在漏失层,探索应力 一 渗流作 用下的缝洞性岩石开裂、裂缝延伸机理,构建漏失诊 断技术框架。 2 . 结合地质、邻井钻井测井等资料,分析其对漏 失层的异常响应,建立漏失概率评价数学模型,对地 层漏失概率进行综合评判,实现钻前预测同区块潜在 漏失层。 3 . 综合室内应力敏感实验和数值模拟开展多尺度 裂缝 的变形规律 ,从而预测漏失通道变形程度 ,优选 堵漏材料, 该方法经现场试验,已取得显著的防漏堵 漏效果 。 4 . 初步建立了井漏诊断技术框架,分析钻时、岩 屑, 钻井液异常特征及漏失压力诊断来识别漏失类型 , 描述井漏过程,为优化漏失控制技术提供理论依据。 5 . 漏失控制是一项系统工程,加强潜在漏失层预 测及漏失通道变形规律研究 ,重视井漏数据采集及钻 井液漏失行为描述 ,完善钻井液漏失动力学理论。 参 考 文 献 [ 1 】 S h a r e r D M ,L o e p p k e G E,Gl o wk a D A,e t a 1 . An e v a l u a t i o n o f fl o wme t e r s f o r t h e d e t e c t i o n o f k i c k a n d l o s t c i r c u l a t i o n d u r i n g d r i l l i n g [ C ] . I A DC / S P E 2 3 9 3 5 . [ 2 ] B e d a G,C a r u g o C.Us e o f mu d mi c r o l o s s a n a l y s i s wh i l e d r i l l i n g t o i mp r o v e the f o r ma t i o n e v a l u a t i o n i n fra c t u r e d r e s e r v o i r s [ C 】 . S P E 7 1 7 3 7 . 【 3 ] Dy k e C G,Wu B,Mi l t o n T a y l e r D. Ad v a n c e s i n c h a r a c t e r i z i n g n a t u r a l - fra c t u r e p e r m e a b i l i t y fro m mu d -- l o g d a t a [ C ] . S P E 2 5 0 2 2 . [ 4 ] O l i v e r L i e t a r d . F r a c t u r e w i d th L WD a n d d r i l l i n g mu d / L C M s e l e c t i o n g u i d e l i n e s [ C 】 . S P E 3 6 8 3 2 . [ 5 ] Ma j i d i R,Mi s k a S Z,Yu M. Qu a n t i t a t i v e a n a l y s i s o f mu d l o s s e s i n n a t u r a l l y f r a c t ure d r e s e rvo i r st h e e f f e c t o f r h e o l o g y [ C ] . S P E 1 1 4 1 3 0 . [ 6 ] 金衍,卢运虎,李再钧 . 一种井漏层位钻前风险预测新 方法 [ J ] .石油钻采工艺,2 0 0 8 ,3 0 3 2 4 . 2 8 . J i n Ya n , L u Y u n h u , L i Z a ij u n . A n e w me t h o d p r e d i c t i n g p r e - d r i l l i n g r i s k i n l o s t c i r c u l a t i o n f o r ma t i o n [ J ] . O i l Dr i l l i n g &P r o d u c t i o n T e c h n o l o g y ,2 0 0 8 ,3 0 3 2 4 - 2 8 . 学兔兔 w w w .x u e t u t u .c o m 第 3 2卷 第 6期 李松等 复杂地层钻 井液漏失诊断技术 系统构建 9 5 [ 7 】 李松, 康毅力 , 刘修善, 等 . 先兆漏失诊断钻井液漏 失控制的一个关键 [ J 】 . 钻井液与完井液, 2 0 1 1 , 2 8 2 8O- 8 3. Li S o n g,Ka n g Yi l i ,Li u Xi u s h a n,e t a 1 . P o r t e n t d i a g n o s t i c o f mu d L o s s - A k e y t o mu d l o s s c o n t r o l l i n g [ J ] .Dr i l l i n g F l u i d&C o m p l e t i o nF l u i d ,2 0 1 1 ,2 8 2 8 0 8 3 . [ 8 】 沈海超 ,胡晓庆,李桂芝.破碎性地层漏失力学机理及 井漏诊断与处理思路 [ J ] . 钻井液与完井液, 2 0 1 3 , 3 0 1 8 5 . 8 8 . S h e n Ha i c h a o,Hu X i a o q i n g,Li Gu i z h i . Re s e a r c h o n me c h a n i c a l me c h a n i s m a n d t r e a t i n g me t h o d o f mu d l o s s i n fr a g i l e f o r ma t i o n [ J ] . Dr i l l i n g F l u i d& C o m p l e t i o n Fl u i d 。 2 0 1 3 ,3 0 1 8 5 . 8 8 . [ 9 ] L i S o n g ,Ka n g Yi l i ,L i D a q i ,e t a 1 .Mo d e l i n g H e r s c h e l B u l k e l y d r i l l i n g fl u i d fl o w i n a v a r i a b l e r a d i a l f r a c t u r e [ J ] . J o u r n a l o f P o r o u s Me d i a ,2 0 1 4 ,1 7 3 2 3 9 - 2 5 4 . [ 1 0 ]康毅力,王海涛 ,游利军,等.基于层次分析法的地层 钻井液漏失概率判定 [ J ] .西南石油大学学报 自然科学 版 2 0 1 3 ,3 5 4 1 8 0 . 1 8 6 . Ka n g Yi l i ,Wa n g Ha i t a o , Yo u L ij u n , e t a 1 .P r o b a b i l i t y d e t e r mi n a t i o n f o r l o s s c i r c u l a t i o n o f d r i l l i n g fl u i d s b a s e d o n a n a l y t i c h i e r a r c h y p r o c e s s[ J 】 . J o u r n a l o f S o u t h w e s t Pe t r o l e u m Un i v e r s i t y S c i e n c e& T e c h n o l o g y Ed i t i o n , 2 0 1 3 ,3 5 4 1 8 0 . 1 8 6 . [ 1 1 ]兰林 ,康毅力,陈一健,等.储层应力敏感性评价实验 方法与评价指标探讨 [ J ] _ 钻井液与完井液, 2 0 0 5 , 2 2 3 1 - 4. L a n L i n ,Ka n g Yi l i ,C h e n Yi j i a n,e t a 1 .Di s c u s s i o n o n e v a l u a t i o n m e t h o d s f o r s t r e s s s e ns i t i v i t i e s o f l o w p e rm e a b i l i t y a n d t i g h t s and s t o n e r e s e r v o i r s [ J ] . D r i l l i n g F l u i d &C o m p l e t i o n F l u i d ,2 0 0 5 ,2 2 3 1 - 4 . [ 1 2 ]李松 ,康毅力,李大奇 ,等 . 考虑孔洞沟通性的碳酸盐 岩储层应力敏感性实验 [ J ] _ 重庆大学学报, 2 0 1 2 , 3 5 5 l 1 5 . 1 2 0. L i S o n g,Ka n g Yi l i ,L i Da q i ,e t a 1 . L a b o r a t o r y i n v e s t i g a t i o n o n s t r e s s s e n s i t i v i ty o f c a r b o n a t e r e s e rvo i r c o n s i d e rin g v u g g y c o n n e c t i v i ty[ J ] .J o u r n a l o f C h o n g q i n g U n i v e r s i t y Na t u r a l S c i e n c e E d i t i o n ,2 0 1 2 ,3 5 5 1 1 5 - 1 2 0 . [ 1 3 】K a n g Yi l i ,Yo u L i j u n ,X u Xi n g h u a ,e t a 1 . P r e v e n t i o n o f f o r ma t i o n d a ma g e i n d u c e d b y mu d l o s t i n d e e p f r a c t u r e d t i g h t g a s r e s e rvo i r i n we s t e rn S i c h u a n B a s i n [ J ] . J C a n P e t T e c h n o 1 . 2 0 1 2 ,5 1 1 4 6 . 5 1 . [ 1 4 】雷鸣,瞿佳,康毅力 , 等 . 川东北裂缝性碳酸盐岩气层钻 井完井保护技术 [ J ] . 断块油气田, 2 0 1 1 ,1 8 6 7 8 3 . 7 8 6 . L e i Mi n g,Qu J i a ,Kang Yi l i ,e t a 1 . P r o t e c t i o n t e c h n o l o g y o f d r i l l i n g a n d we l l c o mp l e t i o n o f f r a c t u r e d c a r b o n a t e g a s r e s e r v o i r i n n o r t h e a s t e r n S i c h u an B a s i n [ J ] . F a u l t - B l o c k O i l & G a s F i e M, 2 0 1 1 ,1 8 6 7 8 3 . 7 8 6 . [ 1 5 】闫丰明,康毅力,孙凯,等.裂缝 - 孔洞型碳酸盐岩储 层暂堵性堵漏机理研究 [ J ] . 石油钻探技术, 2 0 1 1 , 3 9 2 81 . 85. Ya n F e n g mi n g,Kan g Yi l i ,S u n Ka i ,e t a 1 . Me c h a n i s m o f t e mp o r a r y s e a l i n g f o r f r a c t u r e d - v u g g y c a r b o n a t e r e s e r v o i r [ J ] . P e t r o l e u m Dr i l l i n g T e c h n i q u e s ,2 0 1 1 ,3 9 2 8 1 8 5 . [ 1 6 】许成元 ,康毅力 ,游利军,等.裂缝性储层渗透率返排 恢复率 的影响因素 [ J ] . 石油钻探技术,2 0 1 2 ,4 0 6 1 7 . 21 . Xu C h e n g y u a n ,K a n g Yi l i , Yo u L i j u n ,e t a 1 . I n fl u e n t i a l f a c t o r s o n p e r m e a bi l i t y r e c o v e r y d u r i n g f l o w b a c k o f f r a c t u r e d r e s e r v o i r s [ J ] . P e t r o l e u m Dr i l l i n g T e c h n i q u e s , 2 0 1 2 ,4 0 6 1 7 . 2 1 . [ 1 7 】李大奇,康毅力,刘修善,等 . 基于漏失机理 的碳酸 盐 岩 地层 漏 失 压 力模 型 [ J ] . 石 油 学报 ,2 0 1 1 ,3 2 5 9 0 0 . 9 0 4 . Li Da q i , Ka n g Yi l i ,Li u Xi u s h a n, e t a 1 . T h e l o s t c i r c u l a t i o n p r e s s u r e o f c a r b o n a t e f o r ma t i o n s o n t h e b a s i s o f l e a k a g e me c h a n i s ms [ J 1 . Ac t a P e t r o l e i S i n i c a,2 0 1 1 ,3 2 5 9 0 0 . 9 0 4 . [ 1 8 】L i S o n g ,Ka n g Yi l i ,Yo u L i j u n ,e t a 1 . E x p e ri me n t a l a n d nu m e ric a l i nv e s t i g a t i on o f mul t i - s c a l e f r a c t ur e de f o r ma t i on i n fr a c t u r e d - v u g g y c ar b o n a t e r e s e rvo i r s [ J ] . T h e A r a b i a n J o u r n a l _厂 0 r S c i e n c e a n d E n g i n e e r i n g,2 0 1 4 ,3 9 5 4 2 4 1 - 4 2 4 9 . [ 1 9 ]李大奇,康毅力,曾义金,等.缝洞型储层缝宽动态变 化及其对钻井液漏失的影响 [ J ] .中国石油大学学报 自 然科学版 ,2 0 1 1 ,3 5 5 7 6 . 8 1 . L i Da q i ,Ka n g Yi l i ,Z e n g Yij i n , e t a 1 . D y n a mi c v a ri a t i o n o f f r a c t u r e wi d t h a n d i t s e f f e c t s o n d r i l l i n g flu i d l o s t c i r c u l a t i o n i n fr a c t u r e d v u g g y r e s e rvo i r s [ J ] .J o u r n a l o f C h i n a U n i v e r s i ty o f P e t r o l e u m E d i t i o n of Na t u r a l S c i e n c e , 2 0 1 1 ,3 1 5 7 6 . 8 1 . [ 2 0 ]刘向君,罗平亚.岩石力学与石油工程 [ M] .北京 石 油工业出版社,2 0 0 41 1 8 . 1 2 2 . L i u Xi a n g j u n, L u o P i n g y a . R o c k me c h ani c s and p e t r o l e um e n g i n e e r i n g[ M】 . B e i j i n g P e tr o l e um I n d u s t r y P r e s s , 2 0 0 4 1 77 . 1 8 2. 【 2 1 ]L a v r o v ,J o h a n T r o n v o l 1 .Mo d e l i n g mu d l o s s i n fra c t u r e d f o r ma t i o n s [ C ] . S P E 8 8 7 0 0 . [ 2 2 ]B r u e l D,C a c a s MC,L e d o u x E ,e t a 1 .Mo d e l i n g s t o r a g e b e h a v i o r i n a fr a c t u r e d r o c k ma s s [ J ] .J o u r n a l o f H y d r o l o g y , 1 9 9 4, 1 6 2 2 6 7 - 2 7 8 . 收稿 日期2 0 1 5 . 5 . 1 3 ;HG F 1 5 0 5 N4 ;编辑王小娜 学兔兔 w w w .x u e t u t u .c o m