泡沫覆盖的随机粗糙海面的双站散射和热辐射

从双尺度随机粗糙面的散射理论和矢量辐射传输理论出发，将具有白帽泡沫覆盖的风驱海面模拟成一层离散球形粒子下垫双尺度随机粗糙界面的复合模型．利用Ｐｉｅｒｓｏｎ的半经验海洋波高谱和修正的Ｃｏｘ和Ｍｕｎｋ的坡度分布函数，推导了强风驱使下白帽泡沫覆盖的随机起伏海面的双站散射系数和热发射率，得到了极化的双站散射和热发射与海面各物理参数之间定量的函数关系．数值计算的结果与已有的海洋遥感实验数据进行了比较和讨论．     

海浪破碎对海面微波后向散射系数的影响

在海面波破碎区，假定随机粗糙的海面被两层不同散射性质的球形粒子覆盖，上层是由海面飞沫组成的球形Ｒａｙｌｅｉｇｈ粒子层，下层是由气泡组成的球形Ｍｉｅ粒子层，根据大量观测和实验的统计结果，确定了两层散射粒子的平均半径和单位体积的粒子数；并用迭代法求解随机粗糙海面上的两层粒子的矢量辐射传输方程．根据风波破碎时白冠覆盖率和海面飞沫、气泡的统计模型，给出了一个新的风驱海面的微波散射理论模型．利用这一模型定量地给出了海面后向散射系数和风速、视角、极化和白冠覆盖率等特征参量之间的函数关系，并研究了覆盖在波破碎区的散射粒子对海面后向散射系数影响程度．     

植被电磁散射的半空间模型研究

本文给出了地面植被电磁散射的半空间模型研究方法.在以往的相关文献中，均采用自由空间格林函数求解单个叶片散射体的散射场，本文利用半空间并矢格林函数求出了单个散射体的散射场，然后运用Monte Carlo方法模拟生成地面植被层，得到了在半空间下植被层的电磁散射特性，计算出了单、双站雷达散射截面，并与自由空间下的散射场做了相互对比.结果表明，在半空间格林函数下得到的散射场更为适用于描述地面植被的电磁散射特性.     

二维介质粗糙面下方三维金属目标复合电磁散射的快速正演算法

研究了二维(2-D)介质粗糙面下方三维(3-D)金属目标的复合电磁散射问题.将表面积分方程(PMCHW)方程应用到介质粗糙面表面,电场积分方程(EFIE)应用于金属目标表面.基于矩量法,使用三角分域基函数(RWG)和伽略金法将表面积分方程离散为矩阵方程,并采用稳定的双共轭梯度迭代(BICGSTAB)算法对矩阵方程进行求解.针对矩量法(MOM)的高存储量和迭代过程中存在的矩阵向量积耗时的瓶颈,采用基于秩的多层矩阵分解法(MLUV),对矩阵元素进行压缩存储,以节省对计算机内存的需求,并加速迭代过程中的矩阵向量积运算.计算了高斯粗糙面下方球体的双站雷达散射截面积(RCS),并与最陡下降快速多级子算法(SDFMM)结果比较以验证该数值方法的正确性.最后分析了不同粗糙度、目标尺寸和目标位置对双站RCS的影响.     

时变海面与其上方导体平板的复合电磁散射研究

应用当前的数值方法求解动态海面与目标的复合电磁散射,由于海面的变化,在不同时刻需要对阻抗矩阵各元素及海面表面电流重新求解,因而要耗费大量内存和运算量.为了克服这一问题,本文应用物理光学(PO)近似求解了导体海面表面电流及导体平板的一阶散射场,应用基尔霍夫近似给出了海面的后向散射场,同时借助互易性定理降低了求解平板和海面之间二次耦合散射场的难度,讨论了平板尺寸、风速等对后向复合散射场的影响.另外,本文还推导出了耦合散射场多普勒谱频移的理论公式,详细分析了复合后向散射场的Doppler 频谱特性.     

利用海面全极化微波辐射获取毛细重力波谱的研究

