海底油气藏及天然气水合物的时频电磁辨识

提出了伪随机多频海洋电磁法观测方案.采用伪随机多频信号作为激励场源,多偏移距同线偶极-偶极同时观测,相关辨识海底地电系统的频率特性和冲激响应,可以在时间域和频率域同时辨识海底高阻薄层.在时间域,瞬变冲激时刻可以直接指示海底地层电导率的变化;在频率域,利用多个频率的电场响应计算的频散率及其道闻变化量,相对相位道间变化量对高阻薄层有很好的反映.从而实现对海底油气及天然气水合物的多参数辨识.     

Influence of seismicity parameter uncertainty on seismic hazard estimation of cities and towns

Introduction In his Doctor thesis, PAN (2000) has roughly studied some about the influence of parameter uncertainty on seismic hazard analysis. The aim of this paper is especially to study the uncertainty influence extent of seismic hazard analysis result in a spot, owing to the uncertainty of seismicity parameters induced by variational statistical time ranges of various statistical regions in the whole country. The authors have inquired into the uncertain extents of influence of seismic p…     

Seismic preprocessing and amplitude cross-calibration for a time-lapse amplitude study on seismic data from the Oseberg reservoir

The cross‐calibration of different vintage data is an important prerequisite in attempting to determine the time‐lapse seismic effects induced by hydrocarbon production in a reservoir. This paper reports the preprocessing and cross‐calibration procedures adopted to modify the data of four seismic vintages (1982, 1989, 1992 and 1999) from the Oseberg field in the North Sea, for optimal conditions for a time‐lapse seismic amplitude analysis. The final results, in terms of time‐lapse variations, of acoustic impedance and of amplitude‐versus‐offset, are illustrated for selected data sets. The application of preprocessing to each individual vintage data set reduces the effects of the different acquisition and noise conditions, and leads to consistency in the amplitude response of the four vintages. This consistency facilitates the final amplitude cross‐calibration that is carried out using, as reference, the Cretaceous horizon reflections above the Brent reservoir. Such cross‐calibration can be considered as vintage‐consistent residual amplitude correction. Acoustic impedance sections, intercept and gradient amplitude‐versus‐offset attributes and coherent amplitude‐versus‐offset estimates are computed on the final cross‐calibrated data. The results, shown for three spatially coincident 2D lines selected from the 1982, 1989 and 1999 data sets, clearly indicate gas‐cap expansion resulting from oil production. Such expansion is manifested as a decrease in acoustic impedance and a modification of the amplitude‐versus‐offset trends in the apical part of the reservoir.     

The impact of reservoir scale on amplitude variation with offset

To investigate the vertical and horizontal impact of reservoir scale on the amplitude‐versus‐offset characteristics, we conduct seismic numerical simulations on models containing spatially confined lithologic units with different scales. We find that the reservoir scale has a nonlinear effect on the amplitude‐versus‐offset intercepts and gradients. As the reservoir width increases, amplitude‐versus‐offset intercept and gradient both first increase, then decrease, and finally remain stable. The amplitude‐versus‐offset intercept is maximum when the reservoir width is 80% larger than the Fresnel zone radius, whereas the amplitude‐versus‐offset gradient peaks at 1.5 times the Fresnel zone radius. Both amplitude‐versus‐offset intercept and gradient are approximately proportional to the reservoir width prior to reaching their maxima. When the lateral extent of the reservoir is more than three times the Fresnel zone radius, the amplitude‐versus‐offset attributes are constant. Modelling the reservoir thickness shows that intercept and gradient behave in a manner similar to that of tuning of thin beds. Both the amplitude‐versus‐offset intercept and gradient first increase and then decrease with the thickness, peaking at the tuning thickness. The thickness contribution to amplitude‐versus‐offset variations is negligible when the thickness is larger than 1.6 times of the tuning thickness. Considering the magnitude of the changes in amplitude‐versus‐offset intercept and gradient caused by reservoir scale, the width causes a maximum 433% intercept increase and a 344% gradient increase, whereas the thickness causes a maximum 100% intercept increase and a 73% gradient increase. Cross‐plotting the amplitude‐versus‐offset intercept and gradient shows the reservoir scale change gives rise to an anti‐clockwise spiraling effect. In conclusion, the lateral and vertical extents of the reservoir both play an important role in amplitude variation with offset. Our analysis shows that the lateral reservoir extent has a larger impact on the amplitude variation with offset than the vertical tuning effect.     

