利用SPAC法估算地壳S波速度结构

S波速度结构能够反映地球介质的物性差异,是地壳内低速区结构特征判别的重要依据.本文尝试利用空间自相关法(SPAC法)从地震台站微动信号的垂直分量中提取瑞利波相速度频散曲线,通过对频散曲线的反演获得地下介质的S波速度结构.以国家数字测震台网8个宽频带地震台站的实测微动数据为例,采用SPAC方法获得了首都圈地区北京附近约30 km 深度范围内的一维S波速度结构.结果表明,该区结晶基底埋深较浅约2 km;分别在5~8 km 和12~16 km 深处发育S波低速层;8 km 和 20 km 处是S波速度差异较大的速度分界面.这一结果与以往地震学及人工地震探测结果较为吻合,表明SPAC法估算地壳S波速度结构是可行、有效的.     

Relative dispersion in the Liguro-Provençal basin: From sub-mesoscale to mesoscale

Relative dispersion in the Liguro-Provençal basin (a subregion of the Mediterranean Sea) is investigated using clusters of surface drifters deployed during two Marine Rapid Environment Assessment (MREA) experiments covering different months in 2007 and 2008, respectively. The clusters have initial radii of less than 1 km, or an order of magnitude below a typical deformation radius (approximately 10-20 km). The data set consists of 45 original pairs and more than 50 total pairs (including chance ones) in the spatial range between 1 and 200 km. Relative dispersion is estimated using the mean square separation of particle pairs and the Finite Scale Lyapunov Exponents (FSLEs). The two metrics show broadly consistent results, indicating in particular a clear exponential behaviour with an e-folding time scale between 0.5 and 1 days, or Lyapunov exponent ?? in the range of 0.7-1 days⁻¹. The exponential phase extends for 4-7 days in time and between 1 and 10-20 km in separation space. To our knowledge, this is only the third time that an exponential regime is observed in the world ocean from drifter data. This result suggests that relative dispersion in the Liguro-Provençal basin is nonlocal, namely controlled mainly by mesoscale dynamics, and that the effects of the sub-mesoscale motions are negligible in comparison. NCOM model results are used to complement the data and to quantify errors arising from the sparse sampling in the observations.     

S-wave velocity images of the Dead Sea Basin provided by ambient seismic noise

Based on passive seismic interferometry applied to ambient seismic noise recordings between station pairs belonging to a small-scale array, we have obtained shear wave velocity images of the uppermost materials that make up the Dead Sea Basin. We extracted empirical Green’s functions from cross-correlations of long-term recordings of continuous data, and measured inter-station Rayleigh wave group velocities from the daily correlation functions for positive and negative correlation time lags in the 0.1–0.5 Hz bandwidth. A tomographic inversion of the travel times estimated for each frequency is performed, allowing the laterally varying 3-D surface wave velocity structure below the array to be retrieved. Subsequently, the velocity-frequency curves are inverted to obtain S-wave velocity images of the study area as horizontal depth sections and longitude- and latitude-depth sections. The results, which are consistent with other previous ones, provide clear images of the local seismic velocity structure of the basin. Low shear velocities are dominant at shallow depths above 3.5 km, but even so a spit of land with a depth that does not exceed 4 km is identified as a salt diapir separating the low velocities associated with sedimentary infill on both sides of the Lisan Peninsula. The lack of low speeds at the sampling depth of 11.5 km implies that there are no sediments and therefore that the basement is near 10–11 km depth, but gradually decreasing from south to north. The results also highlight the bowl-shaped basin with poorly consolidated sedimentary materials accumulated in the central part of the basin. The structure of the western margin of the basin evidences a certain asymmetry both whether it is compared to the eastern margin and it is observed in north–south direction. Infill materials down to ∼8 km depth are observed in the hollow of the basin, unlike what happens in the north and south where they are spread beyond the western Dead Sea shore.     

交叉偶极子声测井多频率反演地层三阶弹性常数和地应力研究

基于声弹性理论,本文详细论述了利用交叉偶极子多频率相速度反演地层三阶弹性常数和地应力的理论及信噪比对反演结果的影响.定性和定量分析表明在速度应力方程中,不含三阶弹性常数项与含三阶弹性常数项相比对速度的影响可以忽略,由此只利用一次测井资料同时反演地应力和三阶弹性常数是一个严重病态的问题或一个不确定性问题.针对不同的反演模型,详细考察了信噪比对各自简化模型和精确模型反演结果的影响.数值模拟结果表明,利用一次测井采用精确模型反演全部四个参数方法受信噪比影响严重.而利用不同流体压力下的测井资料进行反演时可以采用简化模型取代精确模型以便简化计算和提高反演速度,因为简化模型和精确模型有着几乎相同的精度;同时,分步反演的方法可以较好地反演出原位地应力及三阶弹性常数.     

多极子随钻声波测井波场模拟与采集模式分析

应用三维非均匀交错网格有限差分程序模拟了多极子随钻声波测井的波场,考察了硬地层中单极子随钻声波测井响应及软地层中声源频率对多极子随钻声波测井的影响;计算了接收阵列波形的时域相干谱及频散特征,结果表明随钻测井仪器居中时,在硬地层中可以直接利用单极子声系得到地层横波速度,但地层纵波受到钻铤波的极大干扰.而在软地层中可以利用单极子声系直接测得地层纵波速度,但地层横波速度则需要使用工作在较低频率下的四极子随钻声波仪来测得,此时测井频段又会落入钻井噪声频带而受到影响.为此本文提出数值模拟手段结合物理实验的方法得到钻铤波的影响,在实际测井中通过信号处理方法而不是用在仪器本体上刻槽的方式来更完全地消除钻铤的影响.为了避免钻井噪声的影响,本文提出六极子随钻声波仪更适合软地层的横波测量.本文用数值实验结果来说明这些方法的可行性.     

