第七章 地质数据序列分析.pptVIP

第七章地质数据序列分析*第1页，共19页，星期日，2025年，2月5日引言在地质研究中,有很多地质变量观测值构成的数据序列。如沿地层剖面或钻井岩心剖面上地层的厚度,电性曲线和模拟地震记录离散抽样,数字地震记录等等,这些有序的地质数据描述了地质特征的变化。地质数据序列是描述地质特征的地质变量的观测值按观测顺序排列的数据序列：地质数据序列分析:研究不同的地质数据序列之间的关系及自身性质等的统计分析方法。在此介绍简单的相关分析和滑动平均。X=(x1,x2,…,xn)*第2页，共19页，星期日，2025年，2月5日相关分析本质上是一种线性滤波，是数字地震资料的一种基本处理方法。它不仅在压制随机干扰、提高信噪比中作为一种重要的滤波方法，同时相关分析原理在识别和消除多次波、计算速度谱、确定同相轴的基本参数等方面都有广泛的应用。地震勘探中一种独特的方法—连续振动法也是根据相关分析原理提出的。在此，我们把地质数据序列视为一个波形的离散抽样，并从波形的相似性出发，简单介绍相关分析在地质数据序列研究中的应用。§1相关分析*第3页，共19页，星期日，2025年，2月5日对于波形之间的相似性，有时直观上即可看出两个波形的相似性。如图7-1中的两组波形，容易看出x2(t)和y2(t)是很相似的，而波形x1(t)与y1(t)是不相似的。图7-1波形相似性一、相关函数有时却不一定能直观地做出这样的判断。即使能作出判断，也并不知道它们相似到什么程度。*第4页，共19页，星期日，2025年，2月5日若波形的总能量是一个常数(该条件在地质数据序列中常近似满足，如一张记录上的两道地震波形曲线，其总能量就是近似相等的)，那么式(7-2)中前两项就是常数。而两个波形的相似性就决定于:(7-2)σ能反映两个波形的相似程度。它越大，相似程度越低；反之，相似程度越高。展开上式得：(7-1)设x1,x2,…,xn和y1,y2,…,yn是x(t)和y(t)等间隔的离散采样值，定义它们的均方差为：*第5页，共19页，星期日，2025年，2月5日图中两个波形初看不相似,错动后可发现x3(t)的AB段与y3(t)的CD段相似。因此讨论两个波形的相似性时，不仅要从整体上分析其相似，还要分析其中一个移动到什么位置它们最相似。(7-3)ABCD图7-2波形段的相似的大小了(注意φ可正可负)。φ小则σ大，两个波形相似程度低；反之两个波形相似程度高。另外，在分析波形相似性时，还会遇图7-2的情况：*第6页，共19页，星期日，2025年，2月5日(7-4)式中τ称为数据序列的滞后项数。对一系列的τ值,计算出对应的φ(τ)值,当两个波形相似的段(AB和CD)恰好重合时,φ(τ)最大。因此φ(τ)可以做为衡量两个波形相似程度的指标,用以分析两个波形中的相似段。另外,φ(τ)不仅与两个波形本身的特点有关,还与两个波形之间的相对移动量有关。称φ(τ)为两个波形的相关函数。为解决这个问题，需要计算两个波形错动到不同位置时的各个φ值，即:*第7页，共19页，星期日，2025年，2月5日(7-5)当x(t)与y(t)相同时，φ(τ)记为：(7-6)若x(t)、y(t)是两个不同的波形，φ(τ)记为：称为x(t)的自相关函数。称为x(t)与y(t)的互相关函数。*第8页，共19页，星期日，2025年，2月5日图7-3自相关图示意图τ10周期1.自相关函数的性质2.?自相关图①当τ=0时，有正的最大函数值；②当τ→∞时，φxx(τ)=0；③自相关函数是偶函数，即φxx(τ)=φxx(-τ)。自相关图是以φxx(τ)为纵坐标,以滞后项数τ为横坐标绘制的曲线图,如图7-3所示。该图描述了数据序列的周期性。二、自相关函数的性质及自相关图*第9页，共19页，星期日，2025年，2月5日1.φxy(0)不一定有最大值，最大值在两个波形的最大相似段上;2.一般来说，φxy(τ)不是偶函数;3.τ→∞时，φxy(τ)→0。二、互相关函数的性质*第10页，共19页，星期日，2025年，2月5日1.3项滤波方程§2滑动平均数据序列中的每个观测值均由趋势变化、周期性变化和随机干扰三部分组成。应尽量抑制数序列中随机干扰部分。在此介绍抑制干扰的滑动平均法，它也是一种简单的数字滤波。下面是地质研究中常用的几个数字滤波方程：过