现代功率谱估算在反射地震剖面随机介质参数估算中的应用

由于地壳深部物质呈现出多尺度的非均质性,导致了地震剖面上反射能量弱、同相轴连续性差、多成带状分布.传统的反射地震解释方法难以有效应用到深反射地震资料解释中.统计分析的方法为深反射地震资解释提供了新的工具,地下介质这种随机分布的非均质性可以借助随机介质模型来进行描述.深反射地震剖面继承了地下介质的非均质性,因此可以通过从反射地震剖面估算随机介质参数来研究地下介质非均质性的空间分布.由于地震剖面长度有限,随机介质参数估算中自相关函数计算误差较大,导致了估算的不准确.现代功率谱估算克服了短数据自相关函数难以准确求取的问题,提高了参数估算准确性.模型计算和实际数据测算都证明了该方法能有效提高自相关函数计算精度和随机介质参数剖面的分辨率.     

用地震资料估算珊溪水库岩层孔隙度和饱和度

本文首先从室内岩石力学实验数据出发,利用Biot-Gassmann方程和Pride模型建立的岩石骨架与岩石基质体积模量函数关系式求出珊溪水库区岩石基质模量和固结系数,然后分别利用联合地震波走时方程组和多台和达曲线求得P波速度和波速比,最后以流纹斑岩为例,计算了岩石饱和度、孔隙度,并分析了地震中造成的饱和度和孔隙度的变化。结果表明:(1)地震过程中岩石饱和度均大于99%,且大多为100%,即地震中岩石处于水饱和或接近水饱和状态,波速和波速比的变化主要由岩石孔隙度的变化引起。(2)珊溪水库震中区流纹斑岩孔隙度的上限值为6.6%~15.5%,下限值为3.6%~4.2%。下限值与地震前由实验室测得的孔隙度相近。如果实验室测得的流纹斑岩孔隙度是其初始值,则地震中岩石孔隙度增大了3.2%~12.09%。(3)2008年1月地震以后震中区岩石孔隙度逐渐闭合,直到2014年8月孔隙度开始迅速增大,9月15日起原震区再次发生一组新的地震活动,震群活动前岩石出现了扩容现象。     

利用气枪地震资料研究燕山隆起带南部地区地壳S波速度及泊松比结构

泊松比是了解地球内部介质的一个重要参数,在地震学上可以通过P波和S波速度的比值来确定.但是传统人工地震测深采用的炸药等爆破源产生的S波较弱,难以用来研究泊松比结构.对2006年河北省遵化市上关湖水库大容量气枪试验产生的地震数据进行研究,发现气枪震源能够有效地产生S波.根据偏振分析,认为S波主要是气枪产生的P波在水库底部固液界面转换产生的.波形拟合能够充分利用波形和走时信息,是研究地下结构的一种有效的方法.对气枪产生的S波进行波形拟合,获得了燕山隆起带南部地区地壳的S波速度模型,并对比修正后的P波速度模型,进一步得到了该地区的泊松比结构.研究表明:1)该地区的地壳厚度约为33km.2)该地区地壳整体泊松比值偏低,上(0～14km)中(14～28km)下(28～33km)地壳以及上地幔顶部(33km附近)的泊松比平均值分别为0.25,0.25,0.27,0.27.上地壳和中地壳可能以长英质的酸性岩石为主,下地壳和上地幔顶部可能以铁镁质的基性岩石和中性岩石为主.3)该地区上下地壳各存在一个低速层,它们可能主要是流体作用的结果.     

