P-wave traveltime and polarization tomography of VSP data

The presence of anisotropy requires that tomographic methods be generalized to account for anisotropy. This generalization allows geological structure to be correctly imaged and allows the anisotropic parameters to be estimated. Use of isotropic inversion for imaging anisotropic structures gives systematic trends in the traveltime and polarization residuals. However, due to the limited directional coverage, the traveltimes along may not be sufficient to study the anisotropic properties of the structure. Polarizations can provide independent information on the structure. Traveltime and polarization inversion are applied to synthetic examples simulating VSP experiments. Transverse isotropy and 1-D structure are assumed. Plots of traveltime and polarization residuals are an important tool to detect the anomalies due to the presence of anisotropy. For receivers located in anisotropic layers, polarization residuals display consistent anomalies of several degrees. The synthetic examples show that even the simple 1-D problem is difficult, when using direct arrivals only. Large a posteriori errors in anisotropic parameters are obtained by traveltime inversion in layers where available incidence angles are less than 45°. Resolution of the tomographic image of VSP data is greatly improved by a combination of traveltime and polarization information. In order to obtain accurate inversion results, the measurement error of polarization data should be kept to within a few degrees.     

三维物体上二阶波浪绕射的实时模拟

采用二阶时域理论对非线性波浪在任意三维物体周围的绕射问题进行了研究，对自由表面边界条件进行Taylor级数展开，应用摄动展开可以建立相应的边值问题，而且此边值问题的计算域不随时间变化，运用基于B－样条的边界元方法求解每一时刻的波浪场，二阶自由表面边界条件在时间上进行数值积分，在自由表面加了一个人工阻尼层以避免波浪的反射，速度势分解为已知的入射势和未知的散射势，初始条件采用二阶Stokes波浪场，通过加入物体表面边界条件，得到散射势在时间和空间上的发展，本文对圆柱所受规则波的二阶波浪力和波浪爬高进行了计算，数值结果表明此理论计算准确，效率高，数值稳定。     

Statistical Distribution of Depth-Integrated Local Horizontal Momentum for Second-Order Random Ocean Waves in Finite Water Depth

Based on the second-order random wave solutions of water wave equations in tinite water depth, statistical distributions of the depth-integrated local horizontal momentum components are derived by use of the characteristic function expansion method. The parameters involved in the distributions can be all detemained by the water depth and the wavenumber spectrum of ocean waves. As an illustrative example, a fully developed wind-generated sea is considered and the parameters are calculated for typical wind speeds and water depths by means of the Donelan and Pierson spectrum. The effects of nonlinearity and water depth on the distributions are also investigated.     

Spreading rates,rift propagation,and fracture zone offset histories during the past 5 my on the Mid-Atlantic Ridge; 25°–27°30′ S and 31°–34°30′ S

