基于地震波阻抗的预探井随钻井壁稳定预测.pdf

3 9 0 2 0 1 5年 6月 石油勘探与开发 P ET ROL EUM E XPL ORATI ON AND DEVEL 0P M E NT V b 1 ． 4 2 NO ．3 文章编号 1 0 0 0 ． 0 7 4 7 2 0 1 5 0 3 ． 0 3 9 0 ． 0 6 DO I 1 0 ． 1 1 6 9 8 ／ P E D． 2 0 1 5 ． 0 3 ． 1 8 基于地震波阻抗的预探井随钻井壁稳定预测 吴超 ，刘建华 ，张东清 ，陈小峰 ，赵伟杰 2 1 ．中国石油化工股份有限公 司石油工程技术研究院；2 ．中国石油大学 北京 石油天然气工程学院 基金项目国家科技重大专项 “ 海相油气井井筒环境监测技术” 2 0 1 1 Z X0 5 0 0 5 ． 0 0 6 0 0 1 摘要 常规井壁稳定预测方法应用于预探井时，由于资料缺乏及使用条件 限制、计算参数误差、运算过程复杂等因 素会严重影响其预测效果，针对这些问题提 出了一种适用于预探井的随钻井壁稳定预测方法。通过考察波阻抗与孔 隙压力、构造地应力、岩石强度之间的定量关系，建立包含波阻抗与井壁稳定力学参数的非线性模型，在此基础上 利用神经网络算法识别 波阻抗与地层三压力 孔隙压力、坍塌压力与破裂压力 之 间的映射关系。通过神经网络分 层建模及录井资料实时分析 ，实际钻进 中利用地震反演波阻抗数据随钻预测钻头前方的井壁稳定性 ，及时优化钻井 液密度以控制井壁稳定性。现场预探井的应用结果表明，新方法较之常规技术有效提高 了适应性与运算速度，操作 流程更为便捷，预测精度能满足工程需求。 图5 表 1 参 2 1 关键词 预探井；井壁稳定性；随钻预测；地震波阻抗；神经网络算法 中图分类号 T E 2 2 文献标识码A A pr e di c t i o n o f bo r e ho l e s t a b i l i t y whi l e d r i l l i ng p r e l i mi na r y pr o s pe c t i ng we l l s ba s e d o n s e i s m i c i m p e da nc e Wu C h a o ， L i u J i a n h u a ， Z h a n g D o n g q i n g ， C h e n X i a o f e n g ， Z h a o We 0 i e 1 ． S I NO P E CR e s e a r c h I n s t i t u t e o fP e t r o l e u m E n g i n e e r i n g ， B e ij i n g 1 0 0 1 0 1 ， C h i n a ； 2 ． F a c u l ty of P e t r o l e u m E n g i n e e n i n g ， C h i n a U n i v e r s i t y of P e t r o l e u m， Be ij i n g 1 0 2 2 4 9 ， C h i n a Ab s t r a c t Wh e n t h e c o n v e n t i o n a l me t h o d o f p r e d i c t i n g b o r e h o l e s t a b i l i t y i s a p p l i e d t o t h e p r e l i mi n a r y p r o s p e c t i n g we l l s ， t h e p r e c i s i o n o f p r e d i c t i o n i s o f t e n a ff e c t e d s e v e r e l y b y the l a c k o f d a t a ， l i mi t o f a p p l y i n g c o n d i t i o n s ， e r r o r o f c a l c u l a t e d p a r a me t e r s a n d c o mp l e x i t y o f o p e r a t i n g c o u r s e s ．T h i s p a p e r p r e s e n t s a me t h o d o f p r e d i c t in g b o r e h o l e s t a b i l i ty wh i l e d r i l l i n g p r e l i mi n a r y p r o s p e c t i n g we l l s and e s tabl i s h e s a n o nl i n e a r mod e l i nc l u d i ng wa v e i mpe d a n c e a nd bo r e ho l e s t a bi l i ty me c ha ni c a l pa r a me t e r s b y i n ve s t i g a t i ng t h e q u an t i tat i ve r e l a t i o n s h i p s b e t we e n wa v e i mp e d a n c e a n d p o r e p r e s s u r e ， i n s i t u s t r e s s a s we l l a s r o c k s e n g t h ． Ba s e d o n t h e n o n l i n e ar mo d e l ， n e ur a l n e t w o r k a l g o ri t h m i s u s e d t o i d e n t i f y the r e l a t i o n s h i p b e t w e e n the w a v e i mp e d anc e and t h e t h r e e f o r ma ti o n p r e s s ure s p o r e p r e s s ure ， fr a c t u r e p r e s s u r e a n d c o l l a p s e p r e s s ure ． T h r o u g h t h e e s t a b l i s h me n t o f t h e mo d e l b y l a y e r e d n e ura l n e t wo r k a n d t h e t ime l y a n a l y s i s o f l o g g i n g d a ta， t h e we l l b o r e s t a b i l i ty b e f o r e t h e b i t c a n b e p r e d i c t e d wh i l e d r i l l i n g b y u s i n g t h e d a t a o f s e i s mi c wa v e i mp e d a n c e ． F i e l d a p p l i c a t i o n i n p r e l i mina r y p r o s p e c ti n g we l l s s h o ws t h a t t h e n e w me t h o d h a s h i g h e r a d a p t a b i l i ty, f a s t e r c a l c u l a t i o n s p e e d ， a n d s i mp l e r o pe r a t i o n al p r o c e d ure tha n c o n ve nt i o n a l me t ho d s ， a nd i t s pr e d i c t i o n a c c ur a c y C an me e t t h e r e qu i r e me n t of e n gi ne e rin g． Ke y wo r d s p r e l i mi n a r y p r o s p e c t i n g we l l ； b o r e h o l e s t a b i l i ty； p r e dic t i o n wh i l e d r i l l i n g ； wa v e i mp e d anc e ； n e u r a l n e t wo r k O引言 钻井工程中由于井壁失稳引起的复杂故障会造成 严重的钻井质量与安全事故 ，准确描述井壁围岩力学 环境 、深入探索井壁失稳机理 、建立合理 的井壁稳定 计算模型 ，是实现安全 、优质 、高效钻井 的关键 。目 前井壁稳定钻后定量检测方法 已基本成熟[ 1 】 ，它通过 钻井与测井资料 的综合分析 ，确定井壁失稳的主导机 理 ，计算并建立地层孔隙压力 、坍塌压力及破裂压力 剖面，指导后续钻井 。但从实际工程角度而言 ，要求 对井壁失稳问题尽可能做到提前预知 、及早发现、尽 快处理 ，所以井壁稳定性预测研究意义更为重大。当 前井壁稳定预测方法总体包括两类 ，首先提 出的是钻 前预测方法 ，它利用地震特征参数分析 、地震反演等 方法预测得到待钻井全井段的声波速度等参数 ，利用 其预测井壁稳定性并用于钻井设计[ 6 ] 。对探井而言 ， 由于钻前对其地质情况 了解有限 ，钻前预测精度常不 够理想 ，需要在实际钻进中根据实钻情况对钻前预测 模型进行实时调整与更新 ，井壁稳定随钻预测方法 因 此发展起来[ 1 0 - 1 2 ] 。这类方法基于地震属性与测井数据 之间的相关性 ，综合利用地震 、测井 、录井及钻井资 料 ，对钻遇的每一层系的井壁稳定性进行精细分析评 价 ，实钻 中及时修正钻前预测数据。常规随钻预测方 法虽然有效提高了预测精度及效率 ，但也存在一些问 2 0 1 5 年 6 月 吴超 等 基于地震波阻抗的预探井 随钻井壁稳定预测 3 9 1 题 对于资料相对缺乏 的预探井 ，其井壁稳定计算模 型参数存在较大误差 ，多步运算过程中的误差 累计会 严重降低最终预测结果的准确度 ；另一方面对于待预 测 的预探井尤其是区域探井 ，其与参照井可能位于不 同的物探测 网， 这也会影响常规预测方法的应用效果。 针对上述情况 ，本文提 出了一种新的井壁稳定 随 钻预测方法 ，该方法根据波阻抗与岩石力学参数、孔 隙压力、地应力等井壁稳定参数之间存在的映射关系 ， 综合利用地震反演波阻抗 、智能计算工具及实钻信息 进行井壁稳定随钻预测。 1波阻抗与井壁稳定参数 的关系 1 ． 1 概述 地震反演是 以地震勘探资料为基础 ，以已知地质 特征 以及钻井、测井资料为约束 ，对地下岩层空间的 岩石物理参数进行成像 的过程 。地震资料 中包含丰富 的岩性 、构造 、物性信息 ，经过地震反演可以把地震 信息转换为岩石物理信息 ，建立起适用于勘探开发的 地质、力学和储集层模型。 按所使用的地震资料来划分 ，地震反演主要包括 叠前反演和叠后反演两大类 。叠前反演主要包括基于 旅行时的层析成像技术和基于振幅的 AVO分析技术 ； 叠后反演主要指基于振幅的波阻抗及其他岩石物理参 数反演技术。其 中波阻抗反演技术发展得最为成熟 ， 具有严密 的物理和数学理论基础[ 1 3 - 1 4 ] ，较好地将反演 精度与运算效率结合起来 ，反演结果集地震剖面的横 向连续性 、测井资料垂向高分辨率以及地质剖面的直 观性于一体 ，综合来看最适合钻井工程师使用。 波阻抗是重要 的地层特征参数 ，为地震波在地层 中传播速度和岩石密度的乘积。以目前最常用的叠后 纵波阻抗反演为例 ，纵波阻抗公式为 I p p y p 1 由公式 1 可见 ，地层波阻抗与岩石的性质 ，包 括岩性、岩石物理、岩石力学、埋藏深度、地质年代、 构造特征 、含油气性等均关系密切。因此 ，特定 的波 阻抗分布规律包含丰富的地层信息 ，能不同程度地反 映地层的岩石力学特征 ，这为利用波阻抗预测井壁稳 定性提供了依据 。 岩石力学理论表明，导致井壁垮塌和破裂的根本 原因是井筒 围岩发生剪切和拉伸破坏 ， 岩石强度参数 、 地层孔 隙压力、构造地应力是决定井壁稳定性 的关键 地质参数。当以波阻抗数据为基础预测井壁稳定性时， 必须充分考察井壁稳定参数与波阻抗之间的关系 ，以 反映孔隙压力 、构造应力及岩石强度 包括粘聚力 、 内摩擦角 、抗拉强度等 等参数对井壁稳定性 的影响， 在此基础上优选井壁稳定预测模型。 1 ． 2 波阻 抗与地 层孑 L 隙压 力之 间关系 岩石物理实验表 明，纵波速度与有效应力之 间存 在如下关系㈣ V p S o 蜀 P B 2 √ h f 1 0 0 o “v P p。 一 e l 2 根据岩石物理实验 ，泥质含量和纵 、横波速度之 间存在以下经验关系n 】 a 2 v p a 3 v 3 Y s √ 4 综合考察 以上各式可见 ，地层孔隙压力与地震波 速度、密度之间存在一定函数关系，波阻抗是地震波速 度与密度的乘积， 且地震波速度与密度之间也存在一定 相关性 。考虑到地震反演波阻抗相 比分别反演速度和 密度可信度更高 ，因此建立孑 L 隙压力与波阻抗之间的 映射关系模型更适于预测孔隙压力及井壁稳定性 ，即 f l P ， ， 0 5 1 ． 3 波 阻抗 与构 造应 力之 间关系 勘探开发实践表 明，地震波在岩层 中的传播速度 整体上随地质过程中构造作用力 的增大而增大 ，即波 阻抗 与构造运动关系密切。 目前 国内外已提出多种构 造地应力计算模型 ，其适应性取决于 目标探 区的地质 情况 。以其 中使用最为广泛的黄氏模型为例 “ ] 篇l～ g p p㈩ 2} 1 一 篇1 l 7 I 一 l一 I 由弹性动力学理论可知 ，岩石弹性参数与地震波 速度及 密度 关 系密切 [ 1 7 ] E 1 0 3 v p 一 4 ／ u p 一 8 一 2 ／ 2 一 9 综合考察以上各式可见 ， 构造应力与地震波速度 、 岩石密度之间存在一定的函数关系 ，可建立构造地应 力与波阻抗之间的映射关系 ，即 L M， ， 0 1 0 1 ． 4 波 阻抗 与岩石 强度之 间关 系 通过岩石力学实验数据分析 ，可以建立适用于钻 井工程的多种岩石强度经验计算模型 ，例如岩石抗压 3 9 2 石油勘探 与开发 综合研究 V _0 l | 4 2 NO ． 3 和抗拉强度可 以使用 以下经验公式计算[ “ ] C l C 2 h 1 1 S t ／ C 3 1 2 而基于实验资料统计得到的岩石粘聚力与内摩擦 角的常用计算公式【 为 C - A t1 - 2 ／／ ] 1 1 3 U l 1 4 综合考察以上各式 ，认为岩石强度与波阻抗之间 存在一定 的映射关系，再考虑到岩石强度受围压影响 较大 ，结合上文 中构造应力计算模型 ，可以建立包含 岩石强度、波阻抗等因素的关系模型 ， ， 0 1 5 2井壁稳定随钻预测模型 2 ． 1预 测原理 按照通常的研究思路 ，在通过地震反演得到待钻 井 的波阻抗数据之后 ，分别利用 以上各种理论或经验 公式求取孔隙压力 、构造应力、粘聚力 、内摩擦角 、 抗拉强度等数据 ，再利用剪切及拉伸破坏准则分别求 取岩石坍塌压力与破裂压力 ，通过地层三压力 孔隙 压力 、坍塌压力与破裂压力 综合分析 ，对井壁稳定 性进行预测。但对于预探井而言 ，由于探区前期完钻 井数量少 ，用于统计分析 的测井 、测试 、地质等资料 较为匮乏 ，导致常规预测模型中的大量系数和经验参 数存在较大误差 ；同时 ，常规井壁稳定预测方法计算 步骤繁多 ，致使误差逐步累积 ，最终计算结果的预测 精度也将会受到很大影响。其次 ，由于预探井地质条 件复杂 ，常规方法中复杂的运算流程易导致操作出错 。 此外 ，由于地震属性分析 的需要 ，常规预测方法要求 参照井和待预测井必须位于同一地震测网内，对于预 探井尤其是 区域探井而言，其与 已完钻的参照井常具 有不同的观测系统 ，这也会严重影响常规方法的预测 效果 ，所以需研发新方法用于预探井的井壁稳定预测 。 从另一个角度来看 ，波阻抗和上覆压力分别与孔 隙压力 、构造应力 、岩石强度之间构成隐函数关系 ， 进一步根据井壁稳定力学理论 ，波阻抗和坍塌压力、 破裂压力之间也存在一定的非线性映射关系 ，即 ， ， 口 ， ， ‰ ， P c ， P f 0 1 6 由于钻井地质环境 的复杂性和测量手段的限制 ， 这种函数关 系具有非线性、多解性 、模糊性 、随机性 等特征 ，因此直接确定或应用这种非线性关系存在一 定困难 。可使用智能计算工具处理这类 问题 ，其特点 是无需直接确定非线性映射模型的系数 ，而通过智能 算法 自动识别各类地质参数之间的关联性。 目前在理 论和应用上都 比较成熟的智能算法有神经网络 、遗传 算法、蚁群算法等 ，针对预探井井壁稳定预测问题的 实际要求 ，本文采用神经网络算法实现这一过程 ，即 首先利用完钻井资料进行神经网络学习建模 ，再利用 其预测 目标井的井壁稳定性 。 2 ． 2 操作流 程 2 ． 2 ． 1钻前阶段工作 ①收集 、整理与分析已完钻参照井的测井、钻井 、 地质 、测试等资料 ，对于本区已有完钻井 的预探井 ， 最好选取 同区邻井作为参照井 ；对于同区无完钻井的 区域探井 ，则选取邻近区块具有代表性的完钻探井作 为参照井 。② 以测井数据为基础 ，结合其他资料 ，求 取参照井 的波阻抗、上覆压力 、孑 L 隙压力 、坍塌压力 及破裂压力数据 。需要说 明的是 ，在对区域探井进行 预测时 ，由于资料有限 ，利用参照井密度测井资料建 立 的上覆压力梯度模型也用于待预测井。③利用神经 网络算法 ，分层系建立参照井各层段波阻抗 、上覆压 力 、孔隙压力 、坍塌压力 、破裂压力的网络模型 ，以 供实钻计算使用。④收集待预测井所在工区的地震资 料 ，利用广义线性算法反演得到待预测井的波阻抗数 据 ，以供实钻计算使用。 2 ． 2 ． 2 实钻过程 中的随钻预测工作 ①实时检测分析待预测井的录井资料 ，确定当前 钻遇地层的类型。②从 已建立的分层神经网络模型中 选择对应当前钻遇地层类型的模型。③将待预测井钻 头前方一定深度范围内的波阻抗、上覆压力数据输入 神经网络 ，经过计算可以从 网络输 出待钻地层 的孔 隙 压力 、坍塌压力及破裂压力 ，进而预测待钻地层的井 壁稳定性。在 以上步骤 中，最关键 的是建立参照井 的 神经网络模型并利用其计算待预测井的地层三压力 。 下 面重点 论述 这一过 程 。 神经网络算法种类很多 ，根据井壁稳定预测过程 具有先学习再应用的特点 ，同时随钻预测操作要求学 习速度快 、学习精度高，这种情况下较适宜采用有监 督学习的网络。 有监督学习的神经网络主要包括 B P网 络 、径 向基 函数网络、小波网络等。笔者结合实际问 题对上述几种常用算法进行 了反复试算 比较 ，结果表 明在利用神经 网络建立波阻抗与地层压力之间的非线 性映射关系模型时，运用径 向基函数学习算法能在保 证较高预测精度的同时，还具有较快的计算速度 ，因 此在进行井壁稳定随钻预测时使用径向基函数神经 网 络效率最高 。除了神经 网络 ，笔者在研究中也尝试应 2 0 1 5年 6月 吴超 等 基 于地震波阻抗 的预探井随钻井壁稳 定预测 3 9 5 应用中 ，在钻前进行地震波阻抗反演和神经网络分层 学习建模 ， 实钻中利用建立的神经网络分层模型分层预 测钻头下方待钻地层的井壁稳定性， 可以及时优化钻井 液密度 以控制井壁稳定性 。该方法运算过程简捷 、经 验参数少 ，可 以减少误差 累积 ，较快的运算速度有利 于实现随钻预测 ，在现场应用 中能够达到较高的预测 精度 ，还克服 了传统预测方法对地震资料使用条件的 限制 ， 适用于资料相对缺乏的预探井尤其是 区域探井 。 符号注释 ， D 纵波阻抗 ，1 0 k g ／ IT I S ；1 ； p 纵波速度， k m／ s ；p 岩石密度 ，g ／ c m ； h 泥质含量 ，f ； 上覆地层压力 ，MP a ；P p o 地层孔隙压力 ，MP a ；a 1 ，a 2 ， a 3 ，A1 ，A2 ， Be ， B1 ， B2 ， B3 ， B4 ， C1 ， ， C3 ， VI ， ， Z 3 ，“ 1 ，u 2 ， 1 ， 模型系数 ，无量纲 ；V s 横波速度 ， k m ／ s ；a H ， 水平方向的最大与最小地应力 ，MP a ； 弹性模量 ，MP a ； 泊松比，无量纲 ；a 有效应力系 数 ，无量纲 ； 抗压强度 ，MP a ；S抗拉强度 ，MP a ； C粘聚力 ， MP a ； 内摩擦角 ， 。 ； 构造 应力 ，MP a ； 岩石强度参数 ，MP a ；P 坍塌压力 ， MP a ； P f 破裂压力 ，MP a ； 对应深度待预测的地层 3压力值 ，MP a ；c o 隐层到输出层的连接权值 ，无量纲 ； x 参照井的波阻抗与上覆压力数据组成的输入向量 ； ， 厂一径向基 函数的中心和宽度 ，均为 2维向量 ；ok t 聚类区；a O误差 ， 无量纲 ； D ．『 计算编码 ， 无量纲 ； M 0 | 『 中包含的 数 目，无量纲；i 样本编号 ； H 学习样本对的数 目；， 输 出层节点编号 ； 连接 权 ，无 量纲 ； 基 函数 的个数 ；， 递推 步骤数 ； G 。 上覆地层压力梯度当量钻井液密度 ，g ／ c m ； 井 深，m； f t ， ， ， 表示波阻抗与井壁稳定性参数之间 关系的函数 ，无量纲。 参考文献 Ka r s t a d E， Aa d n o y B S．Op t i mi z a t i o n o f b o r e ho l e s t a bi l i t y u s i n g 3 D s t r e s s o p t i mi z a t i o n s [ R ] ． S P E 9 7 1 4 9 ， 2 0 0 5 ． He mp h i l l T．I n t e g r a t e d m a n a ge me n t o f t he s a f e o pe r a t i n g wi n d o w We l l b o r e s t a b i l i t y i s mo r e t h a n j u s t f l u i d d e n s i t y [ R ] ． S P E 9 4 7 3 2 ， 2 0 0 5 ． Gr e e n wo o d J A， Br e h m A， Va n Oo f t E， e t a 1 ． Ap pl i c a t i o n o f r e a l t i me we l l b o r e s t a b i l i t y mo n i t o r i n g o n a d e e p wa t e r E R D we l l [ R ] ．S P E 9 25 8 8 ， 20 0 5 ． Ne s O M， F j a r E ， T r o n v o l l J ， e t a 1 ． D r i l l i n g t i me r e d u c t i o n t h r o u g h i n t e g r a t e d r o c k me c h a n i c s a n a l y s i s [ R ] ． S P E 9 2 5 3 1 ， 2 0 0 5 ． Br a d f o r d I D R， Al d r e d W A， Co o k J M ， e t a 1 ． W h e n r o c k me c h a n i c s m e t d r i l l i ng Ef f e c t i v e i mp l e me nt a t i o n o f r e a l t i m e we l l b o r e s t a bi l i t y c o n t r o l [ R ] ． S P E 5 9 1 2 1 ， 2 0 0 0 ． 金衍，陈勉，张旭东．钻前井壁稳定预测方法的研究[ J ] ． 石油学 报 ， 2 0 0 1 ， 2 2 3 9 6 9 9 ． J in Yan ，Ch e n M i a n ， Zh a n g Xu d o n g ． S t u d y o n pr e d i c t i o n f o r b o r e h o l e s t a b i l i t y b e f o r e d r i l l in g [ J ] ． A c t a P e t r o l e i S i n i c a ， 2 0 0 1 ， 2 2 3 9 6 9 9 ． 金衍 ，陈勉，杨小奇 ．利用层速度钻前预测安全 泥浆密度 窗 口[ J ] ． 岩石力学与工程学 报， 2 0 0 4 ， 2 3 1 4 2 4 3 0 2 4 3 3 ． J i n Ya n ，Ch e n M i a n ，Ya n g Xi a o q i ．P r e d i c t i o n of s a f e mu d d e n s i t y w i n d o w b y i n t e r v a l v e l o c i t y b e f o r e d r i l l i n g [ J ] ． C h i n e s e J o u r n a l o f R o c k Me c h a n i c s a n d E n g i n e e r i n g ， 2 0 0 4 ， 2 3 1 4 2 4 3 0 2 4 3 3 ． 【 8 ] 金衍，陈勉 ，吴超．探井二开 以下地 层井壁稳定性钻前 预测方法 [ J ] ．石油勘探与开发， 2 0 0 8 ， 3 5 6 7 4 2 7 4 5 ． J i n Ya n ，Ch e n M i a n ，W u Ch a o． Bo r e h o l e s t a b i l i t y p r e d i c t i o n b e f o r e d r i l l i n g for f o r ma t i o n s u n de r t he s e c o n d s e c t i o n i n a n e x p l o r a t i o n we l l [ J ] ． P e t r o l e u m E x p l o r a t i o n a n d De v e l o p me n t ， 2 0 0 8 ， 3 5 6 7 4 2 7 4 5 ． [ 9 ]9 金衍 ，陈勉．利用地震记录钻前 预测井壁稳定性研究 [ J 】 ．石油学 报， 2 0 0 4 ， 2 5 1 8 9 9 2 ． J i n Ya n ，Ch e n M i a n ．P r e d i c t i o n o f b o r e h o l e s t a b i l i t y b y s e i s mi c r e c o r d s [ J ] ． Ac t a P e tr o l e i S i n i c a ， 2 0 0 4 ， 2 5 1 8 9 9 2 ． [ 1 0 】 吴超 ，陈勉，金衍．井壁稳定 实时预测新方法 [ J 1 ．石油勘探与 开 发， 2 0 0 8 ， 3 5 1 8 0 培4 ． W u Ch a o ，Ch e rt M i a n ，J i n Ya h．Re a l t i me p r e d i c t i o n m e t h o d o f b o r e h o l e s t a b i l i t y [ J ] ． P e t r o l e u m E x p l o r a t i o n a n d D e v e l o p me n t ， 2 0 0 8 ， 3 5 1 8 0 8 4 ． [ 1 1 】 吴超，陈勉，金衍 ．基于地震属 性技术 的井壁稳定 随钻预测新方 法【 J 】 _中国石油大学学报自然科 学版， 2 0 0 7 ， 3 1 6 1 4 1 - 1 4 6 W u Ch a o ，Ch e n M i a n ，J i n Ya n．Th e p r e d i c t i o n mo d e l o f b o r e h o l e s t a b i l i t y b a s e d o n s e i s mi c a t t r i b u t e t e c h n o l o g y [ J ] ． J o u r n a l o f C h i n a Un i v e r s i t y o f P e t r o l e u mE d i t i o n o f N a t u r a l S c i e n c e ，2 0 0 7 ， 3 1 6 1 41 - 1 4 6 ． [ 1 2 】 Wu C h a o ， C h e n Mi a n ， J i n Y a n ． A p r e d i c t i o n me t h o d o f b o r e h o l e s t a b i l i t y b a s e d o n s e i s mi c a t t r i b u t e t e c h n o l o g y [ J ] ．J o u r n a l o f P e t r o l e u m S c i e n c e a n d E n g i n e e r i n g ， 2 0 0 9 ， 6 5 3 ／ 4 2 0 8 2 1 6 ． [ 1 3 ] 李 庆忠 ．论地 震约 束反演 的策 略[ J ] ．石 油地球 物理勘 探，1 9 9 8 ， 3 3 4 4 2 3 4 3 8 ． L i Qi n g z h o n g ． O n t h e s t r a t e g y o f s e i s mi c r e s t r i c t e d i n v e r s i o n [ J ] ． O i l Ge o p h y s i c a l P r o s p e c t i n g ， 1 9 9 8 ， 3 3 4 4 2 3 4 3 8 ． [ 1 4 ] T r e i t e l S ， L i n e s L ． P a s t ， p r e s e n t a n d f u t u r e o f g e o p h y s i c a l i n v e r s i o n An e w mi l l e n n i u m a n a l y s i s [ J ] ． Ge o p h y s i c s ， 2 0 0 1 ， 6 6 1 2 1 - 2 4 ． [ 1 5 】 樊洪海． 地层孔隙压力预测检测新方法研究[ D ] ． 北京 石油大学， 2 0 0 1 ． Fa n He n g h a i ． S tud y o n c a l c u l a t i o n a n d p r e d i c t i o n me t h o d s o f f o rm a t i o n p o r e p r e s s u r e [ D] ． Be ij i n g Un i v e r s i t y o f P e t r o l e u m， 2 0 0 1 ． 【 l 6 】 金衍．并壁稳定预测理论及应用研究【 D】 ．北京石油大学， 2 0 0 1 ． J i n “ f a n．S t u d y o n pr e d i c t i n g me t h o d o f b o r e h o l e s t a b i l i t y a n d i t s a p p l i c a t i o n [ D] ． B e i j i n g Un i v e r s i t y o f P e t r o l e u m， 2 0 0 1 ． [ 1 7 】 陈勉， 金衍， 张广清． 石油工程岩石力学[ M】 ． 北京 科学出版社， 2 0 0 8 ． Ch e n M i a n ， J i n Yan ， Z ha n g Gu a n g q i n g． Ro c k m e c ha n i c s i n p e t r o l e u m e n g i n e e r i n g [ M] ． B e ij i n g S c i e n c e P r e s s ， 2 0 0 8 ． [ 1 8 】 楼一珊，金业权 ．岩石力学与石油工程[ M】 ．北京石油工业 出版 社。 2 0 0 6 ． Lo u Yi s h a n ，J i n‘Ye q ua n ．Ro c k me c h a n i c s a n d p e t r ol e u m e n g i n e e r i n g [ M] ． B e ij i n g P e t r o l e u m I n d u s t r y P r e s s ， 2 0 0 6 ． f 1 9 1 El a na y a r S ，S h i n Y C．Ra d i a l b a s i s f u n c t i o n ne u r a l n e t wo r k for a p p r o x i ma t i o n a