地震约束下密度反演.docVIP

地震约束下密度反演 　　摘要：广义线性联合反演对其模型的先验了解依赖性很大，只有在给定初始猜测接近真实地球物理模型时才可以求出有意义的解估计。地震波旅行时可以认为是目前最精确的地球物理探测量，本文在地震波旅行时的限制和纵波速度的约束下，利用常规最小二乘线性反演进行密度反演。利用该方法得到的反演结果比较接近真实值，作为先验信息可初始化重震迭代联合反演的模型参数。并且在速度参数接近准确的情况下，该反演结果比单由速度密度关系式得到的结果更接近真实值。 　　关键字：密度反演；旅行时限制；速度约束 　　Density Inversion under the Restriction of Seismic Data 　　Abstract: Joint inversion depends on much of the prior understanding of the model. Only when the given initial guess is close to the real physical model we can find out meaningful solution. We think that seismic waves travel time is the most accurate geophysical detecting quantity. Based on the restrictions of seismic waves travel time and the constraints of p-wave velocity, we can finish the density inversion using conventional least square linear inversion. The inversion results we get from the method as the priori information can initialize the parameters of the model used in iterative joint inversion of seismic and gravity data. When the velocity parameters are close to the accurate, the inversion results are close to the real value. 　　Key Words: density inversion; restrictions of travel time; constraints of velocity 　　 　　中图分类号：P315文献标识码：A文章编码： 　　 　　0 引言 　　地球物理反演是在地球物理学中利用地球表面观测到的物理现象推测地球内部介质物理状态的空间变化及物性结构的一个分支[1]。地球物理反演可分为单一地球物理现象的反演和多种地球物理现象的联合反演[2]。把单一地球物理资料的反演变为反映各种不同地球物理场或地球物理资料的联合反演, 是当今地球物理资料反演的总趋势, 也是发展的必然[3]。联合反演是地球物理综合研究的主要方法之一[4]。对于缺乏测井信息的深部探测，地震信息与非震信息的联合更是显得尤为重要。 　　广义线性反演是目前联合反演中应用最广泛的反演方法，但是广义线性联合反演对其模型的先验信息依赖性很大，只???在给定初始猜测接近真实地球物理模型时才可以求出有意义的解估计[1]。要想使模型参数良好的收敛于真实值，必须取得接近真实的模型参数。地震波旅行时可以认为是目前最精确的地球物理探测量，本文利用地震波旅行时的限制对地下介质进行时域区块结构划分，结合用资料处理过程中获得的速度信息得到空间域区块结构分布，并在速度-密度关系约束下进行密度反演。本方法可作为重震联合反演的重要组成部分，也可作为联合反演问题研究的重要入口。刘皓等人也曾在模型分块处理的基础上做了地震约束的重力反演，但他是在由速度-密度关系得出重力模型的基础上进行的参数调整[5]。 　　1 多边形重力正演 　　 　　图1 多边形截面 　　Fig.1 Polygon section 　　人们通常使用Parker法进行异常正演计算，但是Parker法依赖于接近层状分布的地下介质模型，而不适用于块体或复杂结构异常计算。因此，本方法利用多边形截面法计算横截面为任意形状二度体的重力异常，如图1所示。首先计算点O与多边形各边组成的多个三角形，，···，所引起的重力异常，并以顺时针给各三角形赋以正、负号。 　　 (1) 　　式中：。整个二度体在O点引起的重力异常，即横截面，，···，等N个二度