小尺度网格层析速度建模方法在低幅构造区时深转换中的应用

渤海M油田区浅层多发育低幅构造,低幅构造区的地震数据对常规层析成像不够敏感,导致常规速度反演方法很难保证低幅构造速度建模的准确性,从而给时深转换工作带来较大误差.本文提出一种小尺度网格层析速度建模方法,该方法基于低幅构造的薄层信息,在有薄层信息的地层加密网格分布,从而实现根据地质条件划分网格的策略,进而通过解析矩阵方程实现速度模型的迭代更新.通过实际数据应用证明,对于渤海M油田利用本文提出的高精度网格层析速度建模方法是有效可行的,该方法能够有效提高低幅构造区速度建模的精度,利用该方法得到的速度模型能够大幅提高时深转换的精度,从而指导构造解释及成图工作.     

Impact of user parameters in inversion velocity analysis

Migration velocity analysis is a method devoted to the evaluation of both reflectivity and background velocity models, associated with the high and low wavenumber components of the model, respectively. Inversion velocity analysis is one of its improved versions, leading to more stable background velocity updates. Still, the impact of the user parameters should be understood for an optimal update of the background velocity. We show that a sign reversal of the background velocity gradient could occur when the selected surface offset range or the space lag range is too small. We derive the theoretical limits and check their consistency through simulations in a simple model with a single interface. These guidelines determine the necessary ranges of surface offsets and space lags for a proper update of the background velocity model. We discuss their applicability on the Marmousi model. Artefacts in the retrieved background velocity model are observed when the guidelines are not satisfied.     

The effectiveness of a pseudo-inverse extended Born operator to handle lateral heterogeneity for imaging and velocity analysis applications

Wave equation–based migration velocity analysis techniques aim to construct a kinematically accurate velocity model for imaging or as an initial model for full waveform inversion applications. The most popular wave equation–based migration velocity analysis method is differential semblance optimization, where the velocity model is iteratively updated by minimizing the unfocused energy in an extended image volume. However, differential semblance optimization suffers from artefacts, courtesy of the adjoint operator used in imaging, leading to poor convergence. Recent findings show that true amplitude imaging plays a significant role in enhancing the differential semblance optimization's gradient and reducing the artefacts. Here, we focus on a pseudo-inverse operator to the horizontally extended Born as a true amplitude imaging operator. For laterally inhomogeneous models, the operator required a derivative with respect to a vertical shift. Extending the image vertically to evaluate such a derivative is costly and impractical. The inverse operator can be simplified in laterally homogeneous models. We derive an extension of the approach to apply the full inverse formula and evaluate the derivative efficiently. We simplified the implementation by applying the derivative to the imaging condition and utilize the relationship between the source and receiver wavefields and the vertical shift. Specifically, we verify the effectiveness of the approach using the Marmousi model and show that the term required for the lateral inhomogeneity treatment has a relatively small impact on the results for many cases. We then apply the operator in differential semblance optimization and invert for an accurate macro-velocity model, which can serve as an initial velocity model for full waveform inversion.     

Attenuation Characteristics of Earthquake Ground Motion for Large Volume Airgun

In order to further deepen the understanding of seismic wave propagation characteristics induced by the large volume airgun source, experimental data from multiple fixed excitation points in Fujian Province were used to obtain the equivalent single excitation high signal-to-noise ratio velocity and displacement records through linear stacking and simulation techniques. Then the peak displacements of different epicentral distances were used to calculate the equivalent magnitude of the airgun source excitation at different fixed excitation points so as to establish the attenuation relationship between equivalent magnitude,epicenter distance and velocity peak. Our results show that:① Within 270 km of epicentral distance,for the large-volume airgun's single shot,the peak velocity range is about 700-4 nm/s,and the peak displacement range is about 200. 0-0. 2 nm;② The equivalent magnitude of the P-wave from the airgun source with a total capacity of 8,000 in 3 is 0. 181-0. 760,and the equivalent magnitude of the S-wave is 0. 294-0. 832. By contrast,the equivalent magnitude of the P-wave from the airgun source with a total capacity of 12,000 in 3 is 0. 533-0. 896,and the equivalent magnitude of the S-wave is 0. 611-0. 946. The S-wave energy is greater than the P-wave energy, and the excitation efficiency varies greatly with different excitation environment;③ The peak velocity increases with the equivalent magnitude,and decreases with the epicentral distance. The vertical component of the P-wave peak velocity is the largest among those three components,while the S-wave has the smallest vertical component and similar horizontal components. Hence,our research can provide an important basis for the quantitative judgment of the seismic wave propagation distance using the airgun and the design of the observation system in deep exploration or monitoring with airgun.     

