Ekman边界层动力学的理论研究

大气边界层及其与自由大气之间的相互作用具有明显的非线性特征,而这些特征是经典Ekman理论所不能描述的,因此,发展中等复杂程度(介于完全模式与经典Ekman模型之间)的大气边界层动力学模式,简称中间模式,对人们从理论上认识大气边界层动力学过程的非线性特征具有重要意义。本文对目前最具代表性的几个中间边界层模型:地转动量近似边界层模型、Ekman动量近似边界层模型以及弱非线性边界层模型进行了总结和分析,阐述了Ekman层主要动力学特征。通过分析上述各模型的理论框架,揭示了各模型的物理意义及其在描述Ekman边界层基本动力特征上的优点和局限性,并指出尽管在细节定量描述上有差异,但各中间模型对Ekman层动力学特征的定性描述具有很好的一致性。对于这些Ekman边界层近似理论模型的进一步应用问题,主要回顾和总结了利用上述模型探讨地形边界层结构、大气锋生过程、低层锋面结构和环流以及边界层日变化、低空急流形成等动力学问题的研究,并对这些研究所揭示的Ekman层动力学特征及其对自由大气低层运动的影响进行了分析,结果表明,这些Ekman边界层近似模型可以较好地揭示大气边界层动力学特征,在大气边界层动力学及其与自由大气相互作用的研究上具有重要价值。另外,还对目前Ekman边界层理论研究中存在的问题进行了一些分析,提出了有待进一步研究的科学问题。     

Recursive evaluation of interaction forces and property matrices from unit-impulse response functions of unbounded medium based on balancing approximation

Starting from the unit-impulse response matrix of the unbounded medium, a discrete-time formulation permitting the recursive evaluation of the interaction forces and a continuous-time formulation yielding property matrices corresponding to a model with a finite number of degrees of freedom are discussed. This is achieved using the balancing approximation method which is easily automated, guarantees stability and leads to highly accurate results. © 1998 John Wiley & Sons Ltd.     

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

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

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.     

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

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

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

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