基于振幅加权的拖缆与海底地震资料联合波形反演方法

南黄海新近系沉积层与其下方中-古生代海相地层之间速度差异大,存在强波阻抗界面,导致接收的浅部和中深部地层的反射波振幅差异极大.一般全波形反演方法未考虑不同界面反射波振幅的差异,对南黄海浅部沉积层速度的反演效果好,而对中深部地层速度的反演效果差.为此,本文提出一种基于振幅加权的拖缆与海底地震资料联合波形反演方法.该方法对中深部弱反射波数据施加较大的权重,对浅部强反射波数据施加较小的权重,从而均衡浅、中和深部地层反射波的振幅,使得自浅至深所有地层速度均能得到较好反演效果.同时,该方法利用海底地震资料和拖缆地震记录积分的低频成分反演低波数背景速度场,利用拖缆资料的高频成分提高反演模型的分辨率,充分发挥了拖缆和海底地震资料的优势互补作用,获得了比仅使用单一资料更好的反演结果.根据地震剖面和钻井资料,建立了南黄海中部隆起的速度模型.利用该理论模型的合成地震记录对本文方法进行了测试和分析,验证了本文方法的正确性和有效性.

     

Geomechanical property estimation of unconventional reservoirs using seismic data and rock physics

An extension of a previously developed rock physics model is made that quantifies the relationship between the ductile fraction of a brittle/ductile binary mixture and the isotropic seismic reflection response. By making a weak scattering (Born) approximation and plane wave (eikonal) approximation, with a subsequent ordering according to the angles of incidence, singular value decomposition analyses are performed to understand the stack weightings, number of stacks, and the type of stacks that will optimally estimate two fundamental rock physics parameters – the ductile fraction and the compaction and/or diagenesis. It is concluded that the full PP stack, i.e., sum of all PP offset traces, and the “full” PS stack, i.e., linear weighted sum of PS offset traces, are the two optimal stacks needed to estimate the two rock physics parameters. They dominate over both the second‐order amplitude variation offset “gradient” stack, which is a quadratically weighted sum of PP offset traces that is effectively the far offset traces minus the near offset traces, and the higher order fourth order PP stack (even at large angles of incidence). Using this result and model‐based Bayesian inversion, the seismic detectability of the ductile fraction (shown by others to be the important rock property for the geomechanical response of unconventional reservoir fracking) is demonstrated on a model characteristic of the Marcellus shale play.     

Dual state–parameter estimation of hydrological models using ensemble Kalman filter

Hydrologic models are twofold: models for understanding physical processes and models for prediction. This study addresses the latter, which modelers use to predict, for example, streamflow at some future time given knowledge of the current state of the system and model parameters. In this respect, good estimates of the parameters and state variables are needed to enable the model to generate accurate forecasts. In this paper, a dual state–parameter estimation approach is presented based on the Ensemble Kalman Filter (EnKF) for sequential estimation of both parameters and state variables of a hydrologic model. A systematic approach for identification of the perturbation factors used for ensemble generation and for selection of ensemble size is discussed. The dual EnKF methodology introduces a number of novel features: (1) both model states and parameters can be estimated simultaneously; (2) the algorithm is recursive and therefore does not require storage of all past information, as is the case in the batch calibration procedures; and (3) the various sources of uncertainties can be properly addressed, including input, output, and parameter uncertainties. The applicability and usefulness of the dual EnKF approach for ensemble streamflow forecasting is demonstrated using a conceptual rainfall-runoff model.     

Characterization of a carbonate reservoir using elastic full-waveform inversion of vertical seismic profile data