晋冀蒙交界地区最小一维速度模型联合反演

文章使用内蒙古地震台网记录到的晋冀蒙交界地区517个定位地震,在参考原始编目报告的基础上,用和达法拟合走时曲线删除错误数据,使用Kissling方法联合反演得到晋冀蒙交界地区最小一维P波速度模型和震源位置。结果显示,研究区P波速度随着深度增加而递增,无明显间断面,联合反演得到的震中位置较原始报告分布更集中,簇状分布更明显。     

利用断裂带首波分析甘孜—玉树断裂带拉张盆地结构特征

对断裂带及附近区域地层精细结构的描述是理解地震孕育和发生过程的基础.不同板块或块体边界在构造或区域应力作用下,常常会形成速度间断面和低速的断裂带,间断面和低速带的存在会产生特殊的断裂带地震波,比如断裂带首波和围陷波,并会影响地震的物理过程、破裂速度和破裂方向等.本文主要利用2010年4月14日M_W6.9玉树地震后布设的流动地震台站,对沿着甘孜—玉树断裂带传播的断裂带首波进行识别和分析.分析结果表明,在甘孜—玉树断裂带的不同区域均观测到了断裂带首波,在西段的结隆拉张盆地附近的3个台站沿断裂带界面的平均纵波速度差异值在5%～8%,而其他区域的平均速度差异为1%～3%.通过台站位置分布和断裂带首波特征关系,初步分析了断裂带拉张盆地的结构,结果显示结隆盆地的长度为～40km,宽度为5.35～5.97km,深度不超过5km,在地表浅层形成了一个低速区,分别与巴颜喀拉块体(NE)和羌塘块体(SW)产生了两个物性差异界面,但没有延伸到主震和余震区震源深度.我们的结果表明结合密集台阵资料,通过断裂带首波特征分析可以为断裂带精细结构及几何特征提供一种新的技术方法和途径.     

鄂尔多斯西缘北段的地壳结构和块体间变形关系

本文对喜马拉雅计划二期部分台站的远震波形数据进行接收函数提取,利用接收函数共转换点叠加方法获得阿拉善地块、鄂尔多斯地块以及银川—河套盆地下方0~80 km深度的速度间断面结构.结果表明:鄂尔多斯地块成层性好,地壳厚度为38~42 km,康拉德界面为18~22 km,阿拉善地区的Moho面深度为38~45 km.河套盆地地壳厚度约52 km,银川断陷盆地和贺兰山下方的Moho面最深为~55 km.鄂尔多斯西缘构造边界下方Moho面变化明显,且黄河断裂为深大断裂直接切割莫霍界面.根据本文的间断面成像结果我们进一步确定阿拉善地块与鄂尔多斯地块分属不同的大地构造单元.与此同时,我们推测贺兰山以西70~80 km范围内和鄂尔多斯地块西缘北段存在地壳增厚变形的可能.     

Input-output Balances of Nutrients and Plankton in a Flooded Area of the Lower Odra

The retention of particulate matter of the Odra River in flooded areas was estimated by determining suspended particulate matter (SPM) elimination and particle-bound nutrient retention in a polder area of the Lower Odra Valley national park. Water and suspended matter samples collected before, during, and after the 1997 summer flood at the inlet and the outlet of the investigated polder (Polder A/B close to Schwedt) offer the opportunity to balance the matter retention inside the floodplain. The maximum level of retained SPM (more than 80%) was calculated for the record flood of summer 1997, while in ordinary winter floods retention differs between 33% and 70%. Basic properties of the flowing particles like settling velocity, density, loss on ignition etc. change on their way through the polder area. In the investigated Polder A/B a retention of at least 50% of particle-bound phosphorus which was independent of the incoming suspended matter load from the Odra River was observed. The retention of particle-bound nitrogen and carbon varied seasonally to different extents. The presented balance demonstrates that particle-bound nutrient retention in polders is an important factor in the self-purification process of the river system. Some additional effort was done to study changes of plankton composition: during a joint field experiment in spring 1998, the authors determined biological parameters (abundance and biomass of phytoplankton and zooplankton) of water samples at polder inflow and outflow points.     

