求解目标地波散射特性的方法研究

本文把求解半空间散射问题的FDTD技术与地波传播理论相结合，研究任意复杂目标的地波散射特性。入射地波设置在FDTD计算区域中的总场边界上并在散射场输出边界面上提取散射近场数据，然后利用等效及镜像原理，计算无地波衰减时的远区散射场，通过引入地波衰减因子，把该远区散射场转换为远区地波散射场。文中给出了方法验证例子以及一个较复杂目标的单站RCS计算结果。

海上舰船地波散射雷达截面的计算与分析

给出了海上目标地波散射雷达截面的计算方法．首先把目标的地波散射问题转化为求解理想导体平面上目标对垂直极化平面波的散射，利用时域有限差分方法得到所关心频段内的目标雷达散射截面；然后由地波传播理论求得包括随机起伏粗糙海面影响在内的地波传播衰减因子．从雷达方程出发，对海上舰船雷达散射截面的计算值与实测值进行对比，表明了该方法的可行性．

平直有耗地面上目标地波散射FDTD近-远场外推

由平直有耗地面上垂直电偶极子产生地波场的解析解和互易原理,本文导出了用FDTD方法处理目标地波散射问题时的近一远场外推公式。在数值结果验证了该公式的正确性之后,对海上一金属半球的地波散射场进行了求解,并以此与李清亮等(1997)提出的工程应用方法相对比,考察了工程应用方法的适用性。

Numerical simulation scattered imaging in deep mines

Conventional seismic exploration,mostly based on reflection theory,hardly has accurate imaging results for disaster geologic bodies which have small scale,steep dip,or complex structure.In this paper,we design two typical geologic models for analyzing the characteristics of scattered waves in mines for forward modeling by finite difference and apply the equivalent offset migration（EOM）and EOM-based interference stack migration methods to mine prospecting.We focus on the analysis of scatted imaging’s technological superiority to reflection imaging.Research shows：1）scattered imaging can improve fold and make the best of weak scattered information,so it shows better results than post-stack migration imaging and 2）it can utilize the diffraction stack migration method-based ray path theory for mine seismic advanced prediction,so it provides an new efficient imaging method for improving resolution of mine seismic exploration.

Physical mechanism of seismic attenuation in a two-phase medium

High-frequency seismic attenuation is conventionally attributed to anelastic absorption. In this paper, I present three studies on high-frequency seismic attenuation and propose that the physical mechanism results from the interference of elastic microscopic multiple scattering waves. First, I propose a new theory on wave propagation in a two-phase medium which is based on the concept that the basic unit for wave propagation is a nano- mass point. As a result of the elasticity variations of pore fluid and rock framework, micro multiple scattering waves would emerge at the wavelength of the seismic waves passing through the two-phase medium and their interference and overlap would generate high- frequency seismic attenuation. Second, I present a study of the frequency response of seismic transmitted waves by modeling thin-layers with thicknesses no larger than pore diameters. Results indicate that high-frequency seismic waves attenuate slightly in a near-surface water zone but decay significantly in a near-surface gas zone. Third, I analyze the seismic attenuation characteristics in near-surface water and gas zones using dual-well shots in the Songliao Basin, and demonstrate that the high-frequency seismic waves attenuate slightly in water zones but in gas zones the 160-1600 Hz propagating waves decay significantly. The seismic attenuation characteristics from field observations coincide with the modeling results. Conclusions drawn from these studies theoretically support seismic attenuation recovery.

Auto-pick first breaks with complex raypaths for undulate surface conditions

