Determination of Q S, spreading factor and S-wave acceleration by inversion of seismic intensity data

In this paper, on the basis of intensity data from 85 earthquakes occurred in China, the method of resolving overdetermined equations by using the damped least squares method is applied to inversing for the mean value of 1 s. S wave quality factor, the spreading factor and the S-wave acceleration at the foci of earthquakes in 7 regions of China. The relation between S-wave acceleration at the foci of earthquakes and magnitude is discussed. As an example, 2-DQ _S distribution in Taiwan Province is obtained by inversion. It is found that there exists some corresponding relation between this distribution and Bouguer gravitational anomaly and seismicity. Preliminary analysis and discussion on results of the inversion, and estimation of error inQ _S are made. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 202–211, 1991. The English version of this paper is improved by Professor Yushou Xie.     

S-wave splitting intensity analysis and inversion

Analysing S-wave splitting has become a routine step in processing multicomponent data. Typically, this analysis leads to determining the principal directions of a transversely isotropic medium with a horizontal symmetry axis, which is assumed to be responsible for azimuthal anisotropy, and to the time delays between the fast and slow S-waves. These parameters are commonly estimated layer-by-layer from the top. Errors in layer stripping occurring in shallow layers might propagate to deeper layers. We propose a method for S-wave splitting analysis and compensation that consists of inverting interval values of splitting intensity to obtain a model of anisotropic parameters that vary with time and/or depth. Splitting intensity is a robust attribute with respect to structural variations and is commutative, which means that it can be summed along a ray (or throughout a sensitivity kernel volume) and can be linearly related to anisotropic perturbations at depth. Therefore, it is possible to estimate anisotropic properties within a geological formation (e.g. the reservoir) by analysing the differences of splitting intensity measured at the top and at the bottom of the layer. This allows us to avoid layer stripping, in particular, for shallow layers where anisotropic parameters are difficult to estimate due to poor coverage, and it makes S-wave splitting analysis simpler to apply. We demonstrate this method on synthetic and real data. Because the splitting intensity attribute shows usefulness in S-wave splitting analysis in transversely isotropic media, we extend the splitting intensity theory to lower symmetry classes. It enables the characterization of tilted transversely isotropic and tilted orthorhombic media, opening new opportunities for anisotropic model building.     

地震数据互相关驱动的多道反演方法

地震反演是储层定量描述和地震油气识别的关键技术,反演结果在复杂构造区域的横向连续性和保真性是影响地震资料定量解释精度的重要因素.基于此,本文发展了地震数据互相关驱动的多道反演方法.考虑地层反射系数与地震数据在结构上具有相似性的特点,基于地震数据互相关描述地层反射系数的结构特征,并将其作为多道地震反演的横向约束条件;此外,为改善地震数据本身横向连续性差对反演结果的影响,在目标泛函的惩罚项中引入局部优化算子,构建了一个易于求解的多道地震反演目标泛函.与常规多道地震反演方法相比,本文方法能够设计更合理、更符合实际情况的横向约束算子,提高反演结果的横向连续性,并且能有效降低地震资料质量对反演结果的影响.模型测试和实际应用验证了本方法的可靠性和稳定性.

     

地震数据互相关驱动的多道反演方法

地震反演是储层定量描述和地震油气识别的关键技术,反演结果在复杂构造区域的横向连续性和保真性是影响地震资料定量解释精度的重要因素.基于此,本文发展了地震数据互相关驱动的多道反演方法.考虑地层反射系数与地震数据在结构上具有相似性的特点,基于地震数据互相关描述地层反射系数的结构特征,并将其作为多道地震反演的横向约束条件;此外,为改善地震数据本身横向连续性差对反演结果的影响,在目标泛函的惩罚项中引入局部优化算子,构建了一个易于求解的多道地震反演目标泛函.与常规多道地震反演方法相比,本文方法能够设计更合理、更符合实际情况的横向约束算子,提高反演结果的横向连续性,并且能有效降低地震资料质量对反演结果的影响.模型测试和实际应用验证了本方法的可靠性和稳定性.     

反射地震数据的逐层波形反演

本文针对层状介质并结合梯度法波形反演,提出逐层波形反演的方法. 首先给出介质扰动响应的概念,并在此基础上分析了梯度法波形反演方法. 波形反演实质上是将实测地震记录和预测地震记录的波形残差信息转化为实际地质模型和预测地质模型的模型残差信息. 波形反演的优点是利用大量振幅相位信息得到高分辨率的反演结果, 其缺点是运行耗时大;当初始模型和实际模型相差较大时,迭代算法容易陷入局部极小点,这是因为目标函数和初始模型同实际模型间的差异是非线性的关系. 逐层波形反演方法是使自上而下每一层的目标函数最小,这样总的目标函数也是最小的. 利用二分法速度扫描确定每一层速度不仅提高了运算速度也避免了迭代算法陷入局部极小点的问题. 结合介质扰动响应和目标函数值变化可以更为准确迅速地确定每一层速度和该层界面位置.     

