星载SSM/Ⅰ微波遥感渤海海冰的辐射特征分析

根据我国渤海区域星载微波遥感ＤＭＳＰＳＳＭ／Ｉ１９９４～１９９６年的观察数据的研究,指出现有的星载微波遥感海冰的数据标定和验证ＴＢＩ和ＯＳＩ算法是不适用的。用２２和８５ＧＨｚ的散射指数ＳＩ和１９ＧＨｚ的极化指数ＰＩ,讨论了渤海区域海冰的主要辐射特征,用矢量辐射传输理论数值模拟了ＳＩ、ＰＩ与海冰覆盖及其厚度变化的定量关系,用ＳＩ－ＰＩ图规划了渤海区域性海冰的识别方法,给出了渤海海冰的空间分布及其辐射特征的时间序列变化。     

Monitoring regional sea ice of the Bohai Sea by SSM/I scattering index

SSM/I data on the Bohai Sea of China from 1994 to 1996 have been studied.The sea ice algorithms of calibration/validation such as the brightness temperature index [AES/YORK TBI (Hollinger et al.,1992)] and ocean scattering index [OSI (Ferraro et al.,1996, The Bulletin of the American Meteorological Society, 77(5), 891-905)] are not applicable to the detection of regional sea ice of mid-latitude Bohai area.This paper presents a method of the scatfeting index SI and polarization index PI to identify regional Bohai Sea ice.Numerical results of vector radiative transfer for a model of a layer of sea ice over ocean simulate the relationship between SI/PI and sea ice depth, and others.It is employed to categorize the SI/PI indexes for the detection of sea ice.     

郯庐断裂带中南段地壳介质散射强度

本研究运用远震扰动场方法,利用134个数字地震观测台记录的10次远震记录对郯庐断裂带中南段地壳介质散射强度进行研究.结果表明,研究区上、下地壳介质散射强度呈现强烈的横向非均匀性特征,上、下地壳介质散射强度空间分布特征基本一致,散射强度高值区主要分布在华北断块区及秦岭—大别褶皱系的东部地区,散射强度低值区主要分布在下扬子断块区,散射强度峰值区主要分布在六安、莒南、嘉山、烟台等地.地壳介质散射强度与构造块体具有较强的相关性,构造块体不同,地壳介质散射强度则存在一定的差异;沿郯庐断裂带地壳介质散射强度呈现分段性特征,可分为潍坊至临沂段、临沂至嘉山段以及嘉山至九江段三段;地震活动与地壳介质散射强度有较强的相关性,地震多发生在地壳介质散射强度梯度带上,郯庐断裂带的潍坊至临沂段、嘉山至九江段地壳介质散射强度梯度较大,此段中小地震活动较强,而临沂至嘉山段地壳介质散射强度梯度小,中小地震活动相对较弱.地壳介质散射强度与构造块体、断裂构造及地震活动之间的相关性,除反映了深部介质物性可能存在差异外,也可能与深部介质物理形态、物质运移以及深部应力应变环境的变化有关.     

稀疏反演多次波去除策略与效果分析

多次波去除是海上勘探中的一项关键技术,在一定程度上决定着最终成像结果的质量.表面相关多次波去除方法(Surface-Related Multiple Elimination,SRME)是一种基于反馈迭代模型的数据驱动类方法,该方法利用地震数据本身预测多次波,进一步对其进行匹配相减获得一次波结果,这往往导致与多次波重叠的一次波被错误地去除.稀疏反演一次波估计方法(Estimation of Primaries by Sparsity Inversion,EPSI)在原理上同样基于反馈迭代模型,通过求解一次波和表面相关多次波总残差最小化问题直接获得一次波.该方法无多次波的匹配相减过程,可减少与多次波发生重叠的一次波产生损失.另外,EPSI可以很好地解决SRME无法解决的近道缺失问题.理论数据和实际资料处理对比表明SRME对数据完整性以及一次波和多次波的分布情况有较强的依赖性,而EPSI几乎不受这两个因素的影响.EPSI相比SRME方法计算成本高,尚未在工业界广泛应用,但随着高性能计算的普及,该方法的实用性有望得到突破.     

SPECTRAL STRUCTURE OF VELOCITY INHOMOGENEITY OF CRUST MEDIUM BELOW THE SOUTHEASTERN MARGIN OF TIBETAN PLATEAU AND ITS ADJACENT REGIONS