Two-dimensional elastic full waveform inversion was applied to two lines extracted from a spiral three-dimensional vertical seismic profile data acquired in an oilfield offshore, Abu Dhabi, in the United Arab Emirates. The lines were selected to be parallel and perpendicular to the plane defined by the deviated borehole. The purpose of the inversion was to derive high-resolution elastic properties of the subsurface. After pre-processing, the data were band-pass filtered with a minimum frequency of 3.5 Hz and a maximum frequency of 30 Hz. A sequential inversion approach was used to mitigate non-linearity. The pre-processing of the data consisted in the removal of bad traces, followed by amplitude and phase corrections. High-resolution P- and S-wave velocity models that show good correlations with the available sonic logs were obtained. The results of the inversion suggest that the oilfield consists of a stack of layers with varying lithology, porosity and possibly fluid content.     

基于地震活动性资料估计海原断裂倾角

本文以海原断裂带区域活动构造为基础,将海原断裂划分为西、中、东3段.基于1999年12月26日至2010年7月26日间的精定位小震目录,估计了海原断裂带各段的倾角.考虑海原断裂沿走向可能存在南倾与北倾两种情况,将倾角的范围设置为.首先运用网格搜索法确定了平面断层模型,其次以特征深度节点为基础数据运用多项式构建了曲面断层模型.结果表明：当进行平面拟合时,海原断裂西段与中段、东段的倾向不同,西段为南倾,其倾角值为71°,而中段、东段为北倾,其倾角值分别为72°、65°,各段的倾角值均由地表以下8 km地震资料确定.当进行曲面拟合时,在8 km深度以内海原断裂西段、中段、东段的倾角均处于80°左右,即接近陡立.西段的倾角在深度为9 km处出现转换,之后倾角接近陡立;中段的倾角在深度为16 km处出现转换,之后倾角逐渐减小,当深度为18 km时倾角为30°;东段的倾角在深度为11 km处出现转换,倾角为42°,在深度为16 km处出现第二次转换,倾角为55°,之后倾角逐渐减小.结合震源机制解和大地测量观测资料反演拟合的合理性,验证了本文所估计倾角的可靠性.     

analyzing the variation characteristics of stress field in hetao seismic belt using focal mechanism data

Based on analysis of background of geological tectonic movement and strong earthquake activity, we first obtained the focal mechanism solutions using amplitude ratio and CAP method, then determined the characteristic of average stress field of the study area by inversion of the stress field. On this basis, we selected the source mechanism consistency parameter as the inspection index to obtain the latest changes of stress field in Hetao seismic zone based on its temporal and spatial analysis. Two methods were used in the stress field inversion for comparison and analysis, which are average stress axis tensor and LSIB(Linear stress inversion bootstrap, LSIB). According to the geological tectonic movement and focal mechanism solutions of M_S≥4.0 earthquakes from 1970, we judge that the stress field evolution process of Hetao seismic belt is controlled jointly by vertical difference movement and horizontal shear movement, resulting in that the normal fault and strike-slip fault mechanisms are dominating.Taking into account the station layout of the study area, and in order to ensure the accuracy of calculation, we calculated 224 earthquakes focal mechanism solutions by using amplitude ratio and CAP method, including 164 earthquakes with 2.8≤M_L<3.5, 42 earthquakes with 3.5≤M_L<4.0, and 18 earthquakes with M_L≥4.0; The statistical results on type of focal mechanisms show that, there are 142 strike-slip earthquakes(63.4%), 50 normal fault earthquakes(22.3%)and 32 thrust fault earthquakes(14.3%). In this study period(from 2001 to 2012), most earthquakes had a strike-slip mechanism in Hetao seismic belt, this is one of the inherent characteristics of the stress field.The result of average stress axis tensor and LSIB shows that, the azimuth of maximum compressional stress is 47°～52°, direction is NE-SW; The azimuth of minimum compressional stress is 313°～322°, direction is NW-SE; This indicates that, the stress field characteristics of Hetao seismic belt and its sub-block are not completely consistent. Linhe Basin exhibits coordinated stress field characteristics with Hetao seismic belt, but Hubao Basin exhibits regional differences, direction of compressive stress has clockwise deflection in Baotou area, and the compressive stress direction is NEE. This heteropical character of stress field is also confirmed by horizontal projection distribution of stress axis of historical strong earthquakes and recent moderate and small earthquakes. Since 2003, the temporal sequence curve of consistency parameter of Hetao seismic belt had a downward trend, this change was caused by focal mechanism consistency parameter of Linhe to Wuhai area, which indicates that this structural position is possible to be a priority area for stress accumulation and accelerated release in future.     