通常利用双尺度模型计算海面微波辐射的原则是, 先计算大波斜率密度函数对电磁扰动权函数的影响, 再引入海面毛细重力波谱的效应. 电磁扰动权函数与海面毛细重力波谱无关, 而且海浪长波谱和毛细重力波谱可以根据分界波数分开计算; 依据这一点, 推导并获得了一个简单的计算海面微波辐射的新模型. 依据该模型, 海面微波辐射可以表示为电磁扰动权函数对海面曲率谱的积分, 其优点是电磁扰动权函数和海面曲率谱对亮温的作用可以进行分解研究. 基于该模型, 进一步讨论了利用多频全极化辐射计反演海面毛细重力波谱的可行性.     

利用全极化微波辐射计资料反演台风境内海面风场

作为一种新兴的被动遥感技术,全极化微波辐射计不仅可以提供海面风速产品,还可以提供海面风向产品.以往利用全极化微波辐射计观测亮温进行海面风场反演仅在晴空条件下进行,本文通过对观测亮温结合台风区域海面风场的分布特征进行分析,验证了全极化微波辐射计具有在台风等恶劣天气条件下进行海面风场观测的能力.基于敏感性分析实验,确定使用6.8 GHz和10.7 GHz等低频通道组合可进行台风区域内海面风场反演.其中,海面风速反演使用基于统计的多元线性回归算法,同时对海面温度、大气水汽含量、云中液态水含量及降水强度等物理量进行反演计算,为海面风向反演做准备.海面风向反演使用物理统计法进行,借鉴散射计风向反演使用的最大似然估计法.通过在全极化辐射传输前向模型中加入降水对大气透过率的影响、设计第三和第四Stokes通道亮温环境影响修正函数,在实现台风区域内海面风向反演的同时减小了反演误差.通过对“云娜”台风境内海面风场进行数值计算,验证了本文反演算法的可行性,并对反演误差的空间分布特征进行了分析.将2004年各台风过程的海面风场反演结果与散射计风场产品进行对比,海面风速和海面风向反演的均方根误差分别为1.64 m·s^-1和18.02°.     

非高斯粗糙面电磁散射特性的谱估计方法

计算粗糙面电磁散射特性的方法有微扰法、基尔霍夫近似法以及双尺度法,这些方法都有各自的适用范围.根据近代谱估计的原理,基于粗糙面高度起伏的谱分布特征,提出了一种用近代谱估计方法计算粗糙面电磁散射的新方法.这种方法只要了解粗糙面的Fourier频域分布,就可计算粗糙面的电磁散射特征.最后对谱估计方法和传统的几种方法进行比较,表明这种方法在一定情况下能够对粗糙面的散射特征进行分析,而且与论述类似问题的文献的结果相当吻合.     

一阶Rytov近似有限频走时层析

传统的波动方程走时核函数(或走时Fréchet导数)多基于互相关时差测量方式及地震波场的一阶Born近似导出,其成立条件非常苛刻.然而,地震波走时与大尺度的速度结构具有良好的线性关系,对于小角度的前向散射波场,Rytov近似优于Born近似.因此,本文基于Rytov近似和互相关时差测量方式,导出了基于Rytov近似的有限频走时敏感度核函数的两种等价形式:频率积分和时间积分表达式.在此基础之上,本文提出了一种隐式矩阵向量乘方法,可以直接计算Hessian矩阵或者核函数与向量的乘积,而无需显式计算和存储核函数及Hessian矩阵.基于隐式矩阵向量乘方法,本文利用共轭梯度法求解法方程实现了一种高效的Gauss-Newton反演算法求解走时层析反问题.与传统的敏感度核函数反演方法相比,本文方法在每次迭代过程中,无需显式计算和存储核函数,极大降低了存储需求.与基于Born近似的伴随状态方法走时层析相比,本文方法具有准二阶的收敛速度,且适用范围更广.数值试验证明了本文方法的有效性.     

埋地目标体矢量电磁散射的一种快速正演算法

利用积分方程方法以及半空间并矢格林函数的快速算法对埋地目标体矢量电磁散射进行正演计算。首先,利用半空间电磁并矢格林函数建立起埋地目标体的体积分方程。然后通过将空间偏导转移至格林函数谱域积分的积分号之外,并采用离散复镜像方法来近似余下的零阶索末菲积分,进一步得到并矢格林函数各个分量的闭合形式。由于避免了对索末菲积分的繁琐数值计算,使得生成反应矩阵和计算散射场时由半空间并矢格林函数计算带来的瓶颈问题得到较好的克服,因而极大地提高了埋地目标体电磁响应正演计算效率,同时也能保证足够的精度。     