In this paper we present results of spectral structure of crustal velocity inhomogeneity beneath the southeastern margin of Tibetan plateau and its adjacent region based on the S wave envelope broadening algorithm. The spectral structure of 8~16Hz band is selected to analyze the special character of crustal inhomogeneity and discuss the correlation between strong earthquakes and inhomogeneities. The result shows that strong and complex inhomogeneities of crustal medium are found in the southeastern margin of Tibetan plateau and its adjacent region. In the upper part of upper crust, the strong and small scale inhomogeneities are imaged in the Longmengshan fault zone and the north of the Anninghe fault zone, the weak and large scale inhomogeneites are imaged in the section from Huolu to Daofu of Xianshuihe fault zone and the south of the Anninghe fault zone. In the lower part of upper crust, strong inhomogeneites are found in the Longmengshan fault zone, Lianfeng fault zone, the north of the Anninghe fault zone and the sections from Huolu to Daofu of the Xianshuihe fault zone, weak inhomogeneites are found in the section from Daofu to Kangding of Xianshuihe fault zone. In the middle crust, strong inhomogeneities are observed in the section of the Baoxing to Dujiangyan, the Baoxing to Kangding, and Kangding to Shimian, and weak inhomogeneities are observed in the northwestern section from Huolu to Kangding, and the Lianfeng fault zone. Comparing the medium inhomogeneities with the location of the strong earthquakes, our results suggest existence of high correlation between them. Strong earthquakes are often located in the transitionary zone between the strong and the weak inhomogeneities. The spatial distribution of the strong and the weak medium inhomogeneities may be related to the broken medium from the strong movement of geological tectonic and the heat flow upwelling along active faults induced by frequent tectonic and volcanic activity.     

The influence of errors in the source wavelet on inversion‐based surface‐related multiple attenuation

We investigate the influence of source wavelet errors on inversion‐based, surface‐related multiple attenuation, in order to address how the inverted primary impulse response, estimated primaries, and predicted multiples are affected by the estimated wavelet. In theory, errors in estimated wavelet can lead to errors in the upgoing waves. Because of smoothness and the band‐limitedness characteristics of the estimated wavelet, errors in the upgoing waves are usually not white and random. Theoretical analysis and two synthetic examples demonstrate that (i) when the overall amplitude scalar of the estimated wavelet is underestimated, the inversion of the primary impulse response suffers from instability, which will distort the estimation of primaries, and (ii) when the wavelet is overestimated, the estimated primaries will simply mimic the recorded upgoing waves. Nevertheless, the quality of the estimated primaries in the region above the first‐order, water‐bottom multiples is independent of the estimated wavelet. Synthetic results illustrate that inversion‐based, surface‐related multiple attenuation with a known wavelet is stable, since slight inaccuracy in amplitude spectrum and/or phase spectrum of the given wavelet or the corresponding upgoing waves will not lead to considerable deviation in the waveforms of the inverted results from those of the references. Furthermore, shot‐to‐shot wavelet variations, with maximum amplitude difference of 5% and maximum phase difference of 10°, create just slight artefacts in both the inverted primary impulse response and the estimated primaries. Moreover, the sensitivity test of estimation of primaries by sparse inversion method involving wavelet estimation shows that this method can stably and alternately update the wavelet and the primary impulse response; however, different choices of the initial wavelet can lead to different final inverted results.     

全球浅源巨大地震序列统计特征

