基于测井曲线数据进行合理自动分层的模型.doc

封一答卷编号（参赛学校填写）：答卷编号（竞赛组委会填写）：论文题目： （C）基于测井曲线数据进行合理自动分层的模型 组 别：本科生参赛学校：吉林大学报名序号：（可以不填）参赛队员信息(必填)： 姓 名所在学院及专业联系电话参赛队员1管炜汽车工程学院车辆工程420906班4209063515143086817参赛队员2曾雷汽车工程学院工业设计420907班 42090735赛队员3刘明强汽车工程学院车辆工程420903班4209033715843139536 封二答卷编号（竞赛组委会填写）： 评阅情况（省赛评阅专家填写）：省赛评阅1：省赛评阅2：省赛评阅3：省赛评阅4：省赛评阅5：38基于测井曲线数据进行合理自动分层的模型基于测井曲线数据进行合理自动分层的模型摘要本文解决的是根据各种测井曲线数据，合理自动划分地层的问题。避免了人工分层的主观性和不确定性，从而达到更为精确的分层结果。这样，在今后的研究中，就能对不同地层确定研究目标，以及确定将要重点研究的地层，统一不同井号的研究范围。对1号井分层时，我们首先建立了模型一和模型二。模型一：变量聚类模型，即对各个测井指标先进行聚类。变量幅值随深度变化大致一样的聚为一类，从每一类中挑选出具有代表性意义的指标，这样就可以有效的减少变量的个数。模型二：褶积求导分层模型，我们首先用褶积法求导代替一般意义上连续函数的求导。然后画出各个代表性指标对深度的导数曲线，找出导数的极值点，作为分层候选点。然后根据这些分层候选点，计算出层间方差和层内方差之比F，当F值达到最大时，就是最佳的分层结果。综合各个指标，我们得到1号井的分层结果如下：人工分层结果自动分层结果294290.175330331.925368366.55410405.935453.8446.3495.3494.8530.2530.8575589.175614.7606.175654.5672.425697691.55736.3738.675771.8766.3814817.8857.8832.175然后，我们以1号井为标准井，用1号井的分层方法分别对2-7号井进行自动分层（具体结果见模型求解部分）。模型三：残差分析模型，计算出我们得出的分层结果和人工分层结果之间的残差，找出与人工分层点的差别以及差别有多大。然后用模型二对8-13号井进行自动分层，并根据1-7号井的分层结果找出8-13号井的分层点（具体结果见模型求解部分）。本模型中，首先对变量进行了聚类分析，从而很有效的减少了变量个数，减小了处理难度。用褶积法求导代替一般意义上连续函数的求导，通过导数的极值点挑选出可能分层点，并用层间方差与层内方差之比来评价、优选这些分层点。关键词：测井曲线；自动分层；聚类分析；褶积求导法；层内与层间方差比；残差分析问题重述在地球物理勘探中需要利用测井资料了解地下地质情况，其中测井曲线分层是首先要完成的基础工作。测井曲线分层的目的是为了在今后的研究中，便于对具有不同特点的地层确定研究目标，以及确定将要重点研究的地层，统一不同井号的研究范围。通常，在一个区域内，通过前期地质研究工作，结合各种测井数据，首先对最早开发的参考井进行详细研究。每一种测井数据，都反映了地质结构的特点和地层的变化，地质人员通过经验，综合各种测井数据反映的地层特点，将井从一定深度开始，对井进行井层划分和命名，如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等地层。井的位置不同可能会导致这口井的每一个层位的深度范围也不同。通常这些工作都是通过人工来进行的，这就是所谓人工分层方法。该方法不仅费时费力，而且分层取值过程中受测井分析人员的经验知识和熟练程度影响较大，主观性较强，也会因为不同的解释人员的个人标准有误差，而造成不同的人员有不同的分层结果。因此，我们需要建立一个自动分层的数学模型。首先建立一个对1号井进行自动分层的数学模型，然后以1号井为标准井，对其他井进行自动分层，并与人工分层结果进行比较分析。根据比较结果，对已建立的模型进行改进，并确定合适的数学模型，以得到更为精确、科学的分层结果，最后分析自动分层的结果。本文要解决的问题有：1.以1号井为标准井，根据此井的各种测井曲线数据，建立数学模型，对第2号至7号井进行自动分层，并与人工分层结果进行比较分析。2.通过前面人工分层与自动分层的比较结果，以及已给的各种测井曲线数据，确定合适的数学模型对第8号井至13号井进行自动分层，并分析结论。问题分析首先考虑到测井指标有66项，变量过多，不利于我们的分析。于是我们利用测井曲线数据通过matlab画出各项指标随深度的变化曲线，发现有很多指标的变化趋势大致相同。这样我们就能把变化趋势大致相同的指标归为一类，从而减少变量个数，使得分析更加方便。对变量分类可以用R型聚类分析，通过matlab软件画出66个变量聚类图。根据聚类结果，分析把指标分成多少类，然后从这些分好类的变量中找出有代表性的指标，这样变量（指标）个就大大减少了，从而有利于我们的分析。针对问题一：我们可以先求出测井曲线数据对深度的导数，导数值越大，说明测井数据随深度变化越明显，作为分层点的可能性越大。但由于所给指标的数据随深度变化是离散的不连续的，故用一般的求导方法不能求出某项指标的数据对深度的导数。为此，我们用褶积法求导代替一般意义下连续函数的求导（下文导数即指褶积法求出的导数）。首先我们利用导数的极值（可以规定一个导数截止值T）找出测井曲线的拐点，这些拐点即为分层候选点。通过查找资料，我们找到了层内方差和层间方差两个评价候选点指标，并定义层间方差和层内方差之比为F。层间方差反映了各层段平均值围绕总平均值的分布情况, 它反映了层与层之间的物性差异, 层间方差大, 说明各层之间的数据差别明显, 层内方差反映了层内测量值围绕该层平均值分布的特征,它反映了层内岩性不均匀程度和随机干扰, 层内差小, 说明层内岩性均匀, 数据分布集中。因此F的值越大，分层越合理。于是我们就可以给出几个合适的导数值，分别计算F值，取F值最大时，截止导数值所分的层作为我们的分层最终结果。对1号井分好层后，以1号井为标准井，用已建立的模型对2-7号井进行自动分层。然后建立残差分析的数学模型（具体方法见模型求解部分），与人工分层结果进行对比分析，用matlab算出残差分析的各项指标,对自动分层结果进行综合评价。最后结合指标的实际含义，对1号井的分层结果进行说明。针对问题二，我们根据问题一中人工分层与自动分层的比较结果，并结合已给的各种测井曲线数据，对已建立的模型进行改进，以确定最优的分层结果，并对8-13号井进行自动分层。最后对分层结果进行总体分析。模型假设与符号说明模型假设假设1：题目所给数据（-9999.000除外）都是正确的，合理的。假设2：多个指标之间相关性大的可以选用一个具有代表性的指标。假设3：以1号井为标准井，1号井中没有的指标，其他井中不予考虑。假设4：不同岩性在曲线上有明显差别，且同一岩层的测井值应相对稳定，波动较小。假设5：相邻两个采样点之间的岩性是均匀变化的。符号说明某项指标第i个采样点所测得的数据f(j)滤波因子N滤波窗口n采样点总数D总体变差层间离差平方和层内离差平方和层间方差层内方差F层间方差与层内方差之比第i层自动分层与人工分层的残差T“导数”截止值模型的建立与模型的求解模型一的建立和求解考虑到指标数量过多，不利于我们的分析，首先我们用matlab画出各个指标随深度变化的曲线（由于数量过多，我们挑选了一些指标图作为示例）如下:从这些曲线（指所有曲线）中，我们发现有些指标的幅值随深度变化是几乎一致的。于是我们建立一个对各个指标进行变量聚类的数学模型。考虑到-9999.000的特殊性(无效数据),我们先把数据分成三类，第1-40列和42、43、45、49-53列为第一类，第55、56列为第二类,第58-66列为第三类，另外第41、44、46、47、48、54、57列除-9999.000数据外，全为0，对1号井分层无任何意义。前三类数据中-9999.000的分布各不相同。为了处理数据方便（matlab中输入的数据须为矩阵），这三类中-9999.000所在层都不予考虑。先定义一个相关系数,它反映了变量和之间的相似性程度。则每类数据所构成的矩阵中每两个指标都可以算出一个相关系数，从式中可以看出，相关系数都小于1。对相关系数进行处理，用1减去相关系数，那么这些数据可以构成一个矩阵，把它称为相关系数矩阵。然后可以根据相关系数矩阵利用MATLAB画出聚类图，如下：其中1-48代表第一类中所对应的指标。其中1-9代表第三类中所对应的指标。结合聚类图，我们得到的第一类的变量（指标）聚类结果如下表类数指标所在列数第1类1 25第2类2 5 7 20 23 24 38 第3类3 21 39第4类4 6 12 13 22 30第5类8 9第6类10 18 35第7类11 19第8类14 26第9类15 17第10类16第11类27 28第12类29 32 37第13类31 40第14类33第15类34 36第16类42 43 50第17类45第18类49第19类51 52 53随后我们通过查找资料，找出测井曲线常用的分类指标，并结合聚类结果，选出第1、3、6、7、8、10、11、14、15、16、27、31、33、34、37、42、45、58、59、61列这22项指标作为代表性指标。模型二的建立和求解在这个模型中，我们建立一个褶积分层模型。原理和基本模型：测井曲线的变化是岩性变化的反映,测井曲线上每一个曲线元可以认为是反映岩性变化的最小单元,即每个曲线元内岩性是“绝对”均匀,不可再分的,因此，地层界面一定对应着曲线元的界面,只要能找出所有曲线元的边界,则地层的界面一定包含在其中了。考虑如何将测井曲线分解成按深度排列的曲线元,若把测井曲线看成是由深度和幅度两个坐标构成的“函数”曲线,那么曲线元素的形态特征及其界面应当与这条曲线的变化率即函数的导数有关。例如,一个单一的峰内应有一个极大值,而一个单一的行内应有一个极小值,峰和谷的界面应是曲线的拐点,对应于导数曲线上特定的数值：极值点的导数为0,拐点的导数是极大值或极小值,极大值对应上升拐点,极值小对应下降拐点。因不可能写出测井曲线的幅度随深度变化的函数表达式,故不能采用一般函数的求导数方法,而只能采用通过褶积求光滑导数的方法,把离散化的数字曲线看成是采样点序号i的函数,记为y(i),为了平滑曲线上小的噪音,将函数y(i)与某一个选择的实用注波因子f(j)在i附近的-N+N范围内进行褶积,以计算曲线的光滑导数y(i),其算式如下: 式 1褶积范围包含2N个采样点（2N称为滤波窗口）。用上式计算的导数, 应当符合导数的数学性质:曲线极值点的导数为0, 上升曲线的导数大于0, 下降曲线的导数小于0, 因此, 滤波因子f(j)应当是一个反对称函数, 即f(-j)=-f(j), 考虑到极值点附近的曲线在(-N,N)范围内是对称的, 滤波因子f(j)还应当是线性的, 据此把滤波因子选为:这时式1中的A为一个常数。有了这个原理，虽然所给数据是离散的，不连续的，但是我们仍然能用光滑导数近似代替导数。上图是一段测井曲线y(i)及其导数曲线y(i)。从图中可以看出，峰和谷的界面是测井曲线的拐点, 并对应于导数曲线的一个极值点, 而终和谷的幅度高低和厚度大小与它们对应的导数值大小有对应关系。一般情况下, 测井曲线上的微小波动对分层是没有意义的, 可以规定一个导数截止值T,分层时只利用导数曲线上导数绝对值大于T的极值点(称为有T意义的极值点), 而不考虑导数绝对值小于T的极值点, 这相当于不考虑测井曲线上某些很小的变化或噪音，这样我们就能找出大致分层候选点。问题分析中已经提到，定义F为层间方差和层内方差之比，其推导如下：设第i层有个数据，该层的第j个数据记为，则测井记录的总体变差总体变差可以分解为层间离差平方和和层内离差平方和两部分：显然，、的大小不仅与数据的离散性质有关，还与数据的个数有关，数据增多，、也增加。为了消除数据个数的影响，写成方差的形式：于是我们得到F=。有了F，我们就有了评价分层合理性的标准了。F值越大，层间差异越大，层内差异越小，分层越合理。这样我们就可以规定不同的导数截止值，从而得到不同的分层结果，然后比较它们的F值，F值最大的结果，分层最合理。模型的求解：我们通过画出挑选出的19个指标的导数随深度变化的曲线，并分别用这19个指标进行分层，发现第7列（自然电位）与人工分层的结果符合的很好，因此我们用这个指标的分层结果作为我们的最后分层结果，其导数曲线如下：我们规定了一系列的截止导数T（滤波窗口选的是200）,得到一系列的分层结果，发现当分为48层时，F的值最大，即分层结果最合理。其结果如下：截止导数：T=0.01分层数：48S1 = 8.1507e+006S2 = 8.5278e+009F=9.5578e-004分层结果：193.300000000000205.675000000000221.050000000000290.175000000000312.175000000000320.925000000000331.925000000000338.925000000000359.050000000000366.550000000000379.675000000000386.175000000000387.675000000000405.925000000000423.550000000000428.050000000000431.925000000000446.300000000000462.300000000000480.550000000000494.800000000000504.050000000000512.175000000000530.800000000000532.800000000000542.300000000000559.425000000000589.175000000000606.175000000000672.425000000000675.175000000000679.050000000000691.550000000000704.300000000000724.550000000000738.675000000000741.425000000000743.675000000000756.175000000000759.925000000000766.300000000000782.175000000000808.800000000000817.800000000000820.925000000000824.925000000000831.175000000000832.175000000000结果分析：我们采用的分层方法所分层数比人工分层的层数更多，结果更精确，更细致，且和人工分层的结果符合的很好，对比如下：人工分层结果自动分层结果294290.175330331.925368366.55410405.935453.8446.3495.3494.8530.2530.8575589.175614.7606.175654.5672.425697691.55736.3738.675771.8766.3814817.8857.8832.175这说明这个指标和人工分层的结果大致符合，因此以1号井为标准井，对2-7号井选用自然电位作为主要分层依据参考指标，并兼顾其他指标（综合不同指标的F值），分别得到如下分层结果：2号井分层结果（48层）160.250000000000174.875000000000209.500000000000211.625000000000216.500000000000243.000000000000272.250000000000277.500000000000280.125000000000290.250000000000303.125000000000324.375000000000347.125000000000352.250000000000389.500000000000406.125000000000421.250000000000438.875000000000450.500000000000465.000000000000471.250000000000502.250000000000517.250000000000523.625000000000533.750000000000557.750000000000576.750000000000605.000000000000620.000000000000620.750000000000626.000000000000636.875000000000658.000000000000668.750000000000679.375000000000685.875000000000690.875000000000714.875000000000738.875000000000756.375000000000763.000000000000773.250000000000810.000000000000830.500000000000857.125000000000890.000000000000909.875000000000930.250000000000人工分层结果自动分层结果257243290290.25326324.375369352.25414421.25458465497.5502.25545.1533.75582576.75612605644.3636.875682679.375717714.8757657638118103号井以后的数据见附录一。模型三的建立和求解在本模型中，我们建立一个残差分析的模型，用我们自动分层的结果，和人工分层的结果进行对比分析。残差分析基本原理：残差是指实际观察值与回归估计值的差，即 显然，有多少对数据，就有多少个残差。残差分析就是通过残差所提供的信息，分析出数据的可靠性、周期性或其它干扰 。而在我们这个问题中，自动分层结果相当于实际观察值，人工分层结果相当于回归估计值，我们就要通过残差分析找出自动分层和人工分层结果的差别，估计这个差别有多大。我们以自动分层结果为纵轴，以人工结果为横轴，画出他们之间的线性回归图，即残差分析图，并计算出残差。模型的求解：模型的求解：通过mbtlab软件，我们画出了1号井自动分层结果和人工分层结果的残差图，如下：从残差图中可以看出自动分层结果与人工分层结果对应的较好，其中第10、15层偏差较大。我们选用自然电位为主要指标（通过比较我们发现自然电位的纵向“分辨”结果比较好），其主要应用为：划分储层，判断岩性，解释层的自然电位异常，确定地层水电阻率。然后我们选用自然电位为主要指标画出8-13号井的导数图，并根据自动分层与人工分层的比较结果，发现在自动分层结果密集处只要选一个分层点，其他的分层点均为虚假层位。于是我们根据导数图，在自动分层结果密集处只选一个分层点，结果如下：8号井：图中的F代表matlab计算的分层候选点，但是很明显有些过于集中，应该结合图形进行人工分析得到最终分层结果。分层结果：247.3296356.8467.6515.6568.9608.3638.5726.5823.4843.9903.5939.1958.19号以后的分层结果见附录二。综合分析2-13号井的分层结果，我们可以看出有些井号有缺失某层的现象，我们认为这是不同井号的地质结构不同造成的。因此，在实际研究中，对于不同井号我们应采用不同的研究方案。通过选取不同的纵向分辨力较强的曲线作为“ 主分层曲线” ,综合不同指标的方差比F值来确定最终的分层结果。模型的评价、改进与推广首先，从自动分层结果来看，自动分层所分的层数明显比人工分层结果的多，有些分层点太过密集，这说明自动分层易出现虚假层位，我们认为这是因为岩性变化剧烈造成的测井曲线在较小的深度范围内呈现强烈锯齿化而最终出现分层过细的结果。针对假层位，我们考虑在做完最优分割自动分层后利用测井指标范围合层的方法减少假层位,利用曲线上的特殊点综合判别合理的界面深度。砂泥岩互层严重或层内非均质性强使岩性测井曲线显示复杂，一个层段包含多个峰谷变化。这就必须合层,一般地可利用泥质上限和泥质下限(上、下限值应以具体地区而定) 进行合层。泥质曲线上的一个谷,如果其最小值都大于泥质上限,就将其合并于泥岩层之中。一个峰它的峰值还小于泥质下限时,它就属于同一层,将其合并于同一个储集层之中。要对其进行判断，要求有合理的判别规则,可以利用统计的方法全段采用一个标准,或是利用专家系统数据库等。在一开始，我们本来打算建立判别分析模型，根据1号井的分层结果，用判别分析的方法，对2-7号井，进行自动分层，然后对结果稍加改进，再用判别分析的方法对8-13号井进行自动分层。但是从2-7号的结果中，我们发现用matlab实现后有些井号只分为一层，这显然不符合实际，我们认为，这是由于测井曲线数据的复杂性的结果。由于不同岩性的地层其测井曲线的特征不同，同一岩性的地层可能对应于多种测井相，因此不同段应采用不同指标，这样分层更合理。另外对于不同的发分层精度要求应该设置不同的滤波窗口N值和导数截止值T。本模型不仅可以应用于测井曲线自动分层，还可以应用于不同地区经济发展，教育程度状况等的分类。参考文献1姜启源，谢金星，叶俊，数学模型第三版，高等教育出版社，2003.8。2苏金明，阮沈勇，王永利，matlab工程数学，2005.8。3江成会，测井曲线自动分层，勘察科学技术，1993年第五期。4鲍晓欢，测井曲线的最优分割法自动分层评价，海洋石油，第25卷第1期。附录附录一：3-7号分层结果如下：3号井（55）170.250000000000178.750000000000187.125000000000195.750000000000210.500000000000223.125000000000259.500000000000274.000000000000276.875000000000277.375000000000284.625000000000306.250000000000321.750000000000333.750000000000337.125000000000414.625000000000454.250000000000463.875000000000475.250000000000484.625000000000501.125000000000507.250000000000514.875000000000528.125000000000530.750000000000531.625000000000581.375000000000588.375000000000588.625000000000592.625000000000597.750000000000599.875000000000617.375000000000646.250000000000658.375000000000669.000000000000676.875000000000686.875000000000688.125000000000690.875000000000730.250000000000740.250000000000744.250000000000752.375000000000760.875000000000770.500000000000796.375000000000806.875000000000831.125000000000852.250000000000860.875000000000865.875000000000867.125000000000876.625000000000886.375000000000人工分层结果自动分层结果213210.5251251259.5297.4306.25340.5333.75377.9337.125422414.625455454.25502501.125540.5531.625578581.375621617.375663658.375700669740740.25788.9796.3754号井（41层）194.625000000000211.375000000000228.125000000000234.625000000000279.000000000000280.625000000000281.625000000000288.125000000000316.375000000000366.875000000000374.875000000000398.750000000000459.625000000000469.750000000000486.500000000000497.625000000000520.250000000000533.000000000000540.000000000000544.875000000000545.375000000000566.500000000000578.625000000000594.750000000000601.250000000000608.750000000000629.625000000000647.500000000000659.500000000000666.625000000000667.625000000000670.625000000000671.875000000000714.125000000000727.625000000000738.250000000000753.125000000000762.625000000000781.250000000000802.500000000000816.875000000000人工分层结果自动分层结果276279311316.375355.7366.875398398.75438459.625483486.5518520.25560.3566.5599.7594.75643.5647.5680671.875720714.125760762.625803802.5846816.8755号井（68层）114.750000000000126.250000000000143.875000000000150.625000000000169.875000000000176.125000000000177.250000000000190.875000000000192.000000000000194.125000000000207.250000000000269.500000000000286.375000000000307.625000000000315.500000000000327.000000000000396.000000000000438.125000000000448.000000000000467.750000000000475.375000000000516.000000000000528.625000000000542.750000000000545.875000000000548.000000000000562.375000000000579.500000000000620.750000000000640.375000000000650.125000000000662.000000000000675.875000000000676.000000000000677.875000000000700.125000000000709.750000000000727.500000000000740.625000000000747.125000000000755.375000000000757.125000000000759.125000000000825.625000000000838.125000000000847.625000000000859.250000000000867.125000000000867.750000000000871.625000000000880.625000000000887.000000000000902.500000000000910.375000000000917.125000000000918.500000000000920.375000000000922.375000000000923.750000000000946.500000000000948.125000000000953.250000000000964.875000000000973.500000000000975.875000000000976.875000000000989.875000000000998.250000000000人工分层结果自动分层结果454.4448479.8475.375492516535.4528.625618620.75661.7662698700.125741740.625780759.125821.9825.625860.4859.25900.4902.5934923.75972973.51010.7998.256号井（50层）211.025000000000212.400000000000213.275000000000239.775000000000247.775000000000257.275000000000272.650000000000273.025000000000276.400000000000283.775000000000289.525000000000298.900000000000304.525000000000332.650000000000334.650000000000335.400000000000356.025000000000363.150000000000395.900000000000435.275000000000437.525000000000440.025000000000453.900000000000461.275000000000474.900000000000488.650000000000489.400000000000491.525000000000501.650000000000511.900000000000528.275000000000545.900000000000560.650000000000586.025000000000591.900000000000615.150000000000626.775000000000638.775000000000647.650000000000654.900000000000664.775000000000670.275000000000681.400000000000687.525000000000689.775000000000706.900000000000717.650000000000741.775000000000762.400000000000790.400000000000人工分层结果自动分层结果261.7257.275310304.525350356.025396395.9433.6435.275478474.9514.5511.9558560.65602591.9641.2638.775682.1681.4717.6717.65763762.47号井（42层）331.850000000000332.475000000000333.850000000000391.100000000000414.350000000000459.975000000000465.225000000000474.475000000000494.975000000000495.975000000000498.475000000000510.975000000000522.100000000000532.600000000000534.350000000000543.100000000000551.475000000000557.975000000000563.350000000000574.225000000000584.850000000000594.100000000000602.100000000000613.100000000000623.975000000000642.850000000000649.975000000000662.225000000000672.600000000000685.100000000000718.225000000000728.475000000000738.100000000000756.350000000000767.350000000000784.350000000000793.850000000000833.725000000000843.600000000000862.100000000000888.725000000000899.725000000000人工分层结果自动分层结果334333.85368391.1410414.35454459.975495495.975528532.6576574.225615613.1652649.975690685.1733728.475770767.35810793.85860.6862.1901899.725附录二：9-13号井的分层结果如下：9号井分层结果：256.9311.2350.7435.7451.9471.7543.8572.2585.7610.7628.3656.6727.8748.0781.7821.710号分层结果：220.9278.8314.3335.1394.9429.8470,4502.9541.8563.3589.4603.3640.3647.411号井分层结果：141.8183.5220.4284313.3334.1364.1461.1495.3515.1600.8631.612号井分层结果：255.8310.8347408.6459.3471.9546.3572.5601625.1691.6720.8757.577813号井分层结果：210,4255.9281.3349.6397.4508.1553.1594.6624.4754.4769.9806.6附录三：部分计算机程序(1)部分变量的随深度变化曲线程序clear clc%数据的前处理load data.txtx=data;x(:,1)=;x1=x;%对一些-9999的无效数据进行行删除处理 同时由基本数据x1(:,54:66)=;x1(:,46:48)=;x1(:,44)=;x1(:,41)=;%和每一个变量的随深度的变化图分析出明显应x1(1:226,:)=;%该聚为一类的列也删除掉x1(5749:6151,:)=;m,n=size(x1)for i=1:44 figure Y=x1(:,i); X=data(227:5974,1); plot(X,Y,k)end(2)变量聚类分析clear clcload data1.txtx=data1;x(:,1)=;m,n=size(x)x1=x;%对一些-9999的无效数据进行行删除处理 同时由基本数据x1(:,54:66)=;x1(:,46:48)=;x1(:,44)=;x1(:,41)=;%和每一个变量的随深度的变化图分析出明显应x1(1:226,:)=;%该聚为一类的列也删除掉x1(5749:6151,:)=;%数据标准化的处理结果x2=zscore(x1);%将得到的数据进行标准化处理x3=x2;Y=pdist(x3,correlation);J=squareform(Y);figureZ=linkage(Y);dendrogram(Z,66)C=cophenet(Z,Y)I=inconsistent(Z)T= cluster(Z,0.75)k=max(T)for i=1:k r(i)=find(T=i)%r最后得到的为一个向量返回每一类的变量end（3）画出变量的褶积求导导函数图像分析以1号井为例clear clcload data1.txtx=data1;x(:,1)=;m,n=size(x)x1=x;%对一些-9999的无效数据进行行删除处理 同时由基本数据x1(:,54:66)=;x1(:,46:48)=;x1(:,44)=;x1(:,41)=;%和每一个变量的随深度的变化图分析出明显应x1(1:226,:)=;%该聚为一类的列也删除掉x1(5749:6151,:)=;%画出褶积求导的曲线N=input(请输入滤波窗口的大小（必须为一个正整数）)A1=-1:-1:-N;B1=N:-1:1;C=-N:N;F=A1;0;B1;A=dot(F,C)Y1=zeros(1,5748-N,double);%离散数据求导f