Series approximation of the Kirchhoff method for wavefield computation on a rough subsurface

ABSTRACT Computation of the wavefield due to reflection from an irregular surface is carried out for subsurfaces with large radii of curvature. The Kirchhoff approximation is proved to be sufficiently accurate provided that the acoustic wavelength is sufficiently small with respect to the asperities of the rough surface. For cases where the irregular surface does not fulfil this condition, a series solution is proposed. The first term of this series appears to be the result obtained by conventional Kirchhoff approximation. The series, initially developed in the space–wavenumber domain by Meecham, is transformed into the space–time domain, and the general expression for the series is obtained by calculation of the normal derivative of the field function. The series solution, restricted to the first two terms, is illustrated by application to three synthetic examples. Applications show that the series approximation obtained by the Kirchhoff method contributes significantly to the modelling of narrow, steep and deep structures and consequently it appears that the second term in the series cannot be ignored in the computation of the wavefields arising from a rough surface.     

Kirchhoff scattering series for sources and receivers located in an elastic medium: Another insight into the concept of virtual events

The convolution-type and correlation-type representation theorems are building blocks of wave-scattering theory whose usefulness expands in many seismological applications. For example, the Kirchhoff scattering series currently used for attenuating free-surface multiples has been derived from the convolution-type representation theorem. The recently introduced concept of virtual events, which allows us to put virtual sources and virtual receivers inside the subsurface based on the data collected at the sea surface, has been derived by a combined use of the convolution-type and correlation-type representation theorems. The formulation of inverse Kirchhoff scattering series and virtual events has been limited so far to the cases in which sources or receivers, or both, are located in the water. Unfortunately, this assumption is not valid, especially in the context of virtual events, in which both sources and receivers will often be located in a solid. We here redescribe the Kirchhoff scattering series and reformulate the concept of virtual events for the cases in which sources and receivers are in a solid. Moreover, we describe a new form of Kirchhoff series based on the correlation-type representation theorem and new formulae for computing virtual events which do not include the complex renormalization operation of the previous formulation.     

Internal multiple attenuation by Kirchhoff extrapolation

Surface‐related multiple elimination is the leading methodology for surface multiple removal. This data‐driven approach can be extended to interbed multiple prediction at the expense of a huge increase of the computational burden. This cost makes model‐driven methods still attractive, especially for the three dimensional case. In this paper we present a methodology that extends Kirchhoff wavefield extrapolation to interbed multiple prediction. In Kirchhoff wavefield extrapolation for surface multiple prediction a single round trip to an interpreted reflector is added to the recorded data. Here we show that interbed multiples generated between two interpreted reflectors can be predicted by applying the Kirchhoff wavefield extrapolation operator twice. In the first extrapolation step Kirchhoff wavefield extrapolation propagates the data backward in time to simulate a round trip to the shallower reflector. In the second extrapolation step Kirchhoff wavefield extrapolation propagates the data forward in time to simulate a round trip to the deeper reflector. In the Kirchhoff extrapolation kernel we use asymptotic Green's functions. The prediction of multiples via Kirchhoff wavefield extrapolation is possibly sped up by computing the required traveltimes via a shifted hyperbola approximation. The effectiveness of the method is demonstrated by results on both synthetic and field data sets.     

Exact time-domain solutions for acoustic diffraction by a half plane

We derive exact time-domain solutions for scattering of acoustic waves by a half plane by inverse Fourier transforming the frequency-domain integral solutions. The solutions consist of a direct term, a reflected term and two diffraction terms. The diffracting edge induces step function discontinuities in the direct and reflected, terms at two shadow boundries. At each boundary, the associated diffraction term reaches a maximum amplitude of half the geometrical optics term and has a signum function discontinuity so that the total field remains continuous. We evaluate solutions for practical point source configurations by numerically convolving the impulse diffraction responses with a wavelet. We solve the associated problems of convolution with a singular, truncated diffraction operator by analytically derived correction techniques. We produce a zero offset section and compare it to a Kirchhoff integral solution. Our exact diffraction hyperbola exhibits noticeable asymmetry, with higher amplitudes on the reflector side of the edge. Near the apex of the hyperbola the Kirchhoff solution approximates the exact diffraction term symmetric in amplitude about the reflection shadow boundary, but omits the other low amplitude term necessary to ensure continuity at the direct shadow boundary.     

