联合PP、PS波叠前AVA反演方法

从单一pp波数据中提取出来的弹性参数与岩性以及孔隙内流体之间的联系是不充分的,这就造成了在油藏的解释中包含有大量的不确定性.通过加入PS波数据信息,可以有效降低反演的多解性,提高反演的稳定性.笔者所研究的反演方法是基于精确的Zoeppritz方程,利用广义线性反演方法,通过迭代的思想来反演出纵波速度、横波速度和密度.该方法相对单一pp反演可以有效提高横波速度的估计精度.     

Test of the upper mantle low velocity layer in Siberia with surface waves

The existence of the upper mantle low velocity layer (LVL) below 100 km depth in cratonic areas is tested with surface waves dispersion curves. Given the ambient noise we find that a pronounced LVL (80 km thick and 2% velocity reduction or 40 km thick and 5% velocity reduction) can be distinguished from a constant velocity model by comparison of the fundamental mode group velocities, whereas a thin LVL (less than 40 km thick) with small velocity contrast (less than 2%) cannot be resolved. The fundamental modes of Love and Rayleigh waves have similar properties and, in general, the phase velocity differences are smaller than the standard error. Phase velocity alone cannot discriminate between the models, and the group velocity is in general more sensitive to the velocity structure than the phase velocity. The higher modes at short periods could potentially determine a LVL but in reality it is difficult to obtain sufficiently accurate measurements. We invert the synthetic dispersion curves by the non-linear Hedgehog inversion method. A pronounced LVL (more than 40 km thick and with a strong velocity contrast of about 5%) is detectable by the non-linear inversion but for a thin LVL with a strong velocity contrast it is not possible to resolve both velocity and thickness. In the inversions all solutions include a LVL for models with a pronounced LVL, whereas the solution space includes models with and without a LVL for models with a zero or positive gradient velocity–depth structure.We invert also real data with travel path across the Siberian craton with the Hedgehog method. Almost all solutions include a LVL in the depth range of 80–150 km with a velocity contrast up to 2% to the surrounding intervals. Hence, the LVL appears to be a common feature of the Siberian upper mantle, although a constant velocity at the same depth range cannot be totally excluded. Despite low resolution at large depth, a pronounced asthenospheric LVL below a depth of about 225 km is a constant characteristic of the set of solutions.     

基于多目标函数的黏弹性全波形反演

全波形反演是勘探地球物理领域兴起的核心技术之一,不但可以构建地下速度结构,也能够反演衰减参数(品质因子Q)模型,有助于识别地下介质类型和构造(如流体和煤层陷落柱),对煤和油气等自然资源的勘探和开发有重要意义。参数串扰是黏弹性全波形反演的关键难点,受速度误差影响,反演的Q模型会包含非常强的串扰噪声。针对该问题,提出了基于多目标函数的黏弹性全波形反演理论与方法,首先通过旅行时反演速度结构,再通过中心频率目标函数反演Q模型,最终使用波形差目标函数同时反演速度和Q模型。由于中心频率主要受衰减影响,因此,可有效减弱速度误差对Q反演的影响。最后,通过数值模拟验证了算法可有效地反演速度和Q模型。     

刚度缓变介质中瑞利波特性

基于薄层刚度矩阵方法,研究了瑞利波在3种典型刚度缓变介质中传播特性,并对瑞利波在缓变介质及分层介质中传播特性进行比较。结果显示,无论剖面刚度是缓变还是分层变化,对规则剖面,表面波场基阶模态占主导地位,远场表面波有效相速度与基阶模态相速度接近;对不规则剖面,表面波场中高阶模态影响较大,影响程度及频率范围与剖面刚度变化有关。在缓变介质中,模态相速度高频趋势值高于剖面最小剪切波速介质的瑞利波速,有效相速度趋势值也明显偏离表面介质瑞利波速,这些特性与分层介质不同。研究结果显示,当缓变介质用数量有限的虚构分层来模拟时,表面波频散测试数据与理论值匹配程度及剖面分析精度会受虚构层数量影响。     

Lower-crust S-wave velocity beneath western Yunnan Province from waveform inversion of dense seismic observations