系统研究了1976年以来全球58次M_W≥7.8浅源地震序列的统计特征。结果显示:(1)在58次巨大地震中,板间地震45次,板内地震13次,板内地震强度低于板间地震。(2)74.1%的板间地震为逆断层错动,61.5%的板内地震为走滑型错动。(3)58次地震序列中,82.8%为主-余型,17.2%为多震型;与5级以上地震序列不同,巨大地震没有孤立型,其余震比较活跃;板内地震中,多震型占7.7%,而板间地震中多震型占20%。(4)对于主-余型序列,75%的主震与最大余震的震级差为1.0～2.0级;震级差与主震震源错动类型有关,走滑型的震级差明显大于逆冲型;68%的最大余震发生在主震后3天内,其次为10天左右与1个月左右; 49%的D_(max-aft)(最大余震震中与主震震中间的距离)不超过余震区长轴的1/3,31%的D_(max-aft)为余震区尺度的1/3～1/2;最大余震的发生时间、最大余震震中与主震震中间的距离同主震断层错动类型间的关系不明显。(5)应用ETAS模型计算了46个序列参数后发现,b值、p值和a值均呈Beta分布,b值平均为1.164±0.211,p值平均为1.559±0.412,a值平均为1.673±0.911; p值和a值分布分散;对于不同的序列类型、震源错动类型及板内、板间地震,b值差异不显著;逆冲型序列p值明显大于走滑型和正断层型;板间地震序列a值明显小于板内地震;逆冲型序列a值明显小于走滑型和正断层型;这表明,与板内地震相比,板间地震具有较强的"余震激发余震"的能力;逆冲型破裂虽然会导致序列衰减较快,但触发次级余震的能力相对较强。(6)逆冲型巨大地震余震区长轴L的对数与主震震级M_W间的拟合关系式为lg L=(-1.399±0.306)+(0.470±0.037) M_W。     

Bayesian Markov chain Monte Carlo inversion for anisotropy of PP- and PS-wave in weakly anisotropic and heterogeneous media

A single set of vertically aligned cracks embedded in a purely isotropic background may be considered as a long-wavelength effective transversely isotropy (HTI) medium with a horizontal symmetry axis. The crack-induced HTI anisotropy can be characterized by the weakly anisotropic parameters introduced by Thomsen. The seismic scattering theory can be utilized for the inversion for the anisotropic parameters in weakly anisotropic and heterogeneous HTI media. Based on the seismic scattering theory, we first derived the linearized PP- and PS-wave reflection coefficients in terms of P- and S-wave impedances, density as well as three anisotropic parameters in HTI media. Then, we proposed a novel Bayesian Markov chain Monte Carlo inversion method of PP- and PS-wave for six elastic and anisotropic parameters directly. Tests on synthetic azimuthal seismic data contaminated by random errors demonstrated that this method appears more accurate, anti-noise and stable owing to the usage of the constrained PS-wave compared with the standards inversion scheme taking only the PP-wave into account.     

弹性波被动源数据一次波估计方法研究

与声波被动源数据相比,弹性波被动源数据含有更加丰富的地下地质信息.弹性波被动源数据合成的虚拟炮记录中,同样存在表面相关多次波,影响数据的后续处理以及对地下地质构造的准确判断.对弹性波被动源数据合成的虚拟炮记录进行多次波的预测与匹配相减,将比声波数据更加困难.本文提出弹性波被动源数据一次波估计方法,针对各向同性介质条件下的弹性波被动源数据进行研究,依据声波被动源数据一次波估计方法的研究思路,推导弹性波被动源数据一次波估计方法的理论过程,充分考虑弹性波被动源数据所包含的波场成分,包括纵波、横波以及在传播过程中产生的转换波.在传播过程中,由于能量的损失,略去多次传播的高阶项,进而给出相应的目标函数.利用数据的稀疏性,结合L_1范数约束的最优化问题求解方法,避免对弹性波虚拟炮记录进行多次波的预测与匹配相减,解决多次波影响问题.本文通过合成的弹性波被动源数据,进行一次波估计的测试,并与主动源和合成的虚拟炮记录进行对比,验证本文提出方法的有效性与准确性.     

A Method for Calculating Bottom Backscattering Strength Using Omnidirectional Projector and Omnidirectional Hydrophone

A method of obtaining bottom backscattering strength by employing an omnidirectional projector and omnidirectional hydrophone is proposed. The backscattering strength is extracted from monostatic backscattering data. The method was adopted in an experiment conducted in the South Yellow Sea of China. The seafloor surface was relatively smooth and covered by a small quantity of shell fragments, as observed through a digital camera system. Sampling data showed that the main component of the sediment at this experimental site was fine sand. In this paper, we detail the calculation method. Preliminary results of backscattering strength as a function of grazing angle(20?–70?) in the frequency range of 6–24 kHz are presented. The measured backscattering strength increased with the grazing angle and changed more rapidly at large grazing angles(60?–70?). A comparison of the data at different frequencies reveals that the measured backscattering strength substantially rises with the increase of acoustic frequency. A fitting curve of Lambert's law against the measured data shows that the backscattering strength deviates from Lambert's law at large grazing angles.