从爆破烈度看仪器烈度的算法

以某次工程爆破的地震动测试结果为例,对《中国地震烈度表》、日本气象厅JMA系统、美国ShakeMap系统、福建省地方标准以及袁一凡教授提出的烈度算法进行了对比分析,讨论了以地震动参数为基础的仪器烈度计算方法,提出了将袁一凡仪器烈度算法中的PGA改为水平分向合成峰值,忽略持时和主频两个参数并修改个各参数权重的意见。经分析认为,袁一凡仪器烈度算法和福建省地方标准中的仪器烈度算法在处理大中小震时结果较为合理,且袁一凡仪器烈度算法保留了继续研究地震烈度机理的学术空间和扩展应用领域的可能。     

地震密度反演及地层孔隙度估计

地层密度直接与孔隙度、孔隙流体类型、饱和度和骨架矿物成分有关.本文通过理论分析和计算,讨论了油气藏储层物性参数变化引起的密度变化及密度变化对地震波速度、阻抗和振幅的影响,提出了基于完全纵波方程的全波形地震密度反演和孔隙度估计方法,克服了常规地震密度反演对地震数据更多处理引起的信号畸变,提高了地震密度反演和地层孔隙度估计的精度.该方法采用波场导数的时间积分和多炮求和,对地震数据中的噪声具有比较强的压制作用.理论模型研究表明该方法是可行的.通过对我国西部某气田实际数据处理、分析和反演,获得了地层密度和孔隙度,结果与测井基本吻合,证明了预测结果的准确性和方法的有效性,从而为后续的有效储层预测和储量计算提供了可靠的数据.     

Constructing a discrete fracture network constrained by seismic inversion data

Rock fractures are of great practical importance to petroleum reservoir engineering because they provide pathways for fluid flow, especially in reservoirs with low matrix permeability, where they constitute the primary flow conduits. Understanding the spatial distribution of natural fracture networks is thus key to optimising production. The impact of fracture systems on fluid flow patterns can be predicted using discrete fracture network models, which allow not only the 6 independent components of the second‐rank permeability tensor to be estimated, but also the 21 independent components of the fully anisotropic fourth‐rank elastic stiffness tensor, from which the elastic and seismic properties of the fractured rock medium can be predicted. As they are stochastically generated, discrete fracture network realisations are inherently non‐unique. It is thus important to constrain their construction, so as to reduce their range of variability and, hence, the uncertainty of fractured rock properties derived from them. This paper presents the underlying theory and implementation of a method for constructing a geologically realistic discrete fracture network, constrained by seismic amplitude variation with offset and azimuth data. Several different formulations are described, depending on the type of seismic data and prior geologic information available, and the relative strengths and weaknesses of each approach are compared. Potential applications of the method are numerous, including the prediction of fluid flow, elastic and seismic properties of fractured reservoirs, model‐based inversion of seismic amplitude variation with offset and azimuth data, and the optimal placement and orientation of infill wells to maximise production.     

叠前地震数据特征波场分解、偏移成像与层析反演

本文提出了一套叠前地震数据稀疏表达(特征波场合成)、深度偏移成像和层析成像的处理流程.不同于传统的变换域中的数据稀疏表达理论, 本文利用局部平面波的传播方向(慢度矢量), 在中心炮检点处同时进行波束合成, 从而将地震数据投影到局部平面波域(高维空间)中.由于波束合成后的地震数据描述了局部平面波的方向特征, 因此称之为特征波场.然而波束合成算法需要估计局部平面波的慢度矢量.当地震数据受噪声干扰时, 难以在常规τ-p谱中自动估计局部平面波的射线参数(慢度矢量).本文提出了基于反演理论的特征波场合成方法, 可以同时反演局部平面波及其传播方向, 从而提高特征波合成的自动化程度并保持方法的稳健性.通过特征波场合成, 可以将地震数据分解为单独的震相(波形).这样的数据可以直接用来成像及反演.在局部平面波域中, 由于局部平面波的入射与出射射线参数已知, 传统的Kirchhoff叠前深度偏移(PSDM)和高斯束/控制束PSDM可以实现从"沿等时面的画弧"到"向反射点(段)的直接投影"的转变, 叠前偏移的效率以及成像质量可以同时提高.此外, 特征波场与地下反射点(段)的一对一映射关系使得叠前深度偏移与层析成像融为一体, 可以极大地提高速度反演的效率.数值试验证明了特征波场合成、叠前深度成像以及层析反演的有效性.     