The southern Mid-Atlantic Ridge (MAR) is spreading at rates (34–38 mm yr⁻¹) that fall within a transitional range between those which characterize slow and intermediate spreading center morphology. To further our understanding of crustal accretion at these transitional spreading rates, we have carried out analysis of magnetic anomaly data from two detailed SeaBeam surveys of the MAR between 25°–27°30′S and 31°–34°30′S. Within these areas, the MAR is subdivided into 9 ridge segments bounded by large- and short-offset discontinuities of the ridge axis. From two-dimensional Fourier inversions of the magnetic anomaly data we establish the history of spreading within each ridge segment for the past 5 my and the evolution of the bounding ridge-axis discontinuities. We see evidence for the initiation and diminishment of small-offset discontinuities, and for the transition of rigid large-offset transform faults to less stable short-offset features. Individual ridge segments display independent spreading histories in terms of both the sense and amount of asymmetric spreading within each which have given rise to changes through time in the lengths of bounding ridge-axis discontinuities. Over the past 3–5 my, the short-offset discontinuities within the area have lengthened/shortened by approximately the same amount (∼ 10 km). During this same time period, larger-offset transform faults have remained comparatively constant in length. A shift in plate motion at anomaly 3 time may have given rise to change in the length of short-offset second-order discontinuities. However, the pattern of lengthening/shortening short-offset discontinuities we see is not simply related to the geometry of the plate boundary in these regions which precludes a simply relationship between plate motion changes and response at the plate boundary. We document a case of rapid (minimum 60 mm yr⁻¹) small-scale rift propagation, occurring between 2.5 and 1.8 my, associated with transition of the Moore transform fault to an oblique-trending ridge-axis discontinuity. Propagation across the Moore discontinuity and similar propagation within the 31°–34°30’S area may be associated with the reduced age contrast in lithosphere across second-order discontinuities. Total opening rates within our northern survey area decreased from anomaly 4′ to 2 time and rates within both areas have increased since the Jaramillo. Total opening rates measured for anomaly intervals differ along the plate boundary significantly, more than expected with changing distance to the pole of rotation. These differences imply a degree of short-term non-rigid plate behaviour which may be associated with ridge segments acting as independent spreading cells. Magnetic polarity transition widths from our inversion studies may be used to infer a zone of crustal accretion which is 3–6 km wide, within the inner floor of the rift valley. A systematic increase of transition width with age would be expected if deeper crustal sources dominate the magnetic signal in older crust but this is not observed. We present results from three-dimensional analysis of magnetic anomaly data which show magnetization highs located at the intersection of the MAR with both large- and short-offset discontinuities. Within the central anomaly the highs exceed 15 A m⁻¹ compared with a background of approximately 8–10 A m⁻¹ and they persist for at least 2.5 my. The highs may be caused by eruption of fractionated strongly magnetized basalts at ridge-axis discontinuities with both large and small offsets.     

波面位移非线性特征数值研究

基于Longuest-Higgins(1963)非线性海浪模型,在有限水深且存在均匀背景流的条件下,根据Song(2006)给出的波面位移二阶表达式,采用Combi海浪频谱计算了海表面定点波面位移时间序列和波面位移概率统计分布。分析了波面位移统计分布随风速、水深、反波龄和均匀背景流的变化特征和规律以及不同海况条件下二阶非线性项对波面位移统计分布的影响。结果表明:二阶非线性项使波面位移分布偏离正态分布,二阶非线性作用受风速、水深、反波龄和均匀背景流的影响。风速增大、水深降低、反波龄减小或者均匀背景流和风速传播方向相反均使波面位移二阶非线性项的作用加强,无因次波面位移概率密度分布的偏度和峰度随之增大,反之则二阶非线性项作用减弱。当均匀背景流和风速相同时,虽然使非线性项的作用减弱,但平均波面位移反而比静止水平面降低。当均匀背景流和风速相反时,虽然使非线性作用增强,但平均波面位移反而趋于静止水平面。得到如下结论:二阶非线性项对于波面位移有显著影响,数值模拟波面位移需要增加二阶非线性项。通过以上研究,提高了数值模拟波面位移的准确性,而波面位移是海浪最基本的特征量,从而增强了海浪模拟和预报的准确性,对海洋工程、海–气相互作用、上层海洋动力学等具有重要意义。     

VTI介质中的地层弹性参数反演

在VTI介质中，用非双曲线可以很好地近似P－P波中长排列情况下的旅行时；而对P－SV波来说，应用非双曲线却不能很好地近似P－SV波的施行时。对P－P波、P－SV波的旅行时进行三项泰勒展开时，它们的系数却包含了地震波的各向异性信息及其垂直速度信息，这为地震波的弹性参数反演提供了基础。实验分析表明，在中等排列长度下可以精确地得到了P－P波，P－SV波旅行时的三项泰勒展式系数，从而可以精确地进行地层弹性参数反演。     

二维地层结构和速度的走时反演

介绍了在倾斜均匀各向同性层状介质情况下,根据地震反射走时反演二维地层结构和地层速度的地震射线层析成像方法——走时反演。正演模型射线追踪是根据Fermat原理实现的,即求解满足该原理的非线性方程组得到射线与界面的交点,进而计算相应的走时。反演是先假设一初始模型,用最优化方法使射线追踪走时与观察走时的残差极小。最后计算了分辨矩阵和信息密度矩阵,以评价反演结果。对有噪情形也进行了反演。