克里格方法在大地电磁静校正中的应用

合集下载

克里金插值法的详细介绍。kriging。

kriging 插值作为地统计学中的一种插值方法由南非采矿工程师D.G.Krige于1951年首次提出，是一种求最优、线形、无偏的空间内插方法。

在充分考虑观测资料之间的相互关系后，对每一个观测资料赋予一定的权重系数，加权平均得到估计值。

这里介绍普通Kriging插值方法的基本步骤：1.该方法中衡量各点之间空间相关程度的测度是半方差，其计算公式为：h为各点之间距离，n 是由h 分开的成对样本点的数量，z 是点的属性值。

2.在不同距离的半方差值都计算出来后，绘制半方差图，横轴代表距离，纵轴代表半方差。

半方差图中有三个参数nugget（表示距离为零时的半方差）,sill（表示基本达到恒定的半方差值）,range(表示一个值域范围，在该范围内半方差随距离增加，超过该范围，半方差值趋于恒定)。

利用做出的半方差图找出与之拟合的最好的理论变异函数模型（这是关键所在），可用于拟合的模型包括高斯模型、线性模型、球状模型、指数模型、圆形模型。

----球状模型，球面模型空间相关随距离的增长逐渐衰减，当距离大于球面半径后，空间相关消失。

3.用拟合的模型计算出三个参数。

例如球状模型中nugget为c0，range为a，sill为c。

4．利用拟合的模型估算未知点的属性值，方程为：，z0为估计值，zx是已知点的值，wx为权重，s是用来估算未知点的已知点的数目。

假如用三个点来估算，则有这样权重就可以求出，然后估算未知点。

（上述内容根据《地理信息系统导论》（Kang-tsung Chang著；陈健飞等译，科学出版社，2003）第十三章内容进行总结，除球状模型公式外其余公式皆来自此书）下面是本人自己编写的利用海洋中断面上观测站点的实测温度值来估算未观测处的温度的Fortran程序，利用距离未知点最近的五个观测点来估算未知点的温度，选用模型为球状模型。

do ii=1,nxif(tgrid(ii,1)==0.)thendo i=1,dsite(ii)!首先寻找距离最近的五个已知点位置do j=1,nhif(d(mm(ii),j).ne.0.or.j==1)thenhmie(j)=d(mm(ii),j)-dgrid(i)elsehmie(j)=9999end ifhmid(j)=abs(hmie(j))end dodo j=1,nhdo k=j,nhif(hmid(j)<hmid(k))thenelsem1=hmid(j)hmid(j)=hmid(k)hmid(k)=m1end ifend doend dodo j=1,5do k=1,nhif(abs(hmie(k))==hmid(j))thenlocat(j)=kend ifend doend dodo j=1,4do k=j+1,5if(locat(j)==locat(k))thendo i3=1,nhif(abs(hmie(i3))==abs(hmie(locat(j))).and.i3.ne.locat(j))thenlocat(j)=i3exitend ifenddoendifenddoenddo！然后求各点间距离，并求半方差do j=1,5do k=1,5hij(j,k)=abs(d(mm(ii),locat(j))-d(mm(ii),locat(k)))/1000.end doend dodo j=1,5hio(j)=sqrt(hmid(j)**2+(abs(latgrid(ii)-lonlat(mm(ii),2))*llat)**2 $ +(abs(longrid(ii)-lonlat(mm(ii),1))*(1.112e5* $ cos(0.017*(latgrid(ii)+lonlat(mm(ii),2))/2)))**2)/1000.end dodo j=1,5do k=1,5if(hij(j,k).eq.0.)thenrleft(j,k)=0.elserleft(j,k)=sill*(1.5*hij(j,k)/range-0.5*hij(j,k)**3/range**3)end ifif(hio(j).eq.0.)thenrrig(1,j)=0.elserrig(1,j)=sill*(1.5*hio(j)/range-0.5*hio(j)**3/range**3)end ifend doend dorrig(1,6)=1.rleft(6,6)=0.do j=1,5rleft(6,j)=1.rleft(j,6)=1.end dotry=rleftcall brinv(rleft,nnn,lll,is,js)ty1=matmul(try,rleft)！求权重wq=matmul(rrig,rleft)!插值所有格点上t,sdo j=1,5tgrid(ii,i)=tgrid(ii,i)+wq(1,j)*t(mm(ii),locat(j)) sgrid(ii,i)=sgrid(ii,i)+wq(1,j)*s(mm(ii),locat(j))end doenddoendifenddo。

大地电磁测深法基本原理及应用