Determination of a shallow velocity–depth model from seismic refraction data by coherence inversion1

Seismic refractions have different applications in seismic prospecting. The travel- times of refracted waves can be observed as first breaks on shot records and used for field static calculation. A new method for constructing a near-surface model from refraction events is described. It does not require event picking on prestack records and is not based on any approximation of arrival times. It consists of the maximization of the semblance coherence measure computed using shot gathers in a time window along refraction traveltimes. Time curves are generated by ray tracing through the model. The initial model for the inversion was constructed by the intercept-time method. Apparent velocities and intercept times were taken from a refraction stacked section. Such a section can be obtained by appling linea moveout corrections to common-shot records. The technique is tested successfully on synthetic and real data. An important application of the proposed method for solving the statics problem is demonstrated.     

地震有效峰值加速度与地震烈度相关性研究

利用2008年5月12日发生的汶川8·0级以及2008年8月30日凉山州会理县—攀枝花市间6·1级地震中获取的四川省境内的地震加速度记录数据,分析了有效峰值加速度（EPA）与加速度峰值（PGA）对比结果及加速度数据有效值与地震烈度的相关性。得出了有效峰值加速度（EPA）与烈度的相关性较之地震峰值加速度（PGA）与烈度的...     

Determination of porosity and saturation using seismic waveform inversion

Characterization of a reservoir model requires determination of its petrophysical parameters, such as porosity and saturation. We propose a new method to determine these parameters directly from seismic data. The method consists of the computation and inversion of seismic waveforms. A high frequency method is presented to model wave propagation through an attenuative and dispersive poroelastic medium. The high frequency approximation makes it possible to efficiently compute sensitivity functions. This enables the inversion of seismic waveforms for porosity and saturation. The waveform inversion algorithm is applied to two laboratory crosswell datasets of a water saturated sand. The starting models were obtained using travel time tomography. The first dataset is inverted for porosity. The misfit reduction for this dataset is approximately 50%. The second dataset was obtained after injection of a nonaqueous-phase liquid (NAPL), possibly with some air, which made the medium more heterogeneous. This dataset was inverted for NAPL and air saturation using the porosity model obtained from the first inversion. The misfit reduction of the second experiment was 70%. Regions of high NAPL and high air saturation were found at the same location. These areas correlate well with the position of one of the injection points as well as regions of higher NAPL concentrations found after excavation of the sand. It is therefore possible to directly invert waveforms for pore fluid saturation by taking into account the attenuation and dispersion caused by the poroelasticity.     

P-P波和P-SV波联合AVO反演研究

目前转换波AVO(Amplitude Versus Offset)反演或纵波、转换波联合AVO反演是基于Aki等的反射系数近似公式,通常反演结果为上、下层速度比值和速度平均比值.本文在Shuey和郑晓东近似公式的基础上,运用最小平方算法,提出一种p-p波和P-SV波联合AVO反演下层介质纵、横波速度及泊松比的方案,并应用于双层含气砂岩模型上,误差在10%以内.此方案要求上、下层物性差别不大且入射角小于30°,但不需要速度比值为常数2的假设.     

地震与大地电磁测深数据的自适应正则化同步联合反演

在传统的联合反演研究中,地球物理学者往往更多地关注数据拟合,很少涉及正则化理论.本文在电阻率和速度随机分布的大地电磁测深(MT)与地震联合反演研究的基础之上,将正则化思想引入到同步联合反演中,加入先验信息进行模型约束,选取最小模型为稳定泛函,并首次采用自适应正则化算法来确定联合反演的正则化因子.根据以往研究成果,采用非线性模拟退火方法来实现MT视电阻率或相位与地震走时的同步联合反演.此外,为了验证该算法的有效性,在模型对比试验中设计了4种不同方案.通过模型试验的对比分析,我们认为加入有效模型约束的自适应正则化联合反演,可以有效地提高解的稳定性和计算效率,并能在一定程度上解决不同地球物理数据加权系数人为选取问题,模型试验结果也表明了自适应正则化联合反演优于MT单独反演.     

测井、井间地震与地面地震数据联合约束下的地质统计学随机建模方法研究