We use teleseismic body waveforms to explore S-wave layered velocity structures beneath 30 portable digital seismic stations deployed around western Yunnan Province. Results show that the Moho depth in this region is ∼40 km and decreases in general from north to south, consistent with previous geophysical studies. Associated with this lateral variation of the Moho depth, the lower crust above the Moho discontinuity has a 15–25 km thick zone with an S-wave velocity lower than that of the upper crust. This lower velocity zone might be interpreted as a lower crust weak channel, which may mechanically partially decouple the upper-crust deformation from the underlying mantle. Thus, the inverted S-wave velocity structure could provide new evidence for the lateral flow of lower crust in the build-up of the south-eastern Tibetan plateau.     

The Sensitive Properties of Hydrate Reservoirs Based on Seismic Stereoscopic Detection Technology

Higher-precision determinations of hydrate reservoirs, hydrate saturation levels and storage estimations are important for guaranteeing the ability to continuously research, develop and utilize natural gas hydrate resources in China. With seismic stereoscopic detection technology, which fully combines the advantages of different seismic detection models, hydrate formation layers can be observed with multiangle, wide-azimuth, wide-band data with a high precision. This technique provides more reliable data for analyzing the distribution characteristics of gas hydrate reservoirs, establishing velocity models, and studying the hydrate-sensitive properties of petrophysical parameters;these data are of great significance for the exploration and development of natural gas hydrate resources. Based on a velocity model obtained from the analysis of horizontal streamer velocity data in the hydrate-bearing area of the Shenhu Sea, this paper uses three VCs(longitudinal spacing of 25 m) and four OBSs(transverse spacing of 200 m) to jointly detect seismic datasets consisting of wave points based on an inversion of traveltime imaging sections. Accordingly, by comparing the differences between the seismic phases in the original data and the forward-modeled seismic phases, multiangle coverage constraint corrections are applied to the initial velocity model, and the initial model is further optimized, thereby improving the imaging quality of the streamer data. Petrophysical elastic parameters are the physical parameters that are most directly and closely related to rock formations and reservoir physical properties. Based on the optimized velocity model, the rock elastic hydrate-sensitive parameters of the hydrate reservoirs in the study area are inverted, and the sensitivities of the petrophysical parameters to natural gas hydrates are investigated. According to an analysis of the inversion results obtained from these sensitive parameters, λρ, Vp and λμ are simultaneously controlled by the bulk modulus and shear modulus, while Vs and μρ are controlled only by the shear modulus, and the latter two parameters are less sensitive to hydrate-bearing layers. The bulk modulus is speculated to be more sensitive than the shear modulus to hydrates. In other words, estimating the specific gravity of the shear modulus among the combined parameters can affect the results from the combined elastic parameters regarding hydrate reservoirs.     

地表大回线源在任意层状介质中产生磁场的计算

在地面核磁共振数值模拟中,地下介质中任意一点处磁场快速计算非常重要。作者然对水以状介质,推导出利用递推方法求得任意一层中电磁波的透射和反射系数的简单方法,从而为精确、确定地计算位于地表的大回线源在地下任意一点处产生的磁场创造了有利条件。     

考虑龄期分层的固体废弃物填埋场边坡稳定分析

固体废弃物填埋场长期使用过程中,垃圾土的强度参数随着龄期的变化而变化。在进行固体废弃物填埋场边坡稳定分析时考虑按龄期分层,基于极限分析上限理论,提出了转动-平动组合和多块体平动两种破坏机构。通过水平分层土坡算例的对比分析,验证两种破坏机构的合理性和有效性。进而对填埋场算例按龄期分层分析,结果表明,极限分析上限解的安全系数和破坏机制与强度折减有限单元法非常接近,多块体平动机构通过随机搜索所得的最危险滑裂面优于转动-平动组合机构。最后分析西班牙Coll Cardús固废填埋场,此填埋场是按龄期分层的典型复杂填埋场,运用两种破坏机构所得的安全系数和最危险滑裂面与有限单元法相符。算例分析表明,在不考虑水位的情况下,填埋体堆填时间越长越有利于填埋场的稳定。当填埋体按低、中、高龄期分层形成时,填埋场的整体稳定安全系数介于中龄期和高龄期的均质边坡之间。     