第二类磁暴与磁亚暴。这种地磁扰动的特
征是磁场强度变化剧烈，尤其是水平分量变化很大，呈现极不规则形状。
第三类地磁脉动。这是一种具有似周期振
动的特殊的短周期振动，地磁脉动是大地电磁测深最重要的场源。其周期范围一般为0.2～1000秒，振幅一般为百分之几到几十个纳特。
大地电磁场特征
大地电磁测深简介
1、20世纪50年代，法国的Cagniard和前苏联的 Tikhonov提出了大地电磁法（MT ）；
2、20世纪60年代的Berdichevski等（1969），提出了音频大地电磁法（AMT）；
3、1971年和1978年，Goldstein和Strangberg提出了可控源音频大地电磁法（CSAMT）。
4、2000年何继善院士提出广域电磁法。
优点 1、不受高阻层屏蔽、对高导层分辨能力强； 2、横向分辨能力较强； 3、资料处理与解释技术成熟； 4、勘探深度大、勘探费用低、施工方便； 5、资料处理和解释技术成熟。
缺点 1、体积效应，反演的非唯一性较强 2、纵向分辨能力随着深度的增加而迅速减弱 3、信号不稳定、不规则，容易受到工业噪声干扰
5、极化特征
不同周期的场和不同时间的场的极化方式具有明显的差异。为了在测深资料分析处理时获得稳定的阻抗张量元素，需要场源具有多
大地电磁场源
地球强大的磁场是保护人类免于遭受外太空各种致命辐射的生死屏障，然而日前，英美科学家发现，在过去的200年内，地球的磁场正在急剧地衰弱。科学家们
理论基础：Maxwell方程组麦克斯韦方程组描述了电
磁场最根本的规律，在时
间域中的表示式为：
LE dls B tdS lH dlIs D tdS

CSAMT静校正的克里格方法研究

Ｊｉａｎｇｘｉ，Ｃｈｉｎａ）
Ａｂｓｔｒａｃｔ：Ｔｈｅｓｔａｔｉｃｓｈｉｔｆａｆｆｅｃｔｓｔｈｅｉｎｔｅｒｐｒｅｔａｔｉｏｎｏｆｔｈｅｄｅｅｐｅｌｅｃｔｒｉｃａｌｃｏｎｄｕｃｔｏｒ，ＳＯｗｅｍｕｓｔａｐｐｒｏｐｒｉａｔｅｃｏｒｒｅｃｔｉｔ．Ｃｏｎｓｉｄｅｒ－ｉｎｇｏｎｔｈｅｌａｔｅｒａｌｃｏｒｒｅｌａｔｉｏｎｏｆｓｐａｔｉａｌｓｔｒｕｃｔｕｒｅｏｆｔｈｅｇｅｏｌｏｇｉｃｓｅｃｔｉｏｎ，ｗｅｉｎｔｒｏｄｕｃｅｔｈｅｃｏｋｒｉｇｉｎｇｏｆｇｅｏｓｔａｔｉｓｔｉｃａｌｔｏｃｏｒｒｅｃｔｉｏｎｔｈｅｄｉｓｔｏｒｔｉｏｎａｐｐａｒｅｎｔｒｅｓｉｓｔｉｖｉｔｙｓｅｃｔｉｏｎａｎｄｐｈａｓｅｓｅｃｔｉｏｎ．Ｉｎｔｈｉｓｐａｐｅｒ，ａｃｃｏｒｄｉｎｇｔｏｃｏｋｒｉｇｉｎｇＳｒｕｌｅａｎｄｕｎｉｔｅｄｔｈｅｔｗｏｔｙｐｅｍｏｄｅｌｏｆＣＳＡＭＴ，ｗｅｃｏｍｐｉｌｅｄｔｈｅｓｏｆｔｗａｒｅａｂｏｕｔｔｈｅｓｔａｔｉｃｃｏｒｒｅｃｔｉｏｎｏｆＣＳＡＭＴｄａｔａ．Ｗｅｈａｖｅｔｅｓｔｅｄｔｈｅｔｈｅｏｒｙｐａｔｔｅｒｎ．Ｉｔｓｈｏｗｓｔｈｅｍｅｔｈｏｄｉｓｍｏｒｅｅｆｆｅｃｔｉｖｅ．Ｋｅｙｗｏｒｄｓ：ＣＳＡＭＴ；ｃｏｋｒｉｇｉｎｇ；ｓｔａｔｉｃｓｈｉｔｆ

地球物理计算常用的插值方法-克里格法

克里格法（Kriging）是地统计学的主要内容之一，从统计意义上说，是从变量相关性和变异性出发，在有限区域内对区域化变量的取值进行无偏、最优估计的一种方法；从插值角度讲是对空间分布的数据求线性最优、无偏内插估计一种方法。

克里格法的适用条件是区域化变量存在空间相关性。

克里格法，基本包括普通克里格方法（对点估计的点克里格法和对块估计的块段克里格法）、泛克里格法、协同克里格法、对数正态克里格法、指示克里格法、折取克里格法等等。

随着克里格法与其它学科的渗透，形成了一些边缘学科，发展了一些新的克里金方法。

