资源描述：
《综合测井数据采集及处理系统的研制_数据处理技术new》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、第30卷第2期成都理工大学学报(自然科学版)Vol.30No.22003年4月JOURNALOFCHENGDUUNIVERSITYOFTECHNOLOGY(Science&TechnologyEdition)Apr.2003[文章编号]167129727(2003)0220205204综合测井数据采集及处理系统的研制数据处理技术1212王艳梅朱双华童纯菡谭显方(1.成都理工大学应用核技术与自动化工程学院,成都610059;2.中国人民解放军63653部队)[摘要]该数据采集及处理系统具有自动采集测井数据、自动识别钻孔剖面地质构造与岩石岩性、解释结果自动显示等功能。采用数字滤波、曲线平滑等处
2、理技术,消除了测井资料中的随机干扰。采用数理统计等技术,归纳出极值方差分层的数学模型,成功地对综合测井数据进行了自动分层处理。[关键词]测井;模型;数据处理[分类号]P631.81[文献标识码]A综合测井的主要目的是确定钻孔剖面的地质孔测井资料进行了数据采集和处理,取得了满意构造和岩性。以往的测井资料采用人工解释,效率的结果,不仅提高了解释的精度和速度,而且提高低,精度不高,丢掉了许多有用的信息,并且不同了解释的可靠性。系统流程图见图1。的人员或同一人员在不同的时刻解释的结果差异该数据采集软件是在Dos环境下,利用较大。为了克服上述问题,本文利用数据处理、数TurberC编制而成;数据处理
3、软件是在Windows理统计、地质解释模型、计算机等技术,研制出了环境下,利用BorlandC++编制而成。系统的主自动识别钻孔剖面地质构造、岩性,自动处理和解要技术指标为:采样率1个öcm;准确识别出地质释测井资料的一整套方法。利用这套方法对某钻构造、岩性,深度误差(0.5m;测值判读与计算误图1综合测井数据采集及处理系统流程图Fig.1Flowchartofdatacollectionandprocessingsysteminsyntheticallogging[收稿日期]2001210218[作者简介]王艳梅(1964-),女,博士生,核技术与计算机数据处理专业.(E2mail:wy
4、mtx66679@163.com)·206·成都理工大学学报(自然科学版)第30卷差≤5%;并且具有良好的用户交互界面。系长度;c为电极系半长度;a为电极系半径;a1为钻孔半径;Qs为现场实测岩体的视电阻率值1测井资料预处理预处理后的曲线与原始曲线的形状十分相1.1滤波处理似,没有发生崎变,并且去掉了随机干扰,为下一滤波处理的目的是去掉测井信息中与岩石性步解释自动分层及取值奠定了基础。质无关的随机干扰。消除随机干扰的方法有简单2自动分层及取值平滑法、加权平滑法和最小平方曲线拟合平滑法。其基本原理是,将一个长的数字序列分段处理,每2.1分层的数学模型及方法段的长度与滑动段相当,先取序列中第一
5、项到第计算机分层的数学模型是用数理统计等方法2N+1项,构成一个M阶多项式,用这个多项式建立起来的。按照测井曲线的形态,可归纳为与层计算2N+1项中中心项的值(即趋势值);然后用中心对称和不对称两大基本类型。根据常规解释第2项到2N+2项重新构造一个多项式,再用这经验;对于与层中心对称的曲线,利用最大斜率法个多项式计算2～2N+2项的中心值,这样每次和数理统计方法分层效果较好,对于与岩层中心向前滑动一个数据点,便可去掉其中的干扰成分,不对称的曲线,利用极值点和数理统计方法分层最后得到一条较光滑的曲线。通过对比三种平滑效果较好。为了研究效果明显且适用于上述两类方法对原始测井曲线的处理效果,发
6、现简单平滑曲线的分层方法,通过归纳两类曲线分层方法原法对分离出大的趋势变化能力较强,但平滑后的理提出了极值方差分层方法。数据变化大。最小平方曲线拟合平滑法能较好地极值方差分层的基本思想是,先利用“分层指逼近原始数据,但去掉随机干扰能力较差。加权平标函数”进行粗分层,大致确定出层界面位置,再滑处理后的数据与原始数据的逼近程度较高,去对相邻两层用方差分析分层法找出准确的层界面掉随机干扰能力较强,因此选用加权平滑法对原位置。具体方法如下:始曲线进行滤波处理。具体方法如下:2.1.1分层指标函数的确定若取滤波因子h(t)=(h(-N),⋯,h(0),方差分析分层的目的是找出层间方差最大、⋯,h(N
7、)),并且要求:层内方差最小的点作为分层点。对应两层介质,设N6h(n)=1有N个采样点,若层界面在n与n+1采样点间,n=-N则两层的层内方差和Q为:用h(t)对X(t)滤波,得到加权平滑公式为:nNN22Q=6(V1j-V1)+6(V2j-V2)Y(t)=h(t)X(t)=6h(t)X(t-n)(1)j=1j=n+1n=-NnN要求加权滤波因子h(t)满足:h(t)≥0;2222=6V1j-nV1+6V2j-(N-n