The influence of sea water velocity variation on seismic traveltimes, ray paths, and amplitude

The main factors affecting seismic exploration is the propagation velocity of seismic waves in the medium. In the past, during marine seismic data processing, the propagation velocity of sea water was generally taken as a constant 1500 m/s. However, for deep water exploration, the sound velocity varies with the season, time, location, water depth, ocean currents, and etc.. It also results in a layered velocity distribution, so there is a difference of seismic traveltime, ray paths, and amplitude, which affect the migration imaging results if sea water propagation velocity is still taken as constant for the propagation wavefield. In this paper, we will start from an empirical equation of seismic wave velocity in seawater with changes of temperature, salinity, and depth, consider the variation of their values, build a seawater velocity model, and quantitatively analyze the impact of seawater velocity variation on seismic traveltime, ray paths, and amplitude in the seawater velocity model.     

Bayesian inversion of time‐lapse seismic data for the estimation of static reservoir properties and dynamic property changes

Seismic conditioning of static reservoir model properties such as porosity and lithology has traditionally been faced as a solution of an inverse problem. Dynamic reservoir model properties have been constrained by time‐lapse seismic data. Here, we propose a methodology to jointly estimate rock properties (such as porosity) and dynamic property changes (such as pressure and saturation changes) from time‐lapse seismic data. The methodology is based on a full Bayesian approach to seismic inversion and can be divided into two steps. First we estimate the conditional probability of elastic properties and their relative changes; then we estimate the posterior probability of rock properties and dynamic property changes. We apply the proposed methodology to a synthetic reservoir study where we have created a synthetic seismic survey for a real dynamic reservoir model including pre‐production and production scenarios. The final result is a set of point‐wise probability distributions that allow us to predict the most probable reservoir models at each time step and to evaluate the associated uncertainty. Finally we also show an application to real field data from the Norwegian Sea, where we estimate changes in gas saturation and pressure from time‐lapse seismic amplitude differences. The inverted results show the hydrocarbon displacement at the times of two repeated seismic surveys.     

利用中国北部连续重力站资料评估全球潮汐模型精度

为评估全球潮汐模型在我国潮汐改正中的适用性,本文首先对10个重力站2016—2018年的观测数据进行了精度评定,而后基于均方根、和方根、纬度依赖关系以及重力残差等指标对7个全球潮汐模型进行了精度评定.结果表明:10个重力站的一些评价指标达到甚至超越了早期超导重力仪,例如M2波潮汐因子的中误差普遍小于0.00070,其中...     

松辽盆地徐深气田火山岩气藏地质约束-测井地震建模方法探讨

本文在借鉴徐深气田以往地质建模成功经验基础上,针对徐深气田B区块火山体规模小、横向变化快、开发井少、认识程度较低的实际,密切结合火山岩成因特点,探索性地提出体控地质建模的新思路,即以地震反射特征和地震属性特征为依据,结合钻井、测井资料,井震结合从三维空间上精细识别和刻画火山岩体,通过利用神经网络算法建立多种地震属性与测井参数的定量关系,实现了火山岩体内部多种地质属性的定量表征,建立的三维可视化地质模型比较真实地反映了火山机构及其内部属性空间发育和分布特征.     

Reconstruction of the layer anisotropic elastic parameters and high-resolution fracture characterization from P-wave data: a case study using seismic inversion and Bayesian rock physics parameter estimation