基于方位相交的纵波AVA数据运用SVD反演HTI介质裂缝密度

裂缝密度是定量预测裂缝性储层的重要参数之一。Varela利用方位角上的反射系数,运用SVD方法反演HTI介质裂缝密度。笔者将不同方位角测线中两两相交测线的反射系数做差,对反射系数差值运用SVD方法反演HTI介质裂缝密度,各向异性参数在做差后的反射系数中具有更显著的作用。笔者首先针对2种方法的基本原理进行了阐述,然后进行数值计算,并对反射系数存在随机扰动情况下2种方法的反演结果进行对比,最后对笔者提出的改进方法进行稳定性分析。反演结果显示：噪声较小时,2种方法的反演结果都接近真实值,然而当增大噪声时,改进的方法比原方法反演结果更好,因此改进的方法比原方法更具有抗噪能力;随着相交测线夹角的增大,改进的方法受随机扰动影响的稳定性增强,在方位正交时受随机扰动的影响最小;纵波高信噪比AVA数据越多,反演的结果越稳定。     

Evaluating the structural model from ambient vibration recordings and shallow seismic refraction techniques. A case study of Zagazig City,Egypt

This study uses the seismic refraction and noise measurements to investigate the velocity structure of the subsurface and emphasize the advantage of ambient vibration over the conventional seismic refraction technique. Field measurements were carried out at nine sites in and around Zagazig city. Shallow seismic refraction data were interpreted using the delay time method to obtain the two-dimension ground model at each site. Ambient vibration arrays are used to infer the one-dimensional compressional and shear wave velocity profiles through two main steps. The first step is to derive the dispersion curve from the recorded signals using the frequency-wavenumber method. The second is to invert the dispersion curve to obtain the site velocity profiles. The results of the compressional wave velocities obtained from seismic refraction technique showed that the subsurface consists of a number of layers ranging from two to four layers and give a good agreement with the results of the seismic wave velocities obtained from the ambient vibration arrays. The ambient vibration arrays gave a deeper depth of penetration than the other method, providing more information on the subsurface structure without any disturbance to the environment. This work provides reliable estimates of the seismic velocity structures of both shallow and deep sedimentary layers within the area of interest.     

Estimation of maximum saturated depth in two-layered drainage blankets over the barrier in landfill cover system

For accurately calculating the maximum saturated depth over the landfill barrier and analyzing the variation of liquid head with inflow time, a new method is presented that considers the unsteady-state condition. Based on the principle of water balance and the extended Dupuit assumption, the solution of the new method can be calculated. Also, this method is capable of revealing the effect of the multilayered drainage media. One of the advantages of the method presented in this paper is that it can be used to draw a phreatic surface in layered drainage media under unsteady-state condition and estimate the location of the maximum liquid depth, for both homogeneous conditions and two-layered drainage conditions. Based on the developed analysis method, the effects of a parametric variation on the maximum saturated depth and the variation of saturated depth with time are presented. The calculated results show that with the increase of slope gradients of cover system, the saturated depth of cover system is easy to reach a stable value. Moreover, when the drainage layer consists of two layers with different thickness and hydraulic conductivity, the increase of hydraulic conductivity of the upper drainage layer will make it easier for the saturated depth of cover system to reach a stable and lower value.     

The Crustal Structure and Seismic Activity in North China

A layered crustal block model of North China has been constructed based on large amount of data from seismic sounding carried out in recent two decades. Some deep fault zones, such as the Zhangjiakou.Penglai and Tancheng-Lujiang fault zones, divide the upper crust of North China into three upper crustal terranes and nine bolcks. There are distinct differences in velocity and depth distributions, which reflects Cenozoic block faulting in North China in the process of formation of the deep structure. The upper crust shows the features of transition in isostatic adjustment. The existence of a low-velocity layer in the middle crust is characteristic of the crustal structure in North China. There seems to be an increase of rheology of the rocks in the lower crust and a persistence of stable regional stress field. The patterns of the Moho on two sides of the Yanshan-Taihang Mountains are different. The relief of the Moho around Beijing, Shijiazhuang and Guangrao where the deep faults join together shows a quadrantal distribution in some degree. The dynamic sources for seismic activity are the NE-SW horizontal compression and the diapirism of the upper mantle. The middle and upper crust, especially the layered block structure has the most significant effects on seismicity, and the occurrence of earthquakes is more closely related to them than to the Moho.     