利用能够整合测井信息与井间地震信息的地质统计学随机模拟方法,结合传统的地质统计学反演思路,得到了一种能够同时整合测井、井间地震与地面地震三种先验信息的地质统计学反演与储层建模方法.由于井间射线信息、测井信息与地面地震数据在随机反演与建模过程当中都得到了尊重,因此与传统地质统计学反演仅利用了测井与地面地震数据相比,本文的地质统计学反演与建模方法更充分地利用了先验信息,有效提高了反演的精度,降低了随机建模中的多解性.基于理论数据的测试证实了上述观点.     

测井、井间地震与地面地震数据联合约束下的地质统计学随机建模方法研究

利用能够整合测井信息与井间地震信息的地质统计学随机模拟方法,结合传统的地质统计学反演思路,得到了一种能够同时整合测井、井间地震与地面地震三种先验信息的地质统计学反演与储层建模方法.由于井间射线信息、测井信息与地面地震数据在随机反演与建模过程当中都得到了尊重,因此与传统地质统计学反演仅利用了测井与地面地震数据相比,本文的地质统计学反演与建模方法更充分地利用了先验信息,有效提高了反演的精度,降低了随机建模中的多解性.基于理论数据的测试证实了上述观点.     

A fast joint seismic data reconstruction by sparsity‐promoting inversion

Seismic field data are often irregularly or coarsely sampled in space due to acquisition limits. However, complete and regular data need to be acquired in most conventional seismic processing and imaging algorithms. We have developed a fast joint curvelet‐domain seismic data reconstruction method by sparsity‐promoting inversion based on compressive sensing. We have made an attempt to seek a sparse representation of incomplete seismic data by curvelet coefficients and solve sparsity‐promoting problems through an iterative thresholding process to reconstruct the missing data. In conventional iterative thresholding algorithms, the updated reconstruction result of each iteration is obtained by adding the gradient to the previous result and thresholding it. The algorithm is stable and accurate but always requires sufficient iterations. The linearised Bregman method can accelerate the convergence by replacing the previous result with that before thresholding, thus promoting the effective coefficients added to the result. The method is faster than conventional one, but it can cause artefacts near the missing traces while reconstructing small‐amplitude coefficients because some coefficients in the unthresholded results wrongly represent the residual of the data. The key process in the joint curvelet‐domain reconstruction method is that we use both the previous results of the conventional method and the linearised Bregman method to stabilise the reconstruction quality and accelerate the recovery for a while. The acceleration rate is controlled through weighting to adjust the contribution of the acceleration term and the stable term. A fierce acceleration could be performed for the recovery of comparatively small gaps, whereas a mild acceleration is more appropriate when the incomplete data has a large gap of high‐amplitude events. Finally, we carry out a fast and stable recovery using the trade‐off algorithm. Synthetic and field data tests verified that the joint curvelet‐domain reconstruction method can effectively and quickly reconstruct seismic data with missing traces.     

利用地震烈度数据重构历史地震震级和震中位置的方法

采用现代仪器地震记录和烈度资料,我们建立了用于定量估算历史地震震级和震中位置的烈度衰减模型和分析方法,并对我国西南川滇地区的历史和近代地震的震中和震级进行了重新分析.利用20世纪以来该地区十四个有仪器记录的地震（5.9≤M_s≤8.0）及相应的烈度数据,对其烈度-震级-震中距衰减关系进行标定,并建立了用于震中和震级估算的烈度衰减模型.结果表明,当震级一定时,川滇地区烈度随震中距增大而衰减的速率明显小于美国加州地区（～60%）.在衰减模型基础之上,发展了确定震中区域和震级的网格搜索试算方法（GSTSL）,并给出了确定震中位置和震级的等值线置信值.采用所建立的分析方法,对1786年康定地震,1850年西昌地震,1913年峨山地震和1970年通海地震进行了分析,给出了更为精确合理的结果.     

基于地震数据子集的波形反演思路、方法与应用

地震数据与地下介质物性参数之间的复杂关系,决定了地震全波形反演在理论方法上面临着强烈的非线性难题.地下不同物性参数的不同分量在地震数据上具有不同的表现,勘探的不同阶段对地下介质模型的精度也具有不同的要求,这就决定了在地震全波形反演过程中不必时刻追求地震数据全部信息的匹配,部分信息的匹配就有可能解决现阶段的某些问题,还可以一定程度上规避匹配全部地震信息所遇到的强烈非线性难题.基于这样的考虑,我们提出了利用地震数据子集进行波形反演的思路,给出了统一的反演方法,并通过基于包络数据子集以及反射波数据子集的波形反演的理论模型与实际资料反演试验,证明了所提出的波形反演思路和方法的正确性.