如与分形的结合，发展了分形克里金法；与三角函数的结合，发展了三角克里金法；与模糊理论的结合，发展了模糊克里金法等等。

应用克里格法首先要明确三个重要的概念。

一是区域化变量；二是协方差函数，三是变异函数一、区域化变量当一个变量呈空间分布时，就称之为区域化变量。

这种变量反映了空间某种属性的分布特征。

矿产、地质、海洋、土壤、气象、水文、生态、温度、浓度等领域都具有某种空间属性。

区域化变量具有双重性，在观测前区域化变量Z(X)是一个随机场，观测后是一个确定的空间点函数值。

区域化变量具有两个重要的特征。

一是区域化变量Z(X)是一个随机函数，它具有局部的、随机的、异常的特征；其次是区域化变量具有一般的或平均的结构性质，即变量在点X与偏离空间距离为h的点X＋h处的随机量Z(X)与Z(X+h)具有某种程度的自相关，而且这种自相关性依赖于两点间的距离h与变量特征。

在某种意义上说这就是区域化变量的结构性特征。

二、协方差函数协方差又称半方差，是用来描述区域化随机变量之间的差异的参数。

在概率理论中，随机向量X与Y 的协方差被定义为：区域化变量在空间点x和x+h处的两个随机变量Z(x)和Z(x+h)的二阶混合中心矩定义为Z(x)的自协方差函数，即区域化变量Z(x) 的自协方差函数也简称为协方差函数。

一般来说，它是一个依赖于空间点x 和向量h 的函数。

通信网络-可控音频大地电磁法的静态效应原理及其应用

可控音频大地电磁法的静态效应原理及其应用作者姓名：林道博班级专业：地球物理学指导教师:肖宏跃摘要可控源音频大地电磁法(CSAMT)是在大地电磁法（MT）基础上，为克服静态效应问题，利用静态校正来提高处理解释质量一种地球物理方法。

本论文从静校正方法原理出发，结合实际情况，建立地下均匀介质中的静态模型和地下层状介质中的静态模型，尝试使用相关系数和小波分析来区分静态效应现象和异常表象。

通过对静校正方法的研究，选择使用EMAP 滤波和法小波分析法进行静态校正。

在实际资料的处理中，通过静态校正处理，有效地消除了低阻静态效应，避免了将某些静态效应解释为低阻异常，同时还压制了高阻静态效应的影响，提高了中间低阻层的分辨能力，本文静校正方法得到了很好的实用效果。

关键字：CSAMT，静态效应，静态校正，EMAP 滤波法，小波分析法The principle of The Static shift of Controlled-source Audio-frequency Magnetotelluric And ApplicationAbstract:Controlled-sourceAudio-frequency Magnetotelluric (CSAMT) are developing from the Magnetotelluric(MT), for correct the static shift ，which is used to static correction in order to achieve a better quality CSAMT data.Starting from the principles of the static correction, in this paper we establish the static shift model in the underground layered stratum and in the homogeneous stratum, attempt to distinguish the static shift from anomalies by using the correlation coefficient or wavelet analysis. After doing the research studies of some kinds of static correction, two static correction methods are chosen in this paper: EMAP and the wavelet analysis method.In practical CSAMT data processing, we take the two metal diggings for example, there is a serious static shift in the collected data. After static correction, a good quality data was achieved. In conclusion, these static correction methods gained a good applicable effect.Key Words:CSAMT, static shift, static correction,EMAP, waveletanalysis method.目录第1章前言 (4)1．1 研究的目的及意义 (4)1．2研究现状 (5)1．3 研究范围及方法 (6)第2章 CSAMT的理论及静态效应原理 (7)2.1 CSAMT 方法的基本理论 (7)2.2 CSAMT 方法的优缺点 (9)2.3 静态效应的特征及小波识别静态效应 (10)2.3.1 静态效应的特征 (10)2.3.2 小波变换识别静态效应 (11)第3章静态效应校正的方法原理及特征分析 (14)3.1 EMAP 滤波法做静校正 (14)3.2 小波分析压制静态效应 (16)3.2.1 S.Mallat算法 (16)3.2.2 多尺度分离 (17)3.2.3 静态效应的多尺度压制 (20)3.3 静态效应的响应特征分析 (22)第4章静态效应压制实例 (24)4.1 理论模型 (24)4.2 实例 (25)结论 (30)致谢 (31)参考文献 (32)第1章前言1．1 研究的目的及意义可控源音频大地电磁法（CSAMT）是在大地电磁法（MT）的基础上，针对解决大地电磁法场源的随机性强和信号微弱而发展成的一种人工源频率域电磁测深方法。

基于格点型有限体积法的大地电磁二维正演

基于格点型有限体积法的大地电磁二维正演在地球物理领域，大地电磁法作为一种探测地下电阻率分布的重要手段，广泛应用于矿产勘探、地质调查、环境监测等领域。

