测井曲线分层.doc

测井曲线分层问题摘要测井曲线分层是在地球物理勘探中利用测井资料了解地下地质情况，首先要完成的基础工作。本文主要解决了以附件中1号井为标准井，并根据此井的各种测井曲线数据，建立了数学模型，对第2号至7号井进行自动分层，并且通过分析，与人工分层结果进行比较分析，进一步对1号井的分层结果进行说明；对第8号井至13号井进行自动分层，并给出结论两个问题。针对问题一，本文首先通过查资料并观察附表中1号井的数据剔除了一些变化规律不明显的指标如CAL、DEVi等，筛选出了SP、GR、AC 、CNL 、RT、 RILD、 RILM七个显著变化的指标，根据经验又剔除了各个指标中明显错误的数据；然后利用主成分分析法的思想挑选了在主成分中贡献率较大的指标AC、CNL、SP三个指标。接着画出三个指标的综合测井曲线，由于每一层的指标差异性、稳定性，本文采取了层内差异法，结合综合测井曲线，将每一个井进行了大致的粗分层。最后要将相似程度高的层进行合并，而聚类分析是根据某一分类统计量来度量多个观察量之间的相似程度，依相似程度高低决定是聚合为一层，还是划为不同层，本文利用聚类分析法将第2至7号井进行细致的分层，与人工分层进行了比较，判断其精度，结果见表4、表5并对模型进行了改进，进一步提高合理性。针对问题二，本文利用问题一中所得出的规律对第8号井至13号井进行了分层，结果见表6，并进行了分析。关键词： 测井曲线自动分层 主成分分析 层内差异法 聚类分析 一、问题重述测井区县分层是在地球物理勘探中利用测井资料了解地下地质情况最基本也是最重要的问题。目前最常用的人工分层方法不仅费时费力，而且分层取值过程中受测井分析人员的经验知识和熟练程度影响较大，主观性较强，也会因为不同的解释人员的个人标准有误差，而造成不同的人员有不同的分层结果。本文主要解决的问题有：1、以1号井为标准井，根据此井的各种测井曲线数据，建立合适的数学模型，对第2号至7号井进行了自动分层，并且通过分析，与人工分层结果进行比较分析并改进了数学模型，对1号井的分层结果进行说明。2、通过前面人工分层与自动分层的比较结果，以及已给的各种测井曲线数据，对第8号井至13号井进行自动分层，并分析本文的结论。二、问题假设1 、 假设题目附件中所给数据仪器的精度都满足要求2 、 1号井所给的分层都是准确的，其他井依靠1号井为标准。3 、 根据某一条或者某几条测井曲线可以较准确的进行分层符号说明 第i个深度样本的第j个指标 B 待定参数 第i个深度样本的均方根误差 第i层测井均值 非地层因素引起的允许误差 第k 层与第k + 1 层的测井值之间的“马氏”距离 各层间测井值的最小“临界距离”三、问题分析实现井位分层人工智能处理，也就是实现自动分层的方法有很多，比如神经网络、层次分析、层内聚类分析。由于问题中所给数据繁多，本文选取了主成分分析法、层次分析法的思想，客观的反映了系统状况。主成分分析法是一种将多维因子纳入同一系统进行定量化研究、理论成熟的多元统计分析方法。通过分析变量之间的相关性，使得反应信息重叠的变量被某一主成分代替，减少了变量的数目，从而降低了评价的复杂性。在以方差贡献率作为每一主成分的权重，每一个主成分的得分加权即可完成对每一层的评价。层内分析法是一种定量与定性相结合的、系统化的、层次化的系统分析方法。以随机数学为工具，通过大量的观察数据寻求统计规律。一般来说引起测井曲线值变化的因素有两类：一类是地层因素（岩性、空隙流体性质）一类是非地层因素如井壁因素、测量系统测量条件等。非地层因素引起的测井相应一般比较小，本文采剔除了这些因素。对于地质因素本文采用主成分分析的思想选取了贡献率较大AC 、CNL、RT三个指标，然后对1号井的每一层数据进行分析。由于附件中所给数据繁多直接应用主成分分析法比较繁琐，本文根据数据的规律剔除了一些明显错误的点，然后利用层内差异法找出了1号井层与层之间的区别。又利用得出的规律对2号井进行了分层，并与人工分层进行了比较，并对分层规律进行了修正， 进一步提高了分层的精确性。四、模型的建立与求解主成分分析法（1）设每一深度的测量值为一个样品，每一个样品观测有AC 、CNL、GR、SP等指标以表示第i个样品第j个指标的得分，根据附表中的数据得到原始矩阵：X=其中=根据主成分分析的方法，分析不同深度各指标的检测数据。做变换： （）得到标准化的数据矩阵=，其中。(2) 在标准化数据矩阵Y=d的基础上计算p个原始数据指标相关系数矩阵R=,其中：，(i=1n;j=1p)(3) 求相关系数矩阵R的特征值并排序；再求出R的特征向量Z(4) 确定主成分数目。在确定主成分之前，需要先给出一个控制值，令，则对应满足条件的q 的最小值即为保留的主成分的个数m。(5) 计算综合得分。首先计算得到的第i 个样本中第k个主成分的得分为， 用数据矩阵x的p个向量（即p个指标向量）， ,做线性组合（即综合指标向量）为：再以m个主成分的方差贡献率为权重求得第i个样本的综合得分。（6）分别以第一主成分，第二主成分，第三主成分中线性组合中系数的绝对值得的作为评判标准。表1 主成分贡献率 主成分的主要贡献率特征值贡献率累计贡献率第一主成分1195.00.59100.5190第二主成分719.40.24020.8316第三主成分219.80.09540.9270表2 三个主成分的7个指标系数第一主分系数0.33630.1235-0.89670.23760.1001-0.03200.0000第二主分系数0.38350.3424-0.0049-0.82930.21840.0149-0.0000第三主分系数0.21450.53760.34110.47020.57050.04640.0000本文选取了第一到第三主成分中系数的绝对值比较大的作为分层的标准，由于第一主成分所占的权重比较大，本文以第一主成分中绝对之比较大的AC作为分层的主要根据，SP、CNL两个指标为辅助分层依据。为了比较明确的看出分层结果，由他们的贡献率作为他们的权重，得到综合测井曲线为：，则M随深度变化的综合曲线为图一 综合测井曲线由综合测井曲线图联系一号井的标准分层数据可得下图： 图二 一号井的综合测井曲线图与分层图利用综合测井曲线对2到7号井进行分层1、首先运用层内分析法将井进行粗分层层内差异法层内差异分层法的依据是：同一层内的测井值相对稳定，其值的变化不超过某一允许误差(即由非地层因素引起的测井值误差) ，并认为每一层内采样值的均值代表该层的真实测值。一相邻的采样点的值与该均值的差异在允许误差范围内， 则认为该采样点属于这一层，否则，便属于下一层。其一般过程为：先选择具有较强纵向分辨率的测井曲线(称为主动曲线) 进行细分层(本区选择自然电位曲线和微电极曲线) ，然后再兼顾其它测井曲线的特点，将分层结用到其它曲线上。下面对本方法加以详细叙述。首先，为减小因测井仪器标准化程度不高、泥浆、矿化度变化等非地层因素对测井值的影响， 必须对测井曲线作归一化处理，以突出测井值的相对变化，归一化的公式为：X=式中，x 为实际测井值、为本条测井曲线所有采样点中的最大值和最小值，而X 则为归一化后的测井值(0 X 100) 。假设已确定出相邻的几个采样值( j = 1 , 2 n ) 均于第i 层, 该层的允许误差为，均值为,方差为，其中,现在的问题是要判断随后的第n + 1 个采样值 ，是否同属于第i 层。判定规则如下：若则认为属于第i 层，并计算这n + 1 个采样点的均值、方差及允许误差，接着进行下一采样点的处理。若 ，则认为点不属于第i 层，该层划分完毕，把前n 个采样点的测井均值及第n 个采样点对应的深度输出。并从点开始进行下一层的划分。如此进行下去，直至处理完整个井段为止。不难看出，上述分层的关键是确定允许误差函数，这里尝试用误差理论与概率统计相结合的办法建立这一函数。从概率统计角度讲，可以认为测井值X 是一个具有有限方差的随机变量。同一层内各采样值的差异反映了非地层因素引起的随机误差， 并满足切比雪夫不等式:式中均方根误差；该层的某个采样值； 第i 层测井均值； 非地层因素引起的允许误差 。由误差理论可求得均方根误差为:且 式中， B 为待定参数， 其大小的选取应视实际所需分层的详细程度而定。B 值越小， 允许误差越小，就越粗。所分层就越细； 反之所分层就会越粗。在这里本文人为的定义B=3，根据一号标准井可以得到每一层的允许误差。根据层内差异法，1号井的分层范围为：表3 一号井分层范围分层长31长32长33长41长42长61长62长63最小值217.29166.10220.84206.08225.6232.85228.53159.2最大值306.85352.00282.7298.94273.64261.19257.81347.54分层长71长72长73长81长82长91长92最小值180.40202.07106.7159.36212.63184.0743.42最大值302.58271.49414.7326.30259.01300.24487.58由层内分析法分层后分层过细，共统计了7个井15小层的分层数据522个,发现误差小于0. 5m 的分层数据为314 个,符合率为60%。，于是对模型进行改进，由聚类分析法将分的过细的层进行合并，方法如下。2、用层内聚类分析法将井进行细分层：在上面的分层处理中,对于整个处理井段均采用了相同的B 值, 而且在精细分层中B 值取得较小,导致实际操作中可能会将某些层段划分过细,把原本同一地层分成了若干小层。我们在层内差异法分层的基础上再结合聚类分析法, 以避免上述问题的出现。聚类分析是根据某一分类统计量来度量多个样品(观察量) 间的相似程度, 依相似程度的高低决定是聚合为一类, 还是划为不同类。这里我们选用马哈拉诺比斯距离(简称“马氏”距离) 作为分类统计量来表示已划分出的各小层间的相似程度, 把相似程度高、且相邻的小层合并为一层。具体作法如下:若在经层内差异法分层的数条曲线中有L 条测井曲线纵向分辨率较强, 按前述方法共分地层总数为M ,这时计算第k 层与第k + 1 层的测井值之间的“马氏”距离：根据计算结果,若,则认为第k 层与第k + 1 层属于同一小层，故进行并层处理;若 则不作并层处理。其中为各层间测井值的最小“临界距离”,可由经验选择。层内差异法的流程图如下图：测井曲线归一化各变量赋初值（0）N=N+1调入一个采样值计算前N个采样值的的均值的均值和方差计算允许误差当前采样样值均值不大于与允许误差N=1输出最后一层测井均值及起止深度结束是否并层最后一个样本点计算最后一层各测井均值最后一个采样点输出前一层测井均值及起止深度NYNNYNY图三 层内分析法的流程图B与的数值选取分析1、B值和值的选取:非均值地层中的人机联作方式,实际上是一种人工干预,参数是针对某一口井或某一地区而言的,具体在确定分层界限的有关经验规律上及分层处理技巧上。也就是说B值和 值是一个统计数据,适用于统计区的地层特征,它虽然受岩性影响,但其大小取值对岩性而言并无规律性,即在非统计区不适用。2、选取合适的d 、B 值后,一些井尚需做“并层”处理。统计发现,多需删除曲线呈负斜率变化的分层,曲线由大变小处的分层结果。这时若上、下二分层所夹地层厚度 (2. 5/ 3. 0) m ,则不再删除。3、并层处理后若二分层间的地层厚度 (6. 0/6. 5) m ,可从其间分出另一层。4、若二分层之间的地层厚度较小,如小于1 m时,可进行并层处理。需要说明的是:上述技巧是在对本区小层厚度分布进行分析、对比之后提出的,仅在本区有效。图四为：一号井的人工分层和自动分层对比图：（彩线为人工分层，红线为自动分层）图四 一号井人工与自动分层图其柱状图为：图三 一号井人工分层与自动分层的比较利用此测井曲线进行自动分层的框图:表4 2号到7号井的自动分层结果长31长32长33长41长42长61长622号人工257 290 326 369 414 458 497.5 2号自动239.6263.6302.75327.75358.64254483号人工213 251 297.4 340.5 377.9 422 455 3号自动210.75277323.125347.125374.875412.8754484号人工276 311 355.7 398 438 483 518 4号自动280.125347.5383407.54344434995号人工450.4 492 535.4 479.8 618 661.7 698 5号自动396.125450.625504.625541.625580.375622.3756656号人工261.7 310 350 396 433.6 6号自动250.125298.4330.65375.5408.6527号人工334 368 410 454 495 528 576 7号自动230.35260.1318.1371.725434.975474.225529长63长71长72长73长81长82长91长92545.1 582 612 644.3 682 717 765 811 497.125541.875607642.8668710.75741790.125502 540.5 578 621 663 700 740788.9 501553.75586.125608.625661696747.375784.15560.3 599.7 643.5 680 720 760 803 846 518560.125612.125656.625730.25796834860.125741 780 821.9 860.4 900.4 934 972 1010.7 682741.5786829.5867.375901.125964.8751000.5478 514.5 558 602 641.2 682.1 717.6 763 439.775477.525508.275561.4601.4643.15664.125701.275615 652 690 733 770 810 860.6 901 560.225615.85642.1676.975730.35773.35853.85872.35图四 二号井的人工分层与自动分层比较图图五 3号井的人工分层与自动分层比较图图六 4号井的人工分层与自动分层比较图图七 5号井的人工分层与自动分层比较图图八 6号井的人工分层与自动分层比较图图九 7号井的人工分层与自动分层比较图对二、三、四号井进行误差分析：表5 2号到4号井的误差分析长31长32长33长41长42长61长62二号井0.06770.09100.07130.11180.13380.07210.0995三号井0.01060.10360.08650.01950.00800.002160.0154四号井0.01490.11740.07680.02390.00910.08280.0367长63长71长72长73长81长82长91长920.08800.06890.00820.00230.02050.00870.03140.02570.00200.02450.01410.01990.00300.00570.01000.00600.07550.06600.04880.03440.01420.04740.03860.0167对二到七号井的误差分析可以得出，误差不大于0.15，可见精度还是符合要求的。本文用此改进后的模型对八到十三号井进行自动分层。结果如下：图6 8号到13号井的自动分层结果长31长32长33长41长42长61长628号井96.5132.5172.625230.75269.625327.5360.1259号井161.2195.75236.5280.75326.375356.95401.7510号井72.5138.25175.75242.75280.625342.62538911号井224.125260.5309.965355.5387424.375472.7512号井181.5233282337.75381441.25491.37513号井232.225278.85307.6345.725385.475432.225482.6长63长71长72长73长81长82长91长92402.875450.5500550.125602.25630679714442.325492.825534.45575.75629.75672.575731.325782.825428.75461.625506.125546.125599.625639.225682.375700524.25565.125611.375667.25698.25751761.125839.25541.25599.5636.5705749.5775.75815857.125554.6603.1640.35687.6717.415757.85782.725834.475对结果进行分析：1、8到13号井没有缺层现象。2、对于同一层其深度范围并不相同，说明井的位置不同可能会导致每一层的层位范围不同，缺失现象也不相同。五、模型的优缺点及改进模型优点：1、 运用主成分分析法将数据进行了简化，方便了计算。2、 运用层内分析法与聚类分析法结合起来进行粗细分层，较好的对各井分层。3、 采用内部一致性原则，将关键变量构造成综合曲线，不仅能简化计算，而且能够形成一条各层差异显著的曲线，对于后面的分层过程起到促进作用。模型缺陷：1、 数据量较大，为了简化模型，可能没有有效利用数据。2、 依据标准井计算出来的允许误差以及“马氏”最小距离，由于人为经验因素，可能会对结果有一定的影响。模型的改进：由于本文对B以及马氏距离的确定是根据以往的经验人为确定的，所以误差是肯定有的。为了减小这个误差，可以对模型再进行改进。对允许误差的选定可以参考一号标准井的数据，先选定几个B，将结果与一号井进行对比，求出每个B的误差，以此方法，可以选出一号井每一层较为合适的B，以此作为其他井的标准，这样就可以将误差大大的减小。 参考文献1 程玉群、李秀荣、刘铁庄，测井曲线自动分层技术及在杏北地区小层划分中的应用，地质勘探，第7 卷第2 期2 /question/219318664.html?an=0&si=13雍世和、陈钢花、白康生，测井曲线自动分层，测井技术Well Logging Technology 编辑部邮箱1987年06期4周开利、邓春晖，Matlab基础及应用教程,北京大学出版社5韩中庚，数学建模方法及应用，高等教育出版社6 卓金武，Matlab在数学建模中的应用，北京航空航天大学出版社附件 对选出来的指标进行主成分分析,其程序如下：Clearclcload data2.txtX=data2;Z=zeros(size(X);Z(:,1)=(X(:,1)-mean(X(:,1)/std(X(:,1);Z(:,2)=(X(:,2)-mean(X(:,