基于有效邻域波场近似的起伏地表保幅高斯束偏移

随着我国陆上地震勘探向复杂地表探区的转移,高精度、适应性强的地震成像方法在地震资料的处理、解释及后续属性分析、储层预测中具有重要意义.本文基于有效邻域波场近似理论发展了一种成像精度更高且适用于复杂起伏地表条件的叠前保幅高斯束偏移方法.在传统水平地表高斯束偏移的基础上,本文根据中心射线附近有效邻域内高斯束表征的近似波场,导出了起伏地表条件下具有相对振幅保持的高斯束偏移公式,并给出了一种精度更高的旁轴射线传播角度计算方法.同现有的高斯束偏移方法相比,本文方法不仅考虑了起伏地表对高斯束走时的线性影响,而且首次引入了由地表高程差异和近地表速度变化引起的二次时差校正项和振幅校正项,使得成像结果更加准确可靠.两个典型模型算例验证了本文方法的正确性和有效性.     

Common‐shot Fresnel beam migration based on wave‐field approximation in effective vicinity under complex topographic conditions

Gaussian beam migration is a versatile imaging method for geologically complex land areas, which overcomes the limitation of Kirchhoff migration in imaging multiple arrivals and has no steep‐dip limits of one‐way wave‐equation migration. However, its imaging accuracy depends on the geometry of Gaussian beam that is determined by the initial parameter of dynamic ray tracing. As a result, its applications in exploration areas with strong variations in topography and near‐surface velocity are limited. Combined with the concept of Fresnel zone and the theory of wave‐field approximation in effective vicinity, we present a more robust common‐shot Fresnel beam imaging method for complex topographic land areas in this paper. Compared with the conventional Gaussian beam migration for irregular topography, our method improves the beam geometry by limiting its effective half‐width with Fresnel zone radius. Moreover, through a quadratic travel‐time correction and an amplitude correction that is based on the wave‐field approximation in effective vicinity, it gives an accurate method for plane‐wave decomposition at complex topography, which produces good imaging results in both shallow and deep zones. Trials of two typical models and its application in field data demonstrated the validity and robustness of our method.     

Incorporating velocity shear into the magneto-Boussinesq approximation

Motivated by consideration of the solar tachocline, we derive, via an asymptotic procedure, a new set of equations incorporating velocity shear and magnetic buoyancy into the Boussinesq approximation. We demonstrate, by increasing the magnetic field scale height, how these equations are linked to the magneto-Boussinesq equations of Spiegel and Weiss (Magnetic buoyancy and the Boussinesq approximation. Geophys. Astrophys. Fluid Dyn. 1982, 22, 219–234).     

Seismic inversion with generalized Radon transform based on local second-order approximation of scattered field in acoustic media

Sound velocity inversion problem based on scattering theory is formulated in terms of a nonlinear integral equation associated with scattered field. Because of its nonlinearity, in practice, linearization algorisms (Born/single scattering approximation) are widely used to obtain an approximate inversion solution. However, the linearized strategy is not congruent with seismic wave propagation mechanics in strong perturbation (heterogeneous) medium. In order to partially dispense with the weak perturbation assumption of the Born approximation, we present a new approach from the following two steps: firstly, to handle the forward scattering by taking into account the second-order Born approximation, which is related to generalized Radon transform (GRT) about quadratic scattering potential; then to derive a nonlinear quadratic inversion formula by resorting to inverse GRT. In our formulation, there is a significant quadratic term regarding scattering potential, and it can provide an amplitude correction for inversion results beyond standard linear inversion. The numerical experiments demonstrate that the linear single scattering inversion is only good in amplitude for relative velocity perturbation ( \( \delta_{c}/c_{0} \) ) of background media up to 10 %, and its inversion errors are unacceptable for the perturbation beyond 10 %. In contrast, the quadratic inversion can give more accurate amplitude-preserved recovery for the perturbation up to 40 %. Our inversion scheme is able to manage double scattering effects by estimating a transmission factor from an integral over a small area, and therefore, only a small portion of computational time is added to the original linear migration/inversion process.     

Time-domain identification of dynamic properties of layered soil by using extended Kalman filter and recorded seismic data

