测井曲线的自动分层模型.doc
《测井曲线的自动分层模型.doc》由会员分享,可在线阅读,更多相关《测井曲线的自动分层模型.doc(20页珍藏版)》请在咨信网上搜索。
1、摘 要在地球物理勘探中,为了了解地下地质情况,以便于对具有不同特点的地层确定研究目标,以及确定将要重点研究的地层,统一不同井号的研究范围,其中测井曲线分层是首先要完成的基础工作。本文以1号井为标准井,建立数学模型实现了测井曲线的自动分层。在建立模型过程中,对1号井的数据进行了分类:有效值、无效值、过渡值。我们采用零替换的方法处理了题中出现的无效数据,对于其他非正常数据,由于其表现出的无规律性,因此我们采用了中值滤波的处理方法减小了噪声干扰,从而提高了数据质量。鉴于测井曲线中评价指标过多的情况,首先根据数据的特点进行了初步的筛选,剔除了信息含量少的指标。对于留下的36个指标,又根据信息论的思想,
2、计算出每一项指标的信息量,进一步剔除信息量较低的指标,最终得到22个测井曲线评价指标。该模型对这22个指标进行了主成分分析,得到五个主成分,其累计贡献率达到了80%以上,起到了降维的作用。再根据主成分的方差贡献率确定了每一个主成分的权重,然后将所有主成分加权求和得到一个新的综合指标,从而根据这一综合指标将所有的测井曲线综合为一条测井曲线,利于模型的后续处理。对于每号井的综合测井曲线,该模型采用matlab软件编程进行了趋势分析,对测井曲线进行了粗分层,确定了分界点的可能位置,然后进行了层界面归并和加权值法命名,达到了测井曲线自动分层的目的。依据建立的模型对1号井进行了自动分层,根据分层结果论证
3、了该模型的准确性程度,得出该模型有较高的准确性。然后对2至7号井进行了自动分层,通过了人工分层结果进行对比,分析了在测井曲线分层中出现的自动分层模型的准确度问题和人工分层的主观性问题。最后,利用文中建立的模型对8至13号井进行了自动分层,给出了分层结果,并进行了对结果的分析。关键字:中值滤波;主成分分析;趋势分析一、问题重述在地球物理勘探中需要利用测井资料了解地下地质情况,其中测井曲线分层是首先要完成的基础工作。测井曲线分层的目的是为了在今后的研究中,便于对具有不同特点的地层确定研究目标,以及确定将要重点研究的地层,统一不同井号的研究范围。通常,在一个区域内,通过前期地质研究工作,结合各种测井
4、数据,首先对最早开发的参考井进行详细研究。每一种测井数据,都反映了地质结构的特点和地层的变化,地质人员通过经验,综合各种测井数据反映的地层特点,将井从一定深度开始,对井进行井层划分和命名,如1号井从距井口深368米处开始,依次往下,定名为长31、长32、长33、长41、长42、长61、长62、长63、长71、长72、长73、长81、长82、长91、长92等地层。接着在分析随后开发的2号井时,也根据和1号井分层的特点和规律,依次定名为长31、长32、长33、长41、长42、长61、长62、长63、长71、长72、长73、长81、长82、长91、长92等地层。井的位置不同可能会导致这口井的每一个层
5、位的深度范围也不同,甚至有可能会出现缺失中间某层的现象。如第6号井缺长31、长32层。通常这些工作都是通过人工来进行的,这就是所谓人工分层方法。该方法不仅费时费力,而且分层取值过程中受测井分析人员的经验知识和熟练程度影响较大,主观性较强,也会因为不同的解释人员的个人标准有误差,而造成不同的人员有不同的分层结果。自动分层的基本思想、实现手段是一个不断发展变化的过程。由人工分层到自动分层,除了计算机工具的引入,各种数据处理技术也被应用于自动分层。随着一个区域开发井的数量增加,我们希望利用已有分层井点数据与变化特点作为控制点,结合每口井丰富的测井曲线数据,如密度 (DEN)、声波 (AC)、中子 (
6、CNL)、自然伽玛 (GR)、自然电位 (SP) 和电阻率 (RT) 等的变化特点,建立合理的数学模型,实现井位分层人工智能处理,也就是实现自动分层。相对于人工分层,自动分层可以避免人为分层的随意性,并可在很大程度上提高工作效率。进行具体的井位分层人工智能处理,这将极大地提高工作效率。另一方面,希望通过自动分层处理,与人工分层的结果进行比较分析,进一步提高分层精度。下面请完成以下工作:1. 以1号井为标准井,根据此井的各种测井曲线数据,建立数学模型,对第2号至7号井进行自动分层,并且通过分析,与人工分层结果进行比较分析。考虑是否需要利用你所建立的数学模型,对1号井的分层结果进行说明。2. 通过
7、前面人工分层与自动分层的比较结果,以及已给的各种测井曲线数据,确定合适的数学模型对第8号井至13号井进行自动分层,并分析你的结论。二、问题分析 在地球物理勘探中需要利用测井资料了解地下地质情况,其中测井曲线分层是首先要完成的基础工作。传统做法是采用人工分层方法,但是这种方法不仅费时费力,而且受分析人员的经验知识和熟练程度影响较大,主观性较强,存在明显的误差。我们现在希望利用已有分层井点数据与变化特点作为控制点,结合每口井丰富的测井曲线数据,建立合理的数学模型,实现井位分层人工智能处理,也就是实现自动分层。第一问中,要求以1号井为标准井,根据此井的已知数据,建立模型,对第2至7号井进行自动分层,
8、并且与人工分层结果作比较。第一步,因为在1号井的测井曲线数据中有66项指标,6000多个样本点,所以我们必须剔除掉一些无关紧要的指标,否则由于数据量太大而无法进行求解。我们通过查阅资料,首先人工初步剔除掉30项指标,然后对剩余的36项指标根据每一项中所包含信息量的大小,再次剔除掉部分指标,最终确定出22项指标作为研究对象。第二步,为了消除因仪器设备引起的测量数据扰动,我们采用中值滤波的方法处理数据,得到更有效的数据。第三步,将数据标准化处理后,采用主成分分析法将多个指标综合成少数几个综合指标,即确定出几个主成分。第四步,根据每一个主成分对分层结果的贡献率确定出每个主成分的权重,再将所有主成分综
9、合成一个综合指标,这样每一口井可以得到一条关于这个指标的测井曲线。第五步,我们对测井曲线进行趋势分析,通过计算测井曲线上各点的切线斜率,可以判断层界面的大体位置,从而完成初步分层,然后对初步分层中不必要的小薄层利用层界面归并归并为一层。第六步,我们通过计算每口井每一层的测井值,利用距离判别的思想,分别于1号井作比较,完成对井层的命名。最后,我们通过作图将第2至7号井自动分层结果与人工分层结果作比较,对比较结果作出分析,并且对1号井的分层结果进行说明。 第二问中,根据第一问中已经建立的模型,结合已有的第8至13号井的测井数据,对第8至13号井进行自动分层,并且对分层的结果作出分析。三、模型假设1
10、测井数据真实可靠,不考虑人为因素引起的误差;2不考虑地层变化引起的分层波动;3以1号井为标准,它包含所有的地质分层情况。四、符号说明符号含义第项指标的信息量第项指标的标准差第项指标的均值标准化后第项指标第个观测值随机向量的均值随机向量的协方差矩阵主成分分析中相关系数矩阵主成分的特征值主成分的特征向量五、模型建立的准备5.1 数据的预处理我们观察分析每口井的测井数据,发现每项指标测得的数据值可以分为以下几类,第一类,无效值,如-9999;第二类,有效值,如30.885;第三类,过渡值,此类值往往出现在从无效值到有效值的过渡段,如-6196.22。我们考虑到无效值会影响分层计算的结果,而且无效值都
11、是相同的值,为了便于计算,将每口井的测井数据中无效值全部替换为0。而对于过渡值,因为它没有规律可循,无法一一替换,我们将在模型中采用中值滤波的方法消除它带来的噪声干扰。六、模型的建立与求解6.1 问题一模型的建立与求解6.1.1 人工初步剔除指标在1号井的测井数据中,一共有66项指标,由于数据量太过庞大,计算机无法在短时间内算出结果,所以必须简化指标。我们通过查阅资料,人为初步剔除了30项指标,剔除后的36项指标,如下表所示,表1 人工初步剔除后的36项指标DENRILDRILMRLL8SP1R4.0SPDEViAZImGRACRMLRMNCALCNLR4.0%cnl%GR%AC%RILD%R
12、ILM%RLL8%SP1%RMN%SP%den%CAL%RML%PORTSWSHSXORTCIDCIMCIL86.1.2 根据信息量的大小再次剔除指标在信息论中,如果某项指标的变化幅度越大,我们就认为这项指标所携带的信息量越大,对结果的影响越大,也就是说这项指标越重要。而在1号井的测井数据中,测井开始和测井结束时每项指标的数据都是无效的,这些数据对分层不能提供有效地信息,所以我们取从长31的底深(294)到长92的底深(857.8)这一段测井数据作为每项指标的有效观测值,然后我们定义某项指标所有有效观测值的标准差与均值的比值大小为该指标的信息量,即 . (1)其中,表示第项指标的信息量,表示第
13、项指标的标准差,表示第项指标的均值。我们在Excel中利用Average函数和STDEVP函数分别求出1号井人工初步剔除后每项指标的均值和标准差,然后代入(1)式,求得每项指标的信息量,作出如下图:图1 36项指标的信息量图中横坐标表示每一项指标,纵坐标表示每一项指标对应的信息量,从中可以看出有些指标的信息量非常小,完全可以剔除掉,经过再次剔除后剩余22项指标,它们的信息量如下图:图2 22项指标的信息量从图2中可以看出,这22项指标的信息量基本都大于0.5,我们认为这些指标携带着影响分层的主要信息,是比较重要的指标。最终选取的22项指标如下表所示:表2 根据信息量大小再次剔除后的指标RILD
14、RILMRLL8R4.0DEViAZImRMLRMNCNLR4.0%cnl%RTRILD%RILM%RLL8%CIDCIMCIL8SHRMN%RML%PORT6.1.3 中值滤波 为了消除因仪器设备产生的测量数据扰动,我们采用中值滤波对数据进行处理。假设有数据样本序列,当滤波窗口为,其中值滤波的步骤为:(1)取以第个数据为中心的个数据行进行排序(顺序或逆序均可);(2)取排序后的中间值,即第个数据作为第点的滤波值;(3)自上而下迭代计算离散曲线上的各点。中值滤波流程如图3所示图3 中值滤波流程图输入;调用快速排序,得输出开始结束否是 滤波窗口中的的的视目的不同而定。时,滤波器失效;较小(如)时
15、可保持曲线的幅度值;较大时可只保持曲线的变化趋势。在这里,我们没有找到很好的办法确定的大小,只能经过对不同值的尝试,最终确定当时,比较合理。 6.1.4 标准化处理数据为了消除各项指标之间不同量纲的影响,对每项指标的观测值进行标准化处理,即: . (2)式中,表示第项指标的第个观测值,表示第项指标所有观测值的均值,且,表示第项指标所有观测值的标准差,且,表示标准化处理后第项指标的第个观测值。根据式(2),我们对所有指标的观测值进行标准化处理。6.1.5 主成分分析主成分分析法是利用降维的思想,在力保数据信息损失最少的原则下,把多个指标转化为少数几个综合指标的一种对多变量数据进行最佳综合简化的多
16、元统计方法。设它的数据有个指标,,分别用表示,这个指标构成的维的随机向量为。设随机向量的均值为,协方差矩阵为。对进行线性变换,形成新的综合变量,满足下式:其中在上面各式中,表示第个主成分,表示第个主成分中,第个指标的载荷。主成分分析的基本步骤如下:第一步:求协方差矩阵。对于1号井,有22项指标需要处理,则先由原始数据得到协方差矩阵。 第二步:指标的标准化处理。利用5.1.4中标准化处理后的结果,可得到新的数据矩阵。第三步:计算相关系数矩阵。其中,为原来变量与的相关系数,其计算公式为: 因为是实对称矩阵(即),所以只需计算其上三角元素或下三角元素即可。第四步:计算特征值与特征向量。解特征方程求特
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 测井 曲线 自动 分层 模型
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【快乐****生活】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【快乐****生活】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。