Theoretical Study on Stochastic Modeling of Combined Gravity-Magnetic-Electric-Seismic Inversion and Its Application

As gravity field,magnetic field,electric field and seismic wave field are all physical fields,their object function,reverse function and compound function are certainly infinite contiuously differentiable func-tions which can be expanded into Taylor (Fourier) series within domain of definition and be further reduced in-to solving stochastic distribution function of series and statistic inference of optimal approximation,This is the basis of combined gravity-magnetic-electric-seismic inversion built on the basis of separation of field and source gravity-magnetic difference-value(D-value)trend surface,taking distribution-independent fault sys-tem as its unit,depths of seismic and electric interfaces of interests as its corresponding bivariate compound re-verse function of gravity-magnetic anomalies and using high order polynomial(high order trigonometric func-tion)approximating to its series distribution,The difference from current dominant inversion techniques is that,first,it does not respectively create gravity-seismic,magnetic-seismic deterministic inversion model from theoretical model,but combines gravity-magnetic-electric-seismic stochastic inversion model from stochastic model;second,after the concept of equivalent geological body being introduced,using feature of independent variable of gravity-magnetic field functions,taking density and susceptibility related to gravity-magnetic func-tion as default parameters of model,the deterministic model is established owing to better solution to the con-tradictioc of difficulty in identifying strata and less test analytical data for density and susceptibility in newly explored area;third,under assumption of independent parent distribution,a real modeling by strata,the prob-lem of difficult plane closure arising in profile modeling is avoided,This technology has richer and more detailed fault and strata information than sparse pattern seismic data in newly explored area,successfully inverses and plots structural map of Indosinian discontinuty in Hefei basin with combined gravity-magnetic-electric-seismic inversion,With development of high precision gravity-magnetic and overall geophysical technology,it is certain for introducing new methods of stochastic modeling and computational intelligence and promoting the develop-ment of combined gravity-magnetic-electric-seismic inversion to open a new substantial and promoting the develop-ment of combined gravity-magnetic-electric-seismic inversion to open a new substantial path.     

地震波传播的褶积微分算子法数值模拟

为了解决地震波场数值模拟中的速度与精度的匹配及局部信息与全局信息的耦合等问题, 该文借助新推出的基于Forsyte广义正交多项式的褶积微分算子, 将计算数学中的Forsyte多项式褶积微分算子应用到地震波传播的数值模拟中.复杂非均匀介质模型数值模拟结果说明了该方法的可行性和优越性.该方法同时具有广义正交多项式褶积微分算子的高精度和有限差分短算子算法的高速度.通过对算子长度的调节及算子系数的优化, 可同时兼顾波场解的全局信息与局部信息.      

建立地震发射台的建议

目前地球上有成千上万个地震台，但它们都是接收地震波的地震台。建议建立一种新型的地震台，这种地震台用人工的方法不断地向外发射地震波，照亮需要研究地区的地下结构，研究特定的科学问题。建立地震发射台的核心问题是选择合适的人工震源，近年来出现的绿色环保震源，以及现代信息处理技术的发展，使得建立照亮10万km2(深至Moho面)地区地下结构的地震发射台已经完全是可能的。     

P波地震图上的地面反射波与S—P转换波

辐射脉冲约１Ｓ或更短持续时间的震源理论研究表明：对于假定为实际地壳结构的双力偶震源，短周期Ｐ波地震图最突出的震相当是Ｐ，ＰＰ和ＳＰ波。已注意到显示Ｐ波与地面反射的某些地震图，并可以充分模拟，然而，由浅地震观测到的这种地震图数目很小，通常，观测到的地震图要私很简单，即它们显示出突出的初至，而很少有其他信息，要么地震图很复杂，即显示振幅大于或等于初至振幅的３０Ｓ以上的震相，简单和或复杂地震图上缺乏清晰     

The inversion of anelastic coefficient,source parameters and site respond using genetic algorithm

It gradually becomes a common work using large seismic wave data to obtain source parameters, such as seismic moment, break radius, stress drop, with completingof digital seismic network in China (Hough, et al, 1999; Bindi, et al, 2001). These parameters are useful on earthquake prediction and seismic hazard analysis.Although the computation methods of source parameters are simple in principle and the many research works have been done, it is not easy to obtain the parameters accurately. There are two factors affecting the stability of computation results. The first one is the effect of spread path and site respond on signal. According to the research results, there are different geometrical spreading coefficients on different epicenter distance. The better method is to introduce trilinear geometrical spreading model (Atkinson, Mereu, 1992; Atkinson, Boore, 1995; WONG, et al, 2002). In addition, traditional site respond is estimated by comparing with rock station, such as linear inversion method (Andrews, 1982), but the comparative estimation will introduce some errors when selecting different stations. Some recent research results show that site respond is not flat for rock station (Moya, et al, 2000; ZHANG,. et al, 2001; JIN, et al, 2000; Dutta, et al, 2001). The second factor is to obtain low-frequency level and corner frequency fromdisplacement spectrum. Because the source spectrum model is nonlinear function,these values are obtained by eye. The subjectivity is strong. The small change of corner frequency will affect significantly the result of stress drop.     

地震波在湿颗粒介质中的传播

本文利用拉格朗日运动方程建立了湿颗粒介质中的波场方程组,由此研究了平面简谐波在其中的传播。文中还导出了湿颗粒介质中的支承拉梅系数和耦合弹性系数的表示式。     

横向非均匀介质中弹性地震波传播的数值模拟

应用混合变量弹性动力学方程和线性常微分方程组的矩阵指数解法，将层状介质中广泛应用的弹性波传播矩阵解法推广至横向非均匀介质，给出了一种可计算复杂地质体中弹性波传播的广义传播矩阵数值解法。该方法可模拟任意震源及所产生的各种体波、面波，数值结果表明具有很高的计算精度。