A novel time-domain identification technique is developed for the seismic response analysis of soil-structure interaction. A two-degree-of-freedom (2DOF) model with eight lumped parameters is adopted to model the frequency-dependent behavior of soils. For layered soil, the equivalent eight parameters of the 2DOF model arc identified by the extended Kalman filter (EKF) method using recorded seismic data. The polynomial approximations for derivation of state estimators are applied in the EKF procedure. A realistic identification example is given for the layered-soil of a building site in Anchorage, Alaska in the United States. Results of the example demonstrate the feasibility and practicality of the proposed identification technique. The 2DOF soil model and the identification technique can be used for nonlinear response analysis of soil-structure interaction in the time-domain for layered of complex soil conditions. The identified parameters can be stored in a database for use in other similar soil conditions. If a universal database that covers information related to most soil conditions is developed in the future, engineers could conveniently perform time history analyses of soil-structural interaction.     

基于卫星遥感数据的全球表层流场反演重构

利用2000~2008年的卫星高度计资料和QuikSCAT风场资料,反演了全球的海表的地转流和Ekman流,将两者合成后生成了0．5°×0．5°的逐周全球表层流产品。在计算Ekman流的时候,引入了权重函数,改进了Lagerloef方法中Ekman流在25°S和25°N上的不连续问题。分析表明：卫星资料反演的流产品能够反映出海表流场的特征,将其分别于TAO观测和SGUD流产品进行定量化的比较显示,所得流产品具有较高的反演精度和可信度,说明改进的方法是有效的。     

深部咸水层CO2地质储存工程场地选址技术方法

CO₂地质储存作为环保型工程项目,其合理的工程场址是实现长期、安全封存CO₂的首要前提。我国CO₂地质储存工作刚刚起步,尚未形成成熟的选址技术方法体系。CO₂地质储存工程场地选址应遵循目标储层有效储存量大、安全、经济、符合一般建设项目环境保护选址条件、不受外部不良地质因素影响的原则,选址技术宜采用多尺度目标逼近法,选址程序包括规划选址和工程选址两大阶段。规划选址包括国家级、盆地级和目标区级潜力评价3个阶段;工程选址旨在通过目标靶区确定、综合地质调查、钻探及灌注试验和选定场地多因子排序综合评价,最终选出良好的工程场址。深部咸水层CO₂地质储存工程场地多尺度目标逼近选址技术方法对我国批量开展CO₂地质储存工程场地选址具有一定的指导意义。     

一种基于模糊正域类的模糊属性约简方法

为了实现模糊信息系统属性约简问题,探讨一种基于正域类的模糊属性约简方法。利用模糊二元相似关系建立模糊相似类的方法,进而定义决策属性关于条件属性的模糊正域类,得到正域类属性的重要度。结果表明:一方面,实现模糊信息系统属性约简;另一方面,通过取不同的相似精度,得到不同属性约简集。     

Role of Ekman transport versus Ekman pumping in driving summer upwelling in the South China Sea

Relative roles of Ekman transport and Ekman pumping in driving summer upwelling in the South China Sea (SCS) are examined using QuikSCAT scatterometer wind data. The major upwelling regions in the SCS are the coastal regions east and southeast of Vietnam (UESEV), east and southeast of Hainan Island (UESEH), and southeast of Guangdong province (USEG). It is shown that the Ekman transport due to alongshore winds and Ekman pumping due to offshore wind stress curl play different roles in the three upwelling systems. In UESEV, Ekman pumping and Ekman transport are equally important in generating upwelling. The Ekman transport increases linearly from 0.49 Sv in May to 1.23 Sv in August, while the Ekman pumping increases from 0.36 to 1.22 Sv during the same period. In UESEH, the mean estimates of Ekman transport and Ekman pumping are 0.14 and 0.07 Sv, respectively, indicating that 33% of the total wind-driven upwelling is due to Ekman pumping. In USEG, the mean Ekman transport is 0.041 Sv with the peak occurring in July, while Ekman pumping is much smaller (0.003 on average), indicating that the upwelling in this area is primarily driven by Ekman transport. In the summers of 2003 and 2007 following El Niño-Southern Oscillation (ENSO) events, both Ekman transport and Ekman pumping decrease in UESEV due to the abnormally weak southwest monsoon. During the same events, however, Ekman transport is slightly enhanced and Ekman pumping is weakened in UESEH and USEG.