Uncertainty propagation of hydrodispersive transfer in an aquifer: an illustration of one-dimensional contaminant transport with slug injection

It is evident that the hydrodynamic dispersion coefficient and linear flow velocity dominate solute transport in aquifers. Both of them play important roles characterizing contaminant transport. However, by definition, the parameter of contaminant transport cannot be measured directly. For most problems of contaminant transport, a conceptual model for solute transport generally is established to fit the breakthrough curve obtained from field testing, and then suitable curve matching or the inverse solution of a theoretical model is used to determine the parameter. This study presents a one-dimensional solute transport problem for slug injection. Differential analysis is used to analyze uncertainty propagation, which is described by the variance and mean. The uncertainties of linear velocity and hydrodynamic dispersion coefficient are, respectively, characterized by the second-power and fourth-power of the length scale multiplied by a lumped relationship of variance and covariance of system parameters, i.e. the Peclet number and arrival time of maximum concentration. To validate the applicability for evaluating variance propagation in one-dimensional solute transport, two cases using field data are presented to demonstrate how parametric uncertainty can be caught depending on the manner of sampling.     

Multi-Source Emission Determination Using an Inverse-Dispersion Technique

Inverse-dispersion calculations can be used to infer atmospheric emission rates through a combination of downwind gas concentrations and dispersion model predictions. With multiple concentration sensors downwind of a compound source (whose component positions are known) it is possible to calculate the component emissions. With this in mind, a field experiment was conducted to examine the feasibility of such multi-source inferences, using four synthetic area sources and eight concentration sensors arranged in different configurations. Multi-source problems tend to be mathematically ill-conditioned, as expressed by the condition number κ. In our most successful configuration (average κ = 4.2) the total emissions from all sources were deduced to within 10% on average, while component emissions were deduced to within 50%. In our least successful configuration (average κ = 91) the total emissions were calculated to within only 50%, and component calculations were highly inaccurate. Our study indicates that the most accurate multi-source inferences will occur if each sensor is influenced by only a single source. A “progressive” layout is the next best: one sensor is positioned to “see” only one source, the next sensor is placed to see the first source and another, a third sensor is placed to see the previous two plus a third, and so on. When it is not possible to isolate any sources κ is large and the accuracy of a multi-source inference is doubtful.     

槽波层析成像方法在煤田勘探中的应用——以河南义马矿区为例

在地下煤田的开发中,工作面内的小构造、异常体、煤层厚度变化等是需要解决的关键问题,而槽波探测则为这些问题的解决提供了重要的物探方法.本文对河南义马矿区11061工作面进行槽波透射法测量,在巷道显示的煤层厚度变化为1.5~8 m,从理论频散曲线分析速度与厚度关系,确定125 Hz频率槽波主要用于观测厚度约为2~5 m的煤层厚度变化;有效提取了684个频散曲线,并分别拾取了125 Hz时槽波群速度与走时,采用走时层析成像方法获得工作面内煤层速度、厚度以及高应力区分布特征,回采验证了结果的正确性.     

混合ADI-FDTD亚网格技术在探地雷达频散媒质中的高效正演

提出混合ADI-FDTD亚网格技术开展频散介质GPR正演,即在物性参数变化剧烈局部区域采用细网格剖分ADI-FDTD计算,其他的区域采用粗网格剖分常规FDTD计算,ADI-FDTD突破了CFL条件的限制,可选取与粗网格一致的大时间步长,有效地提高了计算效率.本文首先基于Debye方程,推导了粗网格FDTD及细网格ADI-FDTD频散介质差分格式,着重对粗细两种网格结合的场值交换方式进行了深入探讨,给出了该算法的计算流程.然后以一个薄层模型为例,分别应用粗网格、细网格、混合ADI-FDTD亚网格算法对该模型进行正演,计算资源的占用及模拟精度说明了混合ADI-FDTD亚网格算法的优势.最后,建立频散介质与非频散介质的组合模型,应用3种方法对该模型进行正演,对比3种方法优劣,分析雷达剖面中非频散介质及频散介质中波形特征,有效地指导雷达资料的精确解释.     

A Large-Eddy Simulation and Lagrangian Stochastic Study of Heavy Particle Dispersion in the Convective Boundary Layer

Large-eddy simulation and Lagrangian stochastic dispersion models were used to study heavy particle dispersion in the convective boundary layer (CBL). The effects of various geostrophic winds, particle diameters, and subgrid-scale (SGS) turbulence were investigated. Results showed an obvious depression in the vertical dispersion of heavy particles in the CBL and major vertical stratification in the distribution of particle concentrations, relative to the passive dispersion. Stronger geostrophic winds tended to increase the dispersion of heavy particles in the lower CBL. The SGS turbulence, particularly near the surface, markedly influenced the dispersion of heavy particles in the CBL. For reference, simulations using passive particles were also conducted; these simulation results agreed well with results from previous convective tank experiments and numerical simulations.