Vector radiative transfer numerical model of coupled ocean-atmosphere system using matrix-operator method

A vector radiative transfer numerical model of the coupled ocean-atmosphere system is developed based on the matrix-operator method, which is named PCOART. Using the Fourier analysis, the vector radiative transfer equation (VRTE) is separated into a set of equations depending only on the observa-tion zenith angle. Using the Gaussian-Quadrature method, VRTE is finally transferred into the matrix equation solved by the adding-doubling method. According to the reflective and refractive properties of the ocean-atmosphere interface, the vector radiative transfer numerical model of the ocean and at-mosphere is coupled in PCOART. Compared with the exact Rayleigh scattering look-up tables of MODIS (Moderate-resolution Imaging Spectroradiometer), it is shown that PCOART is an exactly numerical model, and the processing methods of the multi-scattering and polarization are correct. Also, validated with the standard problems of the radiative transfer in water, it is shown that PCOART can be used to calculate the underwater radiative transfer problems. Therefore, PCOART is a useful tool for exactly calculating the vector radiative transfer of the coupled ocean-atmosphere system, which can be used to study the polarization properties of the radiance in the whole ocean-atmosphere system and the remote sensing of the atmosphere and ocean.     

Upward continuation of potential fields from a polyhedral surface1

An equivalent source procedure is derived for upward continuation of unevenly spaced gravity and magnetic data. The dipole layer is placed on a topographic relief approximated by a polyhedral surface, the stations being the vertices of the triangular faces. The dipoles have linear magnitudes, being directed along the normal vector over each triangle. The unknown values of the dipole magnitudes at each station are obtained by a suitable modification of the usual integral equation considering the discontinuity of the normal vector at each vertex of the dipole surface. Profile data processing is also studied. A numerical test outlines the accuracy and the limitations of the model for the case of a magnetic field significantly perturbed by a rough topographic relief.     

Reproducing kernel and Galerkin��s matrix for the exterior of an ellipsoid: Application in gravity field studies

In the introductory part of the paper the importance of the topic for gravity field studies is outlined. Some concepts and tools often used for the representation of the solution of the respective boundary-value problems are mentioned. Subsequently a weak formulation of Neumann??s problem is considered with emphasis on a particular choice of function basis generated by the reproducing kernel of the respective Hilbert space of functions. The paper then focuses on the construction of the reproducing kernel for the solution domain given by the exterior of an oblate ellipsoid of revolution. First its exact structure is derived by means of the apparatus of ellipsoidal harmonics. In this case the structure of the kernel, similarly as of the entries of Galerkin??s matrix, becomes rather complex. Therefore, an approximation of ellipsoidal harmonics (limit layer approach), based on an approximation version of Legendre??s ordinary differential equation, resulting from the method of separation of variables in solving Laplace??s equation, is used. The kernel thus obtained shows some similar features, which the reproducing kernel has in the spherical case, i.e. for the solution domain represented by the exterior of a sphere. A numerical implementation of the exact structure of the reproducing kernel is mentioned as a driving impulse of running investigations.     

Understanding the effect of elliptical polarization in surface nuclear magnetic resonance method

Currently, surface nuclear magnetic resonance (SNMR) method is the only geophysical method that detects groundwater directly. In this paper, we investigate the effect of elliptical polarization in the perpendicular excitation magnetic field. The effect of elliptical polarization is clearly visible in our ellipticity calculation and it can cause strong distortion to the excitation field in the presence of high subsurface conductivities. By examining the co-rotating and counter-rotating components of the field, we show that elliptical polarization affects transmitting and receiving processes differently and that a clear phase lag exists between transmitter loop and receiver loop. Finally, we derive the response function of coincident loops and calculate proton tip angles, the kernel function and SNMR response curves of a 1D aquifer model. Based on the simulations, we conclude that the elliptical polarization and phase lag can significantly affect SNMR response and it is essential to include elliptical polarization in SNMR modeling and data interpretation.