Stokes漂流近似公式对海洋表层流场估算的影响

本文采用波浪订正的Ekman模型,研究分析了三种Stokes漂流近似公式(单波公式、e指数公式、Phillips谱近似公式)对海洋表层流场估算的影响。海表总流场由海表面高度(SSH)数据计算的地转流和海浪模式WAVEWATCH Ⅲ输出结果计算的非地转流组成,并采用拉格朗日浮标观测数据对计算结果进行了验证。研究表明,随着Stokes漂流近似公式精度的提高,其计算的拉格朗日流速更接近于谱积分公式的计算结果,更贴近拉格朗日浮标观测数据。与谱积分公式计算的海表拉格朗日流速相比,单波公式的平均相对偏差为0.0834,e指数公式的平均相对偏差为0.0392,Phillips谱近似公式的平均相对偏差为0.0101,说明Phillips谱近似公式在不同风速下均能对谱积分公式有良好的近似效果。在低风速条件下,由Stokes漂流近似公式精度引起的海洋表层流场估算误差可以忽略不计,但随着风速增加,由近似公式精度引起的偏差逐渐变大,此时应该选择Phillips谱近似公式计算Stokes漂流,来减小误差。     

Parametrizing surface wave tomographic models with harmonic spherical splines

We present a mathematical framework and a new methodology for the parametrization of surface wave phase-speed models, based on traveltime data. Our method is neither purely local, like block-based approaches, nor is it purely global, like those based on spherical harmonic basis functions. Rather, it combines the well-known theory and practical utility of the spherical harmonics with the spatial localization properties of spline basis functions. We derive the theoretical foundations for the application of harmonic spherical splines to surface wave tomography and summarize the results of numerous numerical tests illustrating the performance of a practical inversion scheme based upon them. Our presentation is based on the notion of reproducing-kernel Hilbert spaces, which lends itself to the parametrization of fully 3-D tomographic earth models that include body waves as well.     

横向各向同性介质拟声波方程及其在逆时偏移中的应用

地下岩石的速度各向异性影响地震波的传播与成像.横向各向同性(TI)介质为最普遍的等效各向异性模型.引入TI介质拟声波方程可以避免复杂的弹性波方程求解以及各向异性介质波场分离,以满足对纵波成像的实际需要.本文从垂直横向各向同性(VTI)介质弹性波方程出发,推导出正应力表达的拟声波方程以及相应的纵波分量的表达式,进而分析从频散关系得到的拟声波方程的物理意义,而后将拟声波方程扩展到更一般的倾斜横向各向同性(TTI)介质中.波前快照与群速度平面的对比验证了拟声波方程可以很好地近似描述qP波的运动学特征.在此基础上,将拟声波方程应用在逆时偏移中并与其特例声波近似方程进行对比,讨论了计算效率、稳定性等实际问题.数值试验表明VTI介质情况下采用声波近似方程可以提高计算效率,而TTI介质qP-qSV波方程则在效率相当的情况下可以保证稳定性.SEG/HESS模型和逆冲模型逆时偏移试验验证了本文TI介质拟声波方程的实用性.     

Stability and identification for rational approximation of frequency response function of unbounded soil

Exact representation of unbounded soil contains the single output–single input relationship between force and displacement in the physical or transformed space. This relationship is a global convolution integral in the time domain. Rational approximation to its frequency response function (frequency‐domain convolution kernel) in the frequency domain, which is then realized into the time domain as a lumped‐parameter model or recursive formula, is an effective method to obtain the temporally local representation of unbounded soil. Stability and identification for the rational approximation are studied in this paper. A necessary and sufficient stability condition is presented based on the stability theory of linear system. A parameter identification method is further developed by directly solving a nonlinear least‐squares fitting problem using the hybrid genetic‐simplex optimization algorithm, in which the proposed stability condition as constraint is enforced by the penalty function method. The stability is thus guaranteed a priori. The infrequent and undesirable resonance phenomenon in stable system is also discussed. The proposed stability condition and identification method are verified by several dynamic soil–structure‐interaction examples. Copyright © 2009 John Wiley & Sons, Ltd.     

基于第二代Curvelet变换的地震资料随机噪声衰减

噪声衰减是地震资料处理中的关键问题之一。根据Curvelet变换对含有光滑边界的二维二阶连续可微函数所具有的稀疏表示性能,给出了Curvelet变换域地震资料随机噪声衰减的阈值方法;并给出了基于地震资料中随机噪声是独立同分布的高斯白噪假设条件下的阈值估计方法。通过合成数据和叠后实际数据算例,对该方法的有效性进行验证。结果表明,Curvelet变换不仅可以很好地衰减随机噪声,并且能较好地保持有效信号。      