而基于格点型有限体积法的大地电磁二维正演，作为一种计算方法，对于模拟和解释地下介质电阻率结构具有重要意义。

本文将深入探讨基于格点型有限体积法的大地电磁二维正演，并结合实际案例，对其应用进行详细分析。

1. 格点型有限体积法简介格点型有限体积法是一种数值模拟方法，用于求解偏微分方程。

它将求解区域离散化为多个小区域（或称为网格单元），通过在每个小区域内建立有限体积方程，将偏微分方程转化为代数方程组，从而进行数值求解。

在大地电磁二维正演中，格点型有限体积法能够准确描述电磁场在不同介质中的传播特性，为电阻率结构的模拟提供了有效手段。

2. 大地电磁二维正演原理大地电磁二维正演是指根据地球电磁法理论，通过计算建立在地球表面或空间中的电磁场分布情况，以得出地下介质电阻率结构的一种方法。

在二维正演中，主要考虑了地下介质的横向变化情况，将地下介质分割成多个区域进行离散化处理，然后利用格点型有限体积法对电磁场的传播过程进行数值模拟，最终得出在不同位置和时间的电磁场响应。

3. 应用案例分析以矿产勘探为例，通过大地电磁二维正演，可以模拟地下电阻率分布，进而推断出潜在的矿产资源分布情况。

在实际勘探中，通过布设电磁探测器，采集地下电磁场数据，并结合基于格点型有限体积法的二维正演模拟，得出地下电阻率结构图像，为矿产勘探提供了重要的依据。

4. 个人观点和理解基于格点型有限体积法的大地电磁二维正演在地球物理勘探领域具有重要的应用意义。

它不仅可以对地下介质进行精确模拟，还能够帮助解释大地电磁场数据，为地下结构的认识提供重要支持。

在未来的研究中，我认为可以结合机器学习和人工智能等技术，进一步提升大地电磁二维正演的精度和效率，从而更好地应用于实际勘探和监测中。

总结回顾：本文深入探讨了基于格点型有限体积法的大地电磁二维正演，介绍了其原理和应用，并结合个人观点进行了分析。

地质统计学多重指示克里格(MIK)方法原理及运用

图３扎拉金矿区地质特征
矿区内构造活动主要有２个变形阶段：①挤压和褶皱；②高角度逆冲断层。由于后期的强烈挤压构造活动，在２种岩性之间产生了平行于地层片理而近乎垂直的逆冲断层带。杂拉金矿化主要赋存于一个高角度逆冲断层带中，其东部巨大的火山沉积岩体被推覆在西部沉积物上。２种岩性的显著差异为随后的石英斑岩岩脉侵入和热液活动以及再活化提供了渗透空间。近乎ＳＮ走向的高角度逆冲断层构造控制了扎拉主矿化带的分布。金矿化赋存于蚀变石英斑岩内，在杂砂岩与碳酸盐，绢云母蚀变杂砂岩接触带之间。主矿化带南北延伸约５５０ｍ，宽度为２０～３５ｍ，走向３５０°～３６０°，向西陡倾斜（８０°～９０°）。区内与金矿化有关的主要热液蚀变类型是碳酸盐岩、绢云母、硅石、绿泥石或黄铁矿蚀变。矿化带多由呈斜梯平行状石英细脉状和浸染脉状黄铁矿方铅矿黄铜矿脉组成（图４）。矿化类型初步定为热液蚀变石英脉状类型［２］。５４
根据间隔约４０ｍ的横剖面进行金矿化解译和圈定，以０．３ｇ／ｔ为下限截值品位。选择该截值品位是因为它代表了捕获金矿化异常的截值品位。首先进行了概率为５０％的指示估算，然后将这个指示估算结果与钻孔品位和钻孔地质编录的信息相结合，建立最终的金矿化圈定标准，以尽量捕获所有的金矿化异常。图５为Ａｕ品位大于５ｇ／ｔ的平面截图，剖面图也显示矿体在品位上的明显差异性。地质统计学多重指示克里格（ＭＩＫ）方法适合于处理这类局部品位变异性较大的数据。２．４数据处理与统计分析
关键词多重指示克里格（ＭＩＫ）阀值选别开采单元（ＳＭＵ）体积变更技术金矿可采矿石储量
ＤＯＩ：１０．３９６９／ｊ．ｉｓｓｎ．１６７４６０８２．２０１９．０４．０１４
１理论方法
１．１指示克里格基本原理指示克里格是一种常用的非线性估值方法。非

01-静校正基础知识分析