Tomographic studies of the upper mantle beneath the Italian region*

Teleseismic P arrivals at seismological stations are inverted into a model of velocity perturbations down to a depth of about 470 km. Directionally independent average residuals, computed from steeply inciding waves, are transformed into a model of lithospheric thickness. Both models show a good correspondence with the main tectonic features of the Italian Peninsula. Positive velocity perturbations are observed beneath the Alps and in depths over 200 km also beneath the Po Basin. A high-velocity anomaly of the Tyrrhenian subduction is less pronounced, probably due to a directional dependence of P velocities in the mantle. Negative velocity perturbations indicate several low-velocity regions, e.g. beneath the Northern Apennines, the Sicily region and in the upper 100 km beneath the Po Basin. The amplitudes of velocity perturbations beneath the depth of 200 km are smaller on the average than those in the upper two layers. The whole region is characterized by large undulations of the lithosphere base which reaches depths from less than 60 km to more than 150 km. The most prominent lithospheric root beneath the Alps is a product of the collision between the European and the Adriatic plates while the lithospheric thickening beneath the Calabrian coast is likely to be connected with the eastern wing of the Tyrrhenian subduction. The dramatic changes of lithosphere thickness between the northern and the southern Apenninic arcs and northern Calabria as well as the thinnings at the western closure of the Po Basin, indicate important deep-seated boundaries of lithospheric blocks of autonomous geodynamic development.     

纵横波速度反演在气层识别中的应用

提出一种部分叠加资料的反演方法。部分叠加资料保留了原有的大量信息,在测井及构造层位的约束下,建立最符合地下实际情况的初始模型。通过非线性最优化理论,采用逐道外推技术,能同时反演出每一个点更准确的纵波和横波速度,进而可以求得纵横波速度比及泊松比,为储层预测和储集层内的流体识别提供丰富的弹性参数信息。     

Propagation anomalies in Northwestern Italy by inversion of teleseismic residuals

Travel-time residuals of teleseismic P waves were analysed in order to elucidate the crust–upper mantle structure in Northwestern Italy, and the Western Alpine Arc. Using digital data obtained from both fixed seismograph networks operating in NW Italy (notably Liguria–Piedmont) and temporary arrays with the aid of cross-correlation techniques reliable travel-time residuals were calculated which were then inverted to obtain models of propagation anomalies. The reliability of the inversion procedure was tested using synthetic data. The model thus obtained appears to be stable and shows strong lateral heterogeneities at a litho–asthenospheric level; in particular, it confirms the high velocity contrast caused by the ‘Ivrea Body’ in the shallower layers and the presence of Alpine ‘roots’ reaching down to at least 200 km. A statistical analysis performed on the propagation times of rays crossing the resulting four-layered model reveals rms below 0.1 s.     

宜兴抽水蓄能电站试验洞的反分析研究

结合江苏宜兴抽水蓄能电站地下厂房试验洞量测资料的分析,研究了成层各向异性岩体中洞室工程问题的反分析方法。采用分层均质模型进行反分析计算。计算区域的各层岩层均被分别视为均质各向同性体,初始地应力假设呈线性分布。有限元分析中,对倾斜成层的岩层、软弱夹层及断层节理等分别划分了单元,据以模拟成层岩体的各向异性特征。为使计算过程简化,对离洞室较远的岩层拟按地质资料对材料特性参数赋值,反分析计算的目标未知数为初始地应力和试验洞相邻围岩的弹性模量。结果表明,与工程实践符合较好。     

南北构造带天水、武都强震区地壳和上地幔顶部结构

