Improving resolution of gravity data with wavelet analysis and spectral method

Gravity data are the results of gravity force field interaction from all the underground sources. The objects of detection are always submerged in the background field, and thus one of the crucial problems for gravity data interpretation is how to improve the resolution of observed information. The wavelet transform operator has recently been introduced into the domain fields both as a filter and as a powerful source analysis tool. This paper studied the effects of improving resolution of gravity data with wavelet analysis and spectral method, and revealed the geometric characteristics of density heterogeneities described by simple shaped sources. First, the basic theory of the multiscale wavelet analysis and its lifting scheme and spectral method were introduced. With the experimental study on forward simulation of anomalies given by the superposition of six objects and measured data in Songliao plain, Northeast China, the shape, size and depth of the buried objects were estimated in the study. Also, the results were compared with those obtained by conventional techniques, which demonstrated that this method greatly improves the resolution of gravity anomalies. Translated from Progress in Geophysics, 2007, 22(1): 112–120 [译自: 地球物理学进展]     

Seismic attenuation of coda waves in the eastern region of Cuba

Cuba's seismic attenuation had never been studied in detail. In this paper we present the results of the research on the seismic attenuation of Cuba's eastern zone based upon the information collected by the seismological Cuban network from 1998 to 2003. 581 earthquakes were selected from the Cuban catalogue to make this study. All of them, recorded by at least three seismic stations, had their epicenters located in the eastern Cuban region (19.3–22 N, 79–73 W), epicentral distances between 15 km and 213 km, their coda duration magnitudes ranging from 2 and 4.1 and their focal depths reaching up to 30 km. The seismic wave attenuation was studied using coda waves. The single scattering method proposed by Sato in 1977 was applied, the attenuation and frequency dependency for different paths and the correlation of the results with the geotectonics of the region are presented in this paper.The mean Q_c value calculated was Q_c = (64 ± 2)f^{0.84 ± 0.01}. The relatively low Q₀ and the high frequency dependency agree with the values of a region characterized by a high tectonic activity. The Q_c values of seven subregions of eastern Cuba were calculated and correlated with the geology and tectonics of the area.     

基于景观分区的植被类型信息提取

利用陆地卫星多光谱遥感数据,在研究草被、针叶林、阔叶林、耕地等不同植被类型波谱特征的基础上,引入景观分区信息,进行空间特征与波谱特征相结合的遥感植被类型信息提取,其精度和效果优于单纯依赖植被指数的分类方法,能够适用于研究区植被类型的详细划分。     

光谱与纹理信息复合的土地利用/覆盖变化动态监测研究

及时、准确地动态监测地球表面特性对于掌握人类和自然现象之间的关系和相互作用是非常重要的,并为决策的制定奠定基础。传统卫星遥感的土地利用/覆盖变化动态监测方法基本上可分为光谱直接比较法和分类比较法两大类,但两类方法多以光谱信息为基础来提取土地利用变化信息,而忽略纹理等空间信息。本研究中,基于变化向量分析方法,将光谱与纹理两种信息复合计算变化强度,并采用支持向量机法提取变化/非变化信息,通过监督分类确定变化区域内的土地利用/覆盖类型的转移方向,完成土地利用/覆盖动态监测。最后,利用两期TM数据,对海淀区1997—2004年进行土地利用/覆盖变化动态监测,以验证该方法。该方法较分类后比较法在一定程度上减少误差积累,降低了错误类型转化,提取的变化信息总精度达到93.1%,Kappa为0.862,比利用光谱信息双窗口变步长的变化向量分析方法提取出土地利用/覆盖变化信息的精度有一定的提高(总体精度为90.2%,Kappa为0.804)。纹理信息与光谱信息复合,能够更大拉开变化/非变化信息之间的距离,有利于动态变化信息的提取,是该方法能够有效提取变化信息的关键所在。     

Study on Inelastic Attenuation and Source Parameters of the Shidian Swarm in 2001