0
x 地表高速层
CMP点
反射层
基本概念
低速带：存在于地表面的低速介质称为
低速带。低速带一般为近代沉积物或风化层。
降速带：低速带之下速度高于低速带的
介质称为降速带。降速带是低速带与高速层之间的过度带，它的岩性一般与低速带差异不大，主要由于压实作用或含水程度的不同造成的速度差异；但也有与低速带岩性不同的情况。降速带有些地方有而有些地方可能没有。
静校正工作实质上就是从地震角度对表层地质的系统研究
静校正基础知识
静校正工作的重要性工作职责、标准和要求静校正有关基础知识
静校正技术支持岗工作职责
1 负责本探区表层调查和静校正方法的确定及静校正技术设计和有关要求的编写工作； 2 负责组织静校正技术攻关和方法试验工作，解决本探区存在的静校正技术问题； 3 定期或不定期检查静校正基础资料和最终静校正成果；
4 负责确定使用的静校正软件资源，协调静校正软件资源中有关问题的解决；
5 协助技术领导作好本探区静校正方面的科研工作； 6 负责本探区静校正技术人员的培训和考核及聘用工作； 7 负责组织本探区静校正数据库的管理和完善。
静校正技术操作岗工作职责
1 按标准、设计及有关要求完成表层调查原始资料的评价、整理和解释； 2 定期或不定期检查表层调查野外采集质量及施工参数的正确性； 3 完成上级技术部门要求的静校正方法研究的试验工作； 4 按设计和有关要求完成静校正计算工作； 5 根据实际情况提出表层调查和静校正方法的调整和改进建议； 6 按标准和要求收集、整理静校正资料并上交地区经理部； 7 负责静校正所用软、硬件资源的日常维护、保养和管理，项目完成后交回发放部门。
利用反射波
静校正基础知识

克里格插值法