First-break picking is the key step in seismic data processing for surveying undulate surfaces, and directly infl uences the precision of near-surface modeling and effects of static corrections. The current first-break auto-picking methods may fail when the signalto-noise ratio（SNR） is low for seismic data in the undulate area, and require labor and time intensive manual picking. This study develops an improved super-virtual interferometry（SVI） method that combines multichannel and multidomain quality control（MMQC） techniques to achieve auto-picked first breaks. The improved SVI method extends the SVI application to enhance the SNR for near-surface scattered waves for the first time, which allows for the SVI method to adapt to first breaks with complex raypaths by linear combination of refractions and near-surface scattered waves. Methods of inverse and multidomain interferometry are developed to effectively enhance the virtual records extracted by the SVI method. The deconvolution filter for waveforms is used to increase resolution and reduce false picks, while the MMQC technique is designed to auto-correct false picks and increase the stability of auto-picking first breaks. The robust technique developed in this study enables stable processing of large 3D seismic datasets. Higher quality results are obtained using the approach presented in this paper to actual field data from the mountain areas in western China, when compared to some commonly used commercial software.

Computations on the Earth Terrain Return of Pulse Doppler Radar Altimeter

Concerning the PDRA (pulse doppler radar altimeter) designing and evaluation, owing to that the specifications of PDRA should be adaptively fixed according to the ETR (earth terrain return), and that in certain stages of product evaluation of PDRA which means the designations of PDRA are successful or not, the usage of ETR are indispensable, so the terrain return from spherical earth is critically important. A complete analytic derivation of the antenna shot section model of PDRA and the bright section model constrained by pulse emitted from antenna are given. Furthermore, the doppler effect mode and the earth terrain RCF (radar crossing factor) model are constructively analyzed. Then, the computing methodology on PDRA, which are used to compute the scattering power, scattering doppler spectrum, and the scattering signal, is studied. Besides, in order to check the correctness and efficiency of the algorithm, computing examples of ETR (earth terrain return) under the supposing premises are furnished. Finally, the conclusion is drawn that the models and algorithm are rational, the computational precise is satisfactory, the cost of computing time is low.

Seismic wave scattering inversion for fluid factor of heterogeneous media

Elastic wave inverse scattering theory plays an important role in parameters estimation of heterogeneous media.Combining inverse scattering theory,perturbation theory and stationary phase approximation,we derive the P-wave seismic scattering coefficient equation in terms of fluid factor,shear modulus and density of background homogeneous media and perturbation media.With this equation as forward solver,a pre-stack seismic Bayesian inversion method is proposed to estimate the fluid factor of heterogeneous media.In this method,Cauchy distribution is utilized to the ratios of fluid factors,shear moduli and densities of perturbation media and background homogeneous media,respectively.Gaussian distribution is utilized to the likelihood function.The introduction of constraints from initial smooth models enhances the stability of the estimation of model parameters.Model test and real data example demonstrate that the proposed method is able to estimate the fluid factor of heterogeneous media from pre-stack seismic data directly and reasonably.

Contrast patterns and trends of lapse rates calculated from near-surface air and land surface temperatures in China from 1961 to 2014

The near-surface lapse rate reflects the atmospheric stability above the surface.Lapse rates calculated from land surface temperature(γTs)and near-surface air temperature( γTa )have been widely used.However,γTs and γTa have different sensitivity to local surface energy balance and large-scale energy transport and therefore they may have diverse spatial and temporal variability,which has not been clearly illustrated in existing studies.In this study,we calculated and compared γTa and γTs at^2200 stations over China from 1961 to 2014.This study finds that γTa and γTs have a similar multiyear national average(0.53°C/100 m)and seasonal cycle.Nevertheless,γTs shows steeper multiyear average than γTa at high latitudes,and γTs in summer is steeper than γTa ,especially in Northwest China.The North China shows the shallowest γTa and γTs,then inhibiting the vertical diffusion of air pollutants and further reducing the lapse rates due to accumulation of pollutants.Moreover,the long-term trend signs for γTa and γTs are opposite in northern China.However,the trends in γTa and γTs are both negative in Southwest China and positive in Southeast China.Surface incident solar radiation,surface downward longwave radiation and precipitant frequency jointly can account for 80%and 75%of the long-term trends in γTa and γTs in China,respectively,which provides an explanation of trends of γTa and γTs from perspective of surface energy balance.