Based on horizontal-component digital seismograms recorded on 6 stations of the Yunnan Regional Digital Network, we inversed the inelastic attenuation in the source region of the Shidian swarm using the Atkinson method and the site responses of the 6 stations 200km around the Shidian epicenters using the Moya method. The observational seismic waveform data were corrected by removing the propagation, instrument and site effects before the source parameters of the Shidian swarm in 2001 were determined using genetic algorithms. The results are as follows: (1) There is a linear relation between seismic moment and local magnitude. The seismic moment is between 1×10~ 12 and 10~ 14 N·m. The rupture radius of the seismic focus varies from 157m to 973m. The seismic moment and the rupture radius maintain a linear correlation. (2) The corner frequency increases as the seismic moment decreases. Based on the expression between corner frequency and seismic moment using least squares fitting, we can obtain the estimated value of the corner frequency. The time-varying value of the calculated corner frequency minus the estimated corner frequency shows that there were continuous high and low anomalies before the strong aftershocks. (3) The seismic stress drop is in the range of 0.07～1.55MPa. The stress drop seems independent of the local magnitude. The variation of stress drops is high before the occurrence of the strong aftershocks. (4) The depth of aftershocks is mostly in a range from 5km to 10km, which means that energy release of aftershocks is mainly concentrated in this range of depth.     

Attenuation of P and pP waves in Vrancea area – Romania

The seismic attenuation in the Vrancea region (Romania) is investigated from teleseismic recordings of P and pP waves during the four major, intermediate-depth Romanian events that occurred since the onset of digital instrumentation. Most stations are located in Canada and in the United States, being equipped with a variety of sensors, especially short-period ones. The amplitude spectral ratio method is used, assuming no frequency dependence of the Q _P factor in the range 0.2–2 Hz. No apparent correlation between the derived attenuation value and the type of recording sensor is observed. Lateral variations of the attenuation are obtained, with a very low Q _P area (values down to 33) located in the northwestern part of the Vrancea seismogenic volume. For the stations with different azimuth angles in relation to the epicentral area, Q _P values routinely exceed 200. Most likely, the low attenuation values are related to an upwelling mantle material located immediately beneath the crust, but limited in depth to at least 100 km.     

An inversion of site response and Lg attenuation using Lg waveform

Based on spectral ratio method, a joint inversion method was used to obtain parameters of Lg wave attenuation and site response. The inversion method allows simple and direct (two-parameter) determination of Lg wave attenua- tion and site response from sparse spectral data, which are not affected by radiation pattern factor and different response of same instrument after geometrical spreading. The method was used successfully for estimating site re- sponse of stations of Zhejiang Seismic Network and measuring Lg wave attenuation. The study is based on 20 earth- quakes occurred in northeast of Taiwan with magnitude MS5.0~6.7 and 960 seismic wave records from 16 stations in Zhejiang area from 2002 to 2005. The parameters of site response and Lg attenuation were calculated with a fre- quency interval of 0.2 Hz in the range of 0.5 Hz to 10 Hz. Lg wave attenuation coefficient corresponding to U-D, E-W and N-S components are γ ( f )=0.00175 f 0.43485, γ ( f )=0.00145f 0.48467 and γ ( f )=0.0021f 0.41241, respectively. It is found that the site response is component-independent. It is also found that the site response of QIY station is significant above the frequency of 1.5 Hz, and that the site response of NIB station is low for most frequency     

祁连山中东段地区非弹性衰减系数、震源参数和场地响应研究

利用甘肃数字地震台网的波形记录,基于遗传算法,首先用Atkinson方法计算非弹性衰减系数,在此基础上用Moya方法同时反演震源参数和场地响应。反演结果表明,祁连山中东段地区非弹性衰减Q值与频率的关系为:Q(f)=687.9f0.460;除了湟源台(HYQ)外,其余台站的场地响应均有不同程度的放大效应;研究区地震的应力降范围为2~40 bar,优势分布在2~20 bar;震源参数拐角频率与地震矩、地震矩与震级之间存在依赖关系。     

用HV谱比法推测银川市区沉积层基底分布特征

在缺乏其他资料的情况下,地脉动常常用做确定地震动对地表介质产生影响效应的简便手段。本文以银川市区为研究对象,根据地震勘探资料有目的的选取5条测线在银川市区进行地脉动观测,观测点共141个。用Nakamura提出的HV谱比法解析观测数据,在给定地下构造剪切波速度(VS)的基础上,由H/V比值曲线中的基本频率推测盆地深层地下构造。通过理论计算进一步验证了解析结果。结合前人地震勘探资料推测了银川市区脉动基本频率分布特征,并给出2D和3D银川市基底构造变化情况,对比分析表明,HV谱比法得到的脉动基本频率值可以反映银川盆地基底变化情况。