Wide-azimuth seismic data can be used to derive anisotropic parameters on the subsurface by observing variation in subsurface seismic response along different azimuths. Layer-based high-resolution estimates of components of the subsurface anisotropic elastic tensor can be reconstructed by using wide-azimuth P-wave data by combining the kinematic information derived from anisotropic velocity analysis with dynamic information obtained from amplitude versus angle and azimuth analysis of wide-azimuth seismic data. Interval P-impedance, S-impedance and anisotropic parameters associated with anisotropic fracture media are being reconstructed using linearized analysis assuming horizontal transverse anisotropy symmetry. In this paper it is shown how additional assumptions, such as the rock model, can be used to reduce the degrees of freedom in the estimation problem and recover all five anisotropic parameters. Because the use of a rock model is needed, the derived elastic parameters are consistent with the rock model and are used to infer fractured rock properties using stochastic rock physics inversion. The inversion is based on stochastic rock physics modelling and maximum a posteriori estimate of both porosity and crack density parameters associated with the observed elastic parameters derived from both velocity and amplitude versus angle and azimuth analysis. While the focus of this study is on the use of P-wave reflection data, we also show how additional information such as shear wave splitting and/or anisotropic well log data can reduce the assumptions needed to derive elastic parameter and rock properties.     

Extended Kalman filter for material parameter estimation in nonlinear structural finite element models using direct differentiation method

This paper presents a novel nonlinear finite element (FE) model updating framework, in which advanced nonlinear structural FE modeling and analysis techniques are used jointly with the extended Kalman filter (EKF) to estimate time‐invariant parameters associated to the nonlinear material constitutive models used in the FE model of the structural system of interest. The EKF as a parameter estimation tool requires the computation of structural FE response sensitivities (total partial derivatives) with respect to the material parameters to be estimated. Employing the direct differentiation method, which is a well‐established procedure for FE response sensitivity analysis, facilitates the application of the EKF in the parameter estimation problem. To verify the proposed nonlinear FE model updating framework, two proof‐of‐concept examples are presented. For each example, the FE‐simulated response of a realistic prototype structure to a set of earthquake ground motions of varying intensity is polluted with artificial measurement noise and used as structural response measurement to estimate the assumed unknown material parameters using the proposed nonlinear FE model updating framework. The first example consists of a cantilever steel bridge column with three unknown material parameters, while a three‐story three‐bay moment resisting steel frame with six unknown material parameters is used as second example. Both examples demonstrate the excellent performance of the proposed parameter estimation framework even in the presence of high measurement noise. Copyright © 2015 John Wiley & Sons, Ltd.     

Deep exploration using long-offset transient electromagnetics: interpretation of field data in time and frequency domain

In the framework of the Deep Electromagnetic Soundings for Mineral Exploration project, we conducted ground-based long-offset transient-electromagnetic measurements in a former mining area in eastern Thuringia, Germany. The large-scale survey resulted in an extensive dataset acquired with multiple high-power transmitters and a high number of electric and magnetic field receivers. The recorded data exhibit a high data quality over several decades of time and orders of magnitude. Although the obtained subsurface models indicate a strong multi-dimensional subsurface with variations in resistivity over three orders of magnitude, the electrical field step-on transients are well fitted using a conventional one-dimensional inversion. Due to superimposed induced polarization effects, the transient step-off data are not interpretable with conventional electromagnetic inversion. For further interpretation in one and two dimensions, a new approach to evaluate the long-offset transient-electromagnetic data in frequency domain is realized. We present a detailed workflow for data processing in both domains and give an overview of technical obstructions that can occur in one domain or the other. The derived one-dimensional inversion models of frequency-domain data show strong multi-dimensional effects and are well comparable with the conventional time domain inversion results. To adequately interpret the data, a 2.5D frequency-domain inversion using the open source algorithm MARE2DEM (Modeling with Adaptively Refined Elements for 2-D EM) is carried out. The inversion leads to a consistent subsurface model with shallow and deep conductive structures, which are confirmed by geology and additional geophysical surveys.