Control strategy for variable damping element considering near‐future excitation influence

A control strategy is proposed for variable damping elements (VDEs) used together with auxiliary stiffness elements (ASEs) that compose a time‐varying non‐linear Maxwell (NMW) element, considering near‐future excitation influence. The strategy first composes a state equation for the structural dynamics and the mechanical balance in the NMW elements. Next, it establishes a cost function for estimating future responses by the weighted quadratic norms of the state vector, the controlled force and the VDEs' damping coefficients. Then, the Euler equations for the optimum values are introduced, and also approximated by the first‐order terms under the autoregressive (AR) model of excitation information. Thus, at each moment t_k, the strategy conducts the following steps: (1) identify the obtained seismic excitation information to an AR model, and convert it to a state equation; and (2) determine VDEs' damping coefficients under the initial conditions at t_k and the final state at t_k+L, using the first‐order approximation of the Euler equations. The control effects are examined by numerical experiments. Copyright © 2000 John Wiley & Sons, Ltd.     

弹道主动段全射向扰动引力快速逼近方法

目前快速逼近方法都必须在已知主动段弹道的条件下才能进行计算,这对于弹道导弹的机动发射(指随时的打击目标的改变)重力场保障带来困难。这里提出了一种能够实现对某种标准弹道的弹道导弹主动段全射向扰动引力快速逼近方法。试验表明:在3×10-5ms-2的精度要求下,在该实验区域只需要建立6组分频多项式模型就可以实现弹道导弹主动段全射向扰动引力的快速逼近。     

均匀半空间瞬变电磁场直接时域响应数值分析

近源时域电磁场具有信号强、探测深度大和精度高等优点,但传统瞬变电磁场理论中偶极子近似在近源区会引起较大误差,推导瞬变电磁场直接时域解析式是解决这一问题的关键.本文在点电荷微元假设下通过时域格林函数,采用分离变量等方法推导出了上半空间一次有源波动场和反射波的时域解析式和下半空间二次无源波动场的时域解析式,结合均匀半空间瞬变电磁场的边界条件给出了均匀半空间瞬变电磁场的直接时域解析式,进而利用第一型曲线积分,通过沿回线源叠加推导出圆回线源在瞬变电磁场中的直接时域解析式.然后在半空间表面上,与传统的电偶极源假设下的表达式作了比较.数值结果表明两者在远源区的计算结果相差甚微,而近源区则存在很大误差.本文利用真正点元（点电荷）严密推导给出的均匀半空间表面上瞬变电磁场的直接时域解析式适用于全场区探测,克服了偶极子假设下只适用远场区的不足,为瞬变电磁法的进一步发展和实际勘探提供了新的理论基础.     

声波介质一次散射波场高斯束Born正演

Born正演是一种常用的地震波场正演模拟方法,也是线性化地震反演的理论基础.在实际应用时,Born正演通常结合常规的地震射线方法进行实现.为了克服常规地震射线方法的弊端,并且保证地震波场的模拟精度和计算效率,本文提出了一种基于高斯束的一阶散射波场Born正演方法.该方法分为两个环节:首先,我们利用高斯束的走时和振幅信息将地下散射点处的反射率映射为地表束中心位置处的局部平面波;然后,我们利用逆倾斜叠加将局部平面波转化为接收点处的时空域散射波场.在具体的实施过程中,我们提出一种以wavelet-bank方式实现的局部平面波合成方法,同现有的算法相比,可以在保持计算精度的同时,大大减少计算时间;此外,我们还利用最速下降法优化了高斯束的迭代循环过程,进一步提高了Born正演的计算效率.两个模型的应用效果证明,本文所提出的高斯束Born正演方法可以精确、高效的实现声波介质一次散射波场的正演模拟,为三维大规模地震波场的正演问题提供了一种切实可行的实现方案.     

Design of raft–pile foundation using combined optimization and finite element approach

This paper describes the application of a structural optimization approach combined with the finite element method for the optimal design of a raft–pile foundation system. The analysis takes into account the non-linear behaviour of the soil medium and the piles. For the optimization process, the sensitivity analysis is carried out using the approximate semi-analytical method while the constraint approximation is obtained from the combination of extended Bi-point and Lagrangian polynomial approximation methods. The objective function of the problem is the cost of the foundation. The design variables are the raft thickness, cross-section, length and number of piles. The maximum displacement and differential displacement are selected as the constraints. The proposed method is shown to achieve an optimum design of raft–pile foundation efficiently and accurately. Copyright © 1999 John Wiley & Sons, Ltd.