Experimental study on mechanism of sea-dike failure due to wave overtopping

This work, which was largely a fruit of China's national marine hazard mitigation service, explicitly reveals the major mechanism of sea-dike failure during wave overtopping. A large group of wave-flume experiments were conducted for sea dikes with varying geometric characteristics and pavement types. The erosion and slide of the landward slope due to the combined effect of normal hit and great shear from overtopping flows was identified the major trigger of the destabilization of sea dikes. Once the intermittent hydrodynamic load and swash caused any deformation (bump or dent) of the pavement layer, pavement fractions (slabs or rubble) on the slope started to be initiated and removed by the water. The erosion of the landward slope was then gradually aggravated followed by entire failure within a couple of minutes. Hence, the competent velocity would be helpful evaluate the failure risk if as well accounted in standards or criteria. However, the dike top was measured experiencing the largest hydrodynamic pressure with a certain cap while the force on the wall increased rapidly as the overtopping intensity approached the dike-failure threshold. The faster increase of the force on the wall than on the landward slope yielded the sequencing of loads reaching hypothetic limits before failure as: dike top – top-mounted wall – landward slope. Therefore, beside the slide failure, the fatigue damage due to the instantaneous hydrodynamic impact might be another mechanism of the dike failure, which did not appear in the experiment but should be kept in mind. Instead of the widely adopted tolerable overtopping rate, a 0.117–0.424 m³/(m s) range of overtopping discharge and a 10 m/s overtopping velocity for the failure risk of typical sea dikes along China's coastlines were suggested, which enables the possible failure risk prediction through empirical calculations. The failure overtopping rate was identified strongly dependent on the pavement material, the landward slope and the dike-mounted wall but showed little variation with the width of the dike top. The flat concrete pavement and gentle landward slopes are suggested for the dike design and construction. For given configurations and hydrodynamic conditions in the experiment, the dike without the wall experienced less overtopping volume than those with the 1-m top-mounted wall. Meanwhile, the remove of the wall increased the failure overtopping rate, which means a certain increase of the failure criterion. Thus, care must be taken to conclude that the dike-mounted wall seems not an entirely appropriate reinforcement for the stability and safety of coastal protections. This should be further checked and discussed by researchers and engineers in the future.     

The validity of generic trends on multiple scales in rock-physical and rock-mechanical properties of the Whitby Mudstone,United Kingdom

Finding generic trends in mechanical and physical rock properties will help to make predictions of the rock-mechanical behaviour of shales. Understanding the rock-mechanical behaviour of shales is important for the successful development of unconventional hydrocarbon reservoirs.This paper presents the effect of heterogeneities in mineralogy and petrophysical properties on the validity of generic trends on multiple scales in rock-mechanical and rock-physical properties of the Whitby Mudstone. Rock-mechanical laboratory experiments have been performed on Whitby Mudstone samples from multiple outcrops within five kilometres laterally in order to investigate the heterogeneity and possible trends on an outcrop scale. Unconfined compression tests and acoustic measurements have been conducted to obtain the rock-mechanical properties, including rock strength, Young's modulus, Poisson's ratio, and velocity anisotropy. The rock-physical properties, including mineralogy, porosity, and matrix density, were measured using X-ray fluorescence and helium pycnometry. Various methodologies have been applied to the resultant data in order to derive different brittleness indices.Significant heterogeneity in rock-mechanical and rock-physical properties is present on an outcrop scale. There is no obvious correlation between mineral content and rock-mechanical properties on an outcrop scale in the Whitby Mudstone. Comparison with shales from different basins show, however, correlations between composition and elastic properties. The presence of significant heterogeneities on an outcrop scale and between shales from different basins make it difficult to find generic trends in rock-physical and rock-mechanical properties.