工程数学
工程数学
提出了如下的平稳假设及内蕴假设：提出了如下的平稳假设及内蕴假设：
{ 随机函数：随机函数：Z (u ), u ∈ 研究范围} ，其空间分布律不因平移而改变，即若对任一向量h，而改变，即若对任一向量，关系式
F ( z1 , z2 , ⋅⋅⋅; x1 , ⋅⋅⋅) = F ( z1 , z2 , ⋅⋅⋅; x1 + h, x2 + h, ⋅⋅⋅)
D(ξ ) = Var (ξ ) = E[ξ − E (ξ )] = E (ξ ) − E (ξ )2 22来自工程数学工程数学
（3）协方差）协方差是用来刻画随机变量之间协同变化程度的指标，协方差是用来刻画随机变量之间协同变化程度的指标，其大小反映了随机变量之间的协同变化的密切程度。大小反映了随机变量之间的协同变化的密切程度。
σ ij = Cov(ξ1 , ξ 2 ) = E[(ξ1 − E (ξ1 ) (ξ 2 − E (ξ 2 ) ] ））
= E (ξ1ξ 2 ) − E (ξ1 ) E (ξ 2 )
（4）相关系数）协方差是有量纲的量，与随机变量分布的分散程度有关，协方差是有量纲的量，与随机变量分布的分散程度有关，为消除分散程度的影响，提出了相关系数这个指标。消除分散程度的影响，提出了相关系数这个指标。
成立时，则该随机函数成立时，则该随机函数Z(x)为平稳性随机函数。为平稳性随机函数。这实际上就是指，无论位移h多大，两个维向量的随机变量多大，这实际上就是指，无论位移多大两个k维向量的随机变量
{ Z ( x1 ), Z ( x2 ),L , Z ( xk )} 和 { Z ( x1 + h), Z ( x2 + h),L , Z ( xk + h)}

kriging 方法

kriging 方法Kriging方法，又称克里格插值法，是一种常用于空间插值的统计方法。

它的主要目的是通过已知的数据点来估计未知位置的值，并给出估计值的可靠性信息。

在地理信息系统（GIS）和地质学领域，克里格插值法被广泛应用于栅格数据的插值和空间预测。

克里格插值法基于一个重要的假设，即空间上相近的点具有相似的属性值。

根据这个假设，插值方法通过计算距离权重来估计未知位置的属性值。

克里格插值法有多种变体，其中最常用的是简单克里格法和普通克里格法。

简单克里格法是克里格插值法的最简单形式，它假设空间上各点之间的距离权重与其距离成反比。

简单克里格法的估计结果仅依赖于最近邻的数据点，因此插值结果可能会出现较大的变化。

普通克里格法是一种改进的插值方法，它考虑了更多的数据点，并通过计算协方差来确定权重。

普通克里格法对距离较近的点赋予较大的权重，对距离较远的点赋予较小的权重。

通过对协方差进行插值，普通克里格法能够提供更准确的空间预测结果。

在使用克里格插值法之前，我们需要先进行数据的分析和预处理。

首先，我们要检查数据的空间分布情况，了解数据点之间的关系。

其次，我们要检查数据的属性值是否存在异常值或离群点。

如果存在异常值，需要进行数据清洗或者采用合适的处理方法。

最后，我们要选择合适的克里格插值方法和参数，以获得最佳的插值效果。

在进行克里格插值时，我们需要选择合适的变程参数和协方差函数。

变程参数决定了插值结果的平滑程度，较大的变程参数会产生较平滑的插值结果，而较小的变程参数则会产生较崎岖的插值结果。

协方差函数则用于计算不同距离下的权重，常用的协方差函数有指数型、高斯型和球型等。

除了简单克里格法和普通克里格法，还有一些改进的克里格插值方法，如克里格法的泛化版本——逆距离加权插值法（IDW）。

逆距离加权插值法通过计算数据点与插值位置之间的距离倒数来确定权重。

与克里格插值法相比，逆距离加权插值法对最近邻点赋予更高的权重，对较远的点赋予较小的权重。

石油物探2008年第47卷第1期～第6期总目次

唐湘蓉
贺振华
林
柏冠军张
地震资料反演砂岩孔隙度方法 … ……………………… …………………………… ……… 张亚敏
张书法
钱
利（３）１６
｜＿
・
２・
石
油
物
探
第４７卷
三维叠前转换波克希霍夫时间偏移方法研究 …………… …… ……… 王
孝
刘文卿廖前进
雍学善王尚旭
周祺
郑荣才赵正文
谢占安刘建武
李凤杰彭宇（７）７
付兴深刘苗张文林吴健生（３）８（９）８
叠前地震多属性反演在南堡油田火成岩研究中的应用 ……… ……… 刘淑华
克里格方法在大地电磁静校正中的应用 … ………………………… ……………………… 吴长祥
第４卷第６期７
２。年１月Ｏ８１
石
油
物
探
Ｖｏ＿７Ｎｏ６ｌ４，．
Ｎｏ．２Ｏｖ，Ｏ８
ＧＥＯＰＨＹＳＣＡＬＰＲＯＳＥＣＴＩＩＰＮＧＯＲＴＲＯＬＥＵＭＦＰＥ
石油物探
ＳＩＨＹｏＵＵＴＡＮＷ
２００８年第４卷７
准噶尔盆地南缘山前带变速构造成图 ………………………………… 杨江峰ＶＰ资料在转换波层位标定中的应用 ……………………… 程莉莉曹立斌Ｓ
（６１７）
波
臧殿光
李
全宁清香
（７）１２
（７）１９（８）１３
王成彬洪太元梁波吴战培
沈向存张学东杨冬梅凌振翔

静校正方法PPT课件教学提纲.ppt

第二类:信息源来自于正常生产地的初至信息
我们知道，正常生产炮的初至信息一般是直达波和近地表折射波，进人复杂山地以后，初至波信息变得十分复杂，除了上述两种类型波以外，可能还有透射波、反射波、反射折射波、折射反射波，以及多次折射波和多次折射反射波等。利用初至信息估算静校正量的方法为数众多，在生产中应用十分广泛，是十分重要的一类静校正量估算方法。
（6）为了适应多波多分量地震勘探技术发展的需
求，横波静校正方法的研究与应用，已逐渐提到议事日程上来。横波速度和频率均比纵波要低，在表层结构模型传播的过程中，必然表现出与纵波有其不同的特征，应当有适应横波特征的静校正量方法。完全采用与纵波相同的静校正量方法估算静校正量，其精度很难满足横波静校正的要求。另外。目前普遍采用纵波震源激发、在界面上发生波的转换，在地面接收横波。这就是说，射线向下是纵波，向上是横波，炮点是纵波校正量，接收点是横波校正量。我们必须研究与此相适应的静校正方法。
第一节概述
静校正是陆地地震资料常规处理流程中必不可少的一环。在我国西北地区，地表条件比较复杂，静校正问题尤为严重。目前地震勘探的重点主要在我国的西部, 在这些地区静校正问题严重制约着地震勘探的效果，解决好静校正问题具有重要的理论意义和实际意义。
静校正方法从所采用的信息源头来分大致可分为三类：
成果：
（1）近地表模型的调查与研究包括野外观测方法与室内处理方法。着重要解决的是折射波界面极不稳定的、初至波信息十分复杂的情况下，近地表模型的调查与建立的有效方法，提高模型与实际情况的吻合程度。
（2）在算法上，一个是走时层析反演技术的发展与完善；一个是非线ቤተ መጻሕፍቲ ባይዱ系统的引入和对非线性算法的研究；而另一个是迭代运算的效率。在地球物理方法上，探索多波多分量信息的运用，扩大用于静校正量估算所用的波的类型。

大气校正中用到的克里金插值原理

大气校正中用到的克里金插值原理下载提示：该文档是本店铺精心编制而成的，希望大家下载后，能够帮助大家解决实际问题。

文档下载后可定制修改，请根据实际需要进行调整和使用，谢谢！本店铺为大家提供各种类型的实用资料，如教育随笔、日记赏析、句子摘抄、古诗大全、经典美文、话题作文、工作总结、词语解析、文案摘录、其他资料等等，想了解不同资料格式和写法，敬请关注！Download tips: This document is carefully compiled by this editor. I hope that after you download it, it can help you solve practical problems. The document can be customized and modified after downloading, please adjust and use it according to actual needs, thank you! In addition, this shop provides you with various types of practical materials, such as educational essays, diary appreciation, sentence excerpts, ancient poems, classic articles, topic composition, work summary, word parsing, copy excerpts, other materials and so on, want to know different data formats and writing methods, please pay attention!一、引言大气校正是遥感图像处理中的重要环节，它能够消除大气造成的干扰，提高遥感图像的质量和准确度。

ArcGIS教程：克里金法的工作原理(参考模板)

ArcGIS教程：克里金法的工作原理克里金法是通过一组具有 z 值的分散点生成估计表面的高级地统计过程。

与插值工具集中的其他插值方法不同，选择用于生成输出表面的最佳估算方法之前，有效使用克里金法工具涉及 z 值表示的现象的空间行为的交互研究。

什么是克里金法?IDW(反距离加权法)和样条函数法插值工具被称为确定性插值方法，因为这些方法直接基于周围的测量值或确定生成表面的平滑度的指定数学公式。

第二类插值方法由地统计方法(如克里金法)组成，该方法基于包含自相关(即，测量点之间的统计关系)的统计模型。

因此，地统计方法不仅具有产生预测表面的功能，而且能够对预测的确定性或准确性提供某种度量。

克里金法假定采样点之间的距离或方向可以反映可用于说明表面变化的空间相关性。

克里金法工具可将数学函数与指定数量的点或指定半径内的所有点进行拟合以确定每个位置的输出值。

克里金法是一个多步过程;它包括数据的探索性统计分析、变异函数建模和创建表面，还包括研究方差表面。

当您了解数据中存在空间相关距离或方向偏差后，便会认为克里金法是最适合的方法。

该方法通常用在土壤科学和地质中。

克里金法公式由于克里金法可对周围的测量值进行加权以得出未测量位置的预测，因此它与反距离权重法类似。

这两种插值器的常用公式均由数据的加权总和组成：在反距离权重法中，权重λi 仅取决于预测位置的距离。

但是，使用克里金方法时，权重不仅取决于测量点之间的距离、预测位置，还取决于基于测量点的整体空间排列。

要在权重中使用空间排列，必须量化空间自相关。

因此，在普通克里金法中，权重λi 取决于测量点、预测位置的距离和预测位置周围的测量值之间空间关系的拟合模型。

以下部分将讨论如何使用常用克里金法公式创建预测表面地图和预测准确性地图。

使用克里金法创建预测表面地图要使用克里金法插值方法进行预测，有两个任务是必需的：•找到依存规则。

•进行预测。

要实现这两个任务，克里金法需要经历一个两步过程：•创建变异函数和协方差函数以估算取决于自相关模型(拟合模型)的统计相关性(称为空间自相关)值。

克里格方法(Kriging)

03 (h)= C(0) —C(h)
精选完整ppt课件
3
克里格法
01 Z(p)为区域Ω上随机过程，p∈Ω； Ω上有n个测点（样本点），
zi z(pi)在 p i处的测值，则 p 0 处的最优线性估计为
n
zˆ0 i zi i1
02 最小化非测点 p 0 处的估值方差 0 2E[z(0zˆ0)2],可推导出克里
2
基本概念
01 变差函数：Z(p)为一随机过程，Z(p)在p，p+h两点处的值之差的方差之半定义为Z（p）在p方向上的变差函数，记为
(h)1V[az((rp)z(ph)]
2 变差函数描述了区域化变量的空间结构性。 (h)只依赖于h。
02 协方差函数：随机过程Z(p) 在p1、p2处的两个随机变量Z(p1) 和Z(p2)的二阶混合中心矩，即 Cov{Z(p1), Z(p2)}=E[Z(p1)*Z(p2)]-E[Z(p1)]*E[Z(p2)]，记为 C(p1, p2) 整个区域中，Z(p)的协方差函数存在且相同，即只依赖于h Cov{Z(p),Z(p+h)} ≜C(h)；当h=0时，C(0)=Var{Z(x)}，x
n
i1
n j1
1 B'(hij)
B'(hij) k
0
精选完整ppt课件
6
优化测点分布的克里格方程组
由(h)=C(0)B(h)，可得 C(h)=C(0)(1-B(h))
设 ce(h)1B(h) ，则上式可表示为
c(h)c(0)ce(h)
令 c(0)e 将上述式子代入克里格方程组可得与C(0)无关的克里格方程组和克里格方差，如下
g(i)
，表明网格节点上的较大估值方差变大了，

用地震资料预测孔隙度的广义协克里格法

用地震资料预测孔隙度的广义协克里格法
孙树海;刘雪莹
【期刊名称】《地球物理学报》
【年(卷),期】1993(000)006
【摘要】在Ｄｏｙｅｎ工作的基础上，本文对简单协克里格法作了一些改进，提出了用地震资料预测孔隙度广义协克里格法，详细推导证明了其基本公式，并分别用３种方法预测了松辽盆地宋站地区扶余油层含气砂岩储层的孔隙度分布，井孔处的交互检验证明，广义协克里格法比线性回归法和简单协克里格法预测孔隙度的均方根误差分别小４０％和３０％以上。

【总页数】1页(P798)
【作者】孙树海;刘雪莹
【作者单位】不详;不详
【正文语种】中文
【中图分类】P631.44
【相关文献】
1.用地震资料预测孔隙度的简单协克里格法 [J], 孙树海;刘雪莹
2.用地震资料提高储层参数插值精度的新方法--3-D协克里格与主成分下沉法 [J], 马云潜;汪步紧;陈慧蓉
3.泛克里格法在三维地震资料变速成图中的应用 [J], 孙开峰;叶勇
4.利用地震资料对川东地区飞仙关鲕滩储层孔隙度的预测研究及效果分析 [J], 吴大奎;李亚林;陈晓超
5.三维地震资料储层参数建模方法研究--分形协同克里格法 [J], 鲁红英;肖思和;何建军;邓林
因版权原因，仅展示原文概要，查看原文内容请购买。

高分辨率去静态频率域大地电磁法获美国专利

高分辨率去静态频率域大地电磁法获美国专利
佚名
【期刊名称】《地质装备》
【年(卷),期】2011(12)3
【摘要】不久前，天津华北地质勘查局物探中心的科技发明项目——高分辨率去静态频率域大地电磁法，在继2009年获国家发明专利后，又顺利通过了美国专利局的审核，喜获美国国家专利。

【总页数】1页(P7-7)
【关键词】大地电磁法;美国专利;高分辨率;频率域;静态;国家发明专利;地质勘查局;科技发明
【正文语种】中文
【中图分类】P631.325
【相关文献】
1.基于频率-速度域多重信号分类的面波高分辨率频散成像方法 [J], 伍敦仕;孙成禹;林美言
2.频率域高分辨率地震波阻抗直接反演方法研究 [J], 王静波;陈祖庆;蒋福友;苏建龙
3.保鲜新技术带动农民种"开心果"/有机乌龙茶生产加工技术通过鉴定/"重盐碱地耐盐植物筛选及种植应用研究"通过鉴定/美国研制成功水果去核素/日本培育出补铁水稻新品种/羊草育种及扩繁栽培方法获专利 [J],
4.重茬种草莓不犯愁了/杨梅保鲜不再难/辣椒水当农药/变色金银花色变效益高/新式饲喂器获专利/家用碾米机、粉碎机获专利/美国培育出含抗癌物质的西红柿/鸡
蛋大小由人控 [J],
5.频率域线源大地电磁法有限元正演模拟 [J], 王若;王妙月;底青云
因版权原因，仅展示原文概要，查看原文内容请购买。

地球物理资料数字处理(第七讲、第八讲)

第二节静校正（Statics Correction）理论上，我们假定激发接收是一个平面，且地下介质是理想均匀的。

实际上的激发接收并不是在一个平面，而是一个起伏不平的，且地下介质也不是均匀的，特别是近地表的风化层形成了一个低速带。

因此，静校正的目的是消除地表起伏和低速带横向变化不均匀对地震记录产生的时差。

本节要点：●一次静校正原理；●剩余静校正地表一致性假设条件；●自动统计静校正原理；●自动统计静校正步骤。

一、基准面校正(Datum plane Statics Correction)如图(2-7)所示,h0为激点高程,hH为激发井深,hs为实际激发高程.hR为接收点高程,h1为激发点基准面以下低速带厚度,h2为接收点基准面以下低速带厚度。

图２－７地震波传播示意图设地震波在低速带中按垂直路径传播，则1．基准面以上的校正量(Above Datum plane Statics)00s Rh h V V + (2-1-1)2．基准面以下的校正量(Below Datum plane Statics)1122000()()h h h hV V V V -+- (2-1-2) 3．总静校正量(Total Statics)11220000()()s R s h h h h h hT V V V V V V∆=++-+- (2-2-3) 分析上式，若想求出静校正量，则必须测出激发点和接收点的高程以及低速带的速度和厚度等数据。

因此，静校正的精度取决于上述数据的精度。

二、剩余静校正(Residual Statics Correction)经过基准面校正（也称为一次静校正）后，所有的数据并未精确地校正到基准面上，仍然存在有一定的剩余静校正量。

这是因为我们很难获得精确的地下低速带的速度及厚度数据，由于存在剩余静校正量，使得多次迭加的结果质量下降，因此必须进行剩余静校正。

1．剩余静校正地表一致性假设(Surface Consistency Assumption Residual Statics Correction)１）同一炮点，地震波在低速带中的旅行时间与入射角无关，即认为地震波在低速带中是垂直传播的。