测井曲线自动分层模型设计与实现[J].docx

５ ｏ．５ ｏ２ 　　０７３ 　ＥＲＬＵ ＤＩＬＮ　ＥＨＱＵＳ Ｍｒ２０　 　 钻井与完井 　 测井曲线自动分层模型设计与实现 纪荣艺　樊洪海　杨雄文　杨皆平 中国石油大学北京石油工程教育部重点实验室北京昌平　０２９ 摘　要：从实用、有效的角度，阐述了测井曲线自动分层的步骤和多种计算模型。提出加入前处理和后处理的具 体方法，以此来消除因测量误差对模型的干扰。阐述了数理统计和聚类分析在自动分层模型设计中的应用，并开发了 相应的软件，实现了测井曲线的自动分层，使理论得到有效的应用，同时也验证了模型的可行性和有效性。 关键词：测井曲线；自动分层；中值滤波；聚类分析；趋势分析 中图分类号Ｐ３８４　文献标识码：Ａ　　文章编号１００９２０００２０ 　　测井数据的计算机自动分层是在数理统计等方 法的基础上建立起来，通过实际应用得到了有效检 验和发展。由于测井数据误差的影响，多年来自动 分层相对于人工经验分层效果并不令人满意，所以 为有效消除测井数据误差的干扰，笔者按前处理、 模型计算和后处理３步进行自动分层。前处理的主 要任务是采用数据滤波器的方法消除因仪器设备产 生的测量数据扰动。模型计算是自动分层的核心， 目前常用的有：方差统计、极值方法、聚类分析和 趋势分析等几种。后处理的任务是使自动分层的结 果更加清晰化和实用化具体流程见图１ 当滤波窗口为２＋，则中值滤波的步骤为 １取以第个数据为中心的２＋１个数据并 进行排序顺序或逆序均可； ２取排序后的中间值即第狀１个数据作为 第犻点的滤波值； ３自上而下迭代计算离散曲线上的各点 中值滤波流程如图２所示。滤波窗口中狀的取 值视目的不同而定。狀＝０时，滤波器失效；狀较 小 狀１时可保持曲线的幅度值狀较大时可 只保持曲线的变化趋势。 图２中值滤波流程 处理起始和终止的狀个值时即当＜ 犻 犖－时让狀分别等于－１和犖－，这样可保 证滤波前后的样本数犖相等。 图１分层计算总流程 １自动分层前处理 １１中值滤波 假设有数据样本序列狓 ＝１２…，犖， 收稿日期２００２；改回日期２０１４ 作者简介：纪荣艺１７ ）男１９年毕业于大连理 工大学２０年获硕士学位，在读博士研究生。 联系电话：００８７３２　　第３卷第２期　　　　　　　　　　　 石　　油　　钻　　探　　技　　术　 　　　　　　　　Ｖｌ３，Ｎ．　年月 　　第３卷第２期　　　　　　　　　　　 石　　油　　钻　　探　　技　　术　 　　　　　　　　Ｖｌ３，Ｎ．　 年月 ，ａ．０７２ Ｐ 　 ＴＴＯＥＭ ＲＬＩＧ ＣＮＩＥ　 （ （ ） ， １２４） ： ．＋ ：０１８０（０７）２４４６１ 　 。 。 ） 犻 狀 （ ） ） ， ＋ ） 。 （如＝） ； 　 ， 犻狀或＞ 　 　 　 ：０６８３ ： ００６２ （ —， ，９７９４ ，０２ （１）９３２１ 第３卷第２　　　　　　　　纪荣艺等：测井曲线自动分层模型设计与实现 ·２· 非线性的中值滤波有如下特性 大犙越小相反狇狀 ＋１越小犙越大 １可以消除曲线的尖峰干扰 且因计算 狀狀１比计算犙容易故称 狀 ２可以使阶梯函数通过且阶梯位置不变 狀１为分层指标函数 ３可以通过斜坡函数且位置和斜率均不变 ２１２分层指标函数在层界面两侧的变化规律 特性３是其它滤波方法难以达到的作为前 通过分析分层指标函数狇在层界面的变化情 处理工作的一部分中值滤波必不可少 况就能够得出准确的层面位置 １２　归一化处理 若正确的层界面点在 狀１间的犅点则两 狀 狀 狀 ＝１ ＝＋ 生的影响。归一化公式为： 狀 狀１ 犖 狀 犻 狀 狀１ 犖 狀 ＝１ ＝＋ 狓 １２ 狓 ），所以狀 狀 ，因狇狀狀１＞ （ 。 比 狓 １２ 狓 狀１ 。 １＜２ 狀 １ 狓 狓 狀 狀 狀 界面处仍为极大值 狓＜２狓＞＋１情况下或 ２ ＝１ 犖 犖 ＝１ ．．　 犖 犖 狀 ２ 先计算犖－１个狇值狇１２狇２３， ＝１ 犖１ ， ２ 极值点处作为粗分层的层界面位置 狀 犖 将每个粗分层界面两侧的犿个测井数据利用 方差分析进行二次分层。具体方法是：将第一个采 犖 犖 样点的数据算作一组其余的测井数据为一组进 个测井数据划为一组，余下的为另一组，又计算出 犖－狀犼＝１ 一个犚值这样继续下去