利用两条相互垂直的高分辨地震折射/宽角反射剖面和相应的非纵观测的多个扇形剖面取得的人工地震资料, 研究天水和武都8级大震区的地壳和上地幔顶部结构和构造.二维剖面结果显示, 地壳沿垂向可分为上地壳和下地壳两大层.上地壳中部存在低速层, 层内介质速度比背景值低0.3~0.5km/s.莫霍面深度大约为46~48km.NE向的天水-武都剖面下地壳速度在横向上变化剧烈, NW向的成县-武山剖面, 在礼县以西, Moho面和C界面有被上涌物质改造过的迹象.三维速度成像显示, 在105°E附近, 从7至11km的深度范围内, 存在一条近NS向的断裂带, 在该带的两侧速度结构有明显的差异, 西侧为低速异常, 而东侧为高速异常, 这一近NS向的断裂带与二维剖面的下地壳深断裂在位置上很接近.该地区的几个8级大震均发生在105°E附近, 并且呈一近NS条带.      

Thermoelastic and thermoplastic response of a double-layer porous space containing a decaying heat source

Solutions are presented for the behaviour of a layered porous space which contains a decaying heat source. Such a problem arises when high-level nuclear waste is placed in deep underground depositories in deep clayey formations of sedimentary basins. The geometry of the problem is one dimensional and the porous space is constituted by two layers: a deep low permeability layer which contains the nuclear waste disposal and a superficial layer. The solution is used to examine the effects of contrasts of permeability, thermal conductivity and specific heat capacities between the two layers on the large-scale behaviour of the porous space. Results are presented, using realistic data, for the pore pressure and temperature evolution at the heat source centre, and for the vertical displacement of the ground level. The superficial layer has no significant effects on pore pressure, temperature and stress evolution near the heat source centre. The vertical displacement of the ground level is mainly due to the thermal dilatation of the pore water, so it decreases with an increasing of permeability of the superficial layer. The solution of the time-dependent problem is carried out by applying Laplace transforms to the field variables, obtaining solutions and then using numerical methods to invert the transformed solutions. Comparisons with numerical simulations taking into account the non-linear and non-reversible behaviour of the rock mass are presented. © 1998 John Wiley & Sons, Ltd.     

Upper bound analysis of tunnel face stability in layered soils

The working face of tunnel constructions has to be kept stable during tunneling to prevent large soil deformations or fatal failure. In layered soils with lower cohesion, failures happen more often and more abrupt than in cohesive soils. Therefore, the maintenance of a proper support pressure at the tunnel working face is of high importance. In this paper, an upper bound analysis is introduced to investigate the minimum support pressure for the face stability in layered soils. A three-dimensional kinematically admissible mechanism for the upper bound analysis is improved to model potential failure within different soil layers. An analytical solution for the support pressure assessment is achieved. The influence of the crossing and cover soil on the face stability is analyzed, respectively. This solution provides an analytical estimation of the minimum support pressure for the face stability. It may be used as a reference for projects under similar conditions.     

Tidal signature recorded in burrow fill

The arrangement of sediment couplets preserved in Thalassinoides shafts suggests that tides regulated the passive filling of these trace fossils and, thus, represent tubular tidalites. The thickness variation in individual layers and couplets implies a mixed diurnal, semi‐diurnal tidal signature where packages of either thick‐layered or thin‐layered couplets alternate. Calcarenitic sediment accumulated when tidal current velocity was too high to allow deposition of mud, whereas a marly mud layer is interpreted to have formed during more tranquil times of a tidal cycle (in particular, low‐tide slack water). The tidal record within the burrows covers a few weeks and the corresponding spring–neap cycles. The fill of the Thalassinoides shafts is the only known record to decipher the tidal signature from otherwise totally bioturbated sediments. These deposits accumulated in a lower‐shoreface to upper‐offshore setting during the late Miocene on a shallow shelf extending from the Atlantic Ocean to the west into northern Patagonia. The fill of all investigated burrows started around spring tide and, thus, the behaviour of the burrow producers – probably crustaceans – is speculated to have been affected by tides or the high water level because all studied burrows became abandoned around the same period of a tidal cycle.