随钻方位电磁波测井多参数快速反演

大斜度井/水平井随钻测井技术已被广泛的应用于各种复杂油气藏中,但传统随钻电磁波测井仅提供仪器轴向磁场分量,井斜和各向异性对测井响应影响耦合在一起,难以剥离,无法基于反演同时获取多个参数,进行相应的地层评价,且不能提供方位信息,难以实现地质导向功能.随钻方位电磁波测井在提供仪器轴向磁场分量的基础上,还提供了仪器横向磁场分量,其提取的地质导向信号,对地层界面非常敏感,受井斜、地层电阻率对比度和各向异性影响弱,可以反演得到准确的界面位置.此外,额外的磁场分量信息,有效降低了井斜与各向异性的耦合程度,可基于反演算法实现多参数联合反演,满足电阻率测井精细解释评价需求.

     

基于方位各向异性的裂缝密度反演方法及应用

裂缝预测是页岩气甜点预测的重要内容.页岩气储层高角度裂缝发育,具有明显的HTI介质各向异性特征,基于各向异性参数与裂缝密度的相关性,提出了一种裂缝密度反演的新方法.首先推导了一种能反映裂缝密度的方位AVO近似公式及其方位弹性阻抗方程（EVAAz方程）;然后应用贝叶斯理论,形成了组稀疏方位弹性阻抗反演目标泛函,在此基础上通过阻尼最小二乘算法求解弹性矩阵,稳健得到了剪切模量、法向弱度与切向弱度;最后应用裂缝弱度与裂缝密度的关系,实现了裂缝密度的定量预测.经实际数据试验,预测结果与页岩压裂参数吻合度高,证实该方法是有效的.

     

Linearized rays,amplitude and inversion

In this paper, ray theoretical amplitudes and travel times are calculated in slightly perturbed velocity models using perturbation analysis. Also, test inversions using travel time and amplitude are computed. The pertubation method is tested using a 3-D velocity model for NORSAR having velocity variations up to 8.0 percent. The perturbed amplitudes are found to be in excellent agreement with the calculated ray amplitudes. Velocity inversions based on travel time and amplitude are next investigated. Perturbation analysis using linearized ray equations is efficiently used to compute amplitude derivatives with respect to model parameters. The trial linearized inversions use smaller velocity variations of 1.7 percent to avoid possible effects due to ray shift, even though the perturbation analysis is valid for larger variations. The trial 2-D inversion results show that linearized amplitude inversions are complementary and not redundant to travel time inversions, even in smoothly varying models.     

单斜介质中方位NMO速度Thomsen参数反演方法研究

采用启发式共轭梯度法,即随机爬山法 + 共轭梯度法,利用单斜介质中P波方位NMO速度椭圆轴向偏转角度接近于零这一特性,简化P波方位NMO速度公式,并利用多方位P波NMO速度,反演出某一初始CMP观测线与自然坐标系之间的夹角,作为进一步进行Thomsen各向异性参数反演的基础. 根据各向异性介质中方位NMO速度与Thomsen参数之间的关系,建立了利用三种波的多方位NMO速度及垂直传播速度反演单层单斜各向异性介质Thomsen各向异性参数的目标函数. 对计算的理论值添加具有一定标准差的正态分布的随机噪声,用以模拟实际观测存在的误差,通过对加噪后的数据进行多次反演的误差分析,表明了所建立的目标函数及选用的反演方法是有效可行的,而且相对稳定.     

Design spectra including effect of rupture directivity in near-fault region

In order to propose a seismic design spectrum that includes the effect of rupture directivity in the near-fault region, this study investigates the application of equivalent pulses to the parameter attenuation relationships developed for near-fault, forward-directivity motions. Near-fault ground motions are represented by equivalent pulses with different waveforms defined by a small number of parameters （peak acceleration, A, and velocity V; and pulse period, Tv）. Dimensionless ratios between these parameters （e.g., ATv/V, VTv/D） and response spectral shapes and amplitudes are examined for different pulses to gain insight on their dependence on basic pulse waveforms. Ratios of ATv/V, VTv/D, and the ratio of pulse period to the period for peak spectral velocity （Tv-p） are utilized to quantify the difference between rock and soil sites for near-fault forward-directivity ground motions. The ATv/Vratio of recorded near-fault motions is substantially larger for rock sites than that for soil sites, while Tvp/Tv ratios are smaller at rock sites than at soil sites. Furthermore, using simple pulses and available predictive relationships for the pulse parameters, a preliminary model for the design acceleration response spectra for the near-fault region that includes the dependence on magnitude, rupture distance, and local site conditions are developed.     

数据空间磁异常模量三维反演

强剩磁的存在通常导致了总磁化强度方向未知,进而影响了磁异常的反演和解释.磁异常模量是一种受磁化方向影响小的转换量,可以在强剩磁条件下通过反演三维磁化强度大小分布来推测场源分布状态.我们提出了一种数据空间磁异常模量反演算法来减少剩磁的影响.与标准的模型空间L2范数正则化反演方法相比,我们的方法有两个优点:一是无需搜索正则化参数(需要反复求解非线性反演问题),因而可以减少计算时间;二是反演结果更加聚焦,深度分辨率更高,我们对此进行了原因分析.通过模型和实测数据测试证明了该算法的有效性和更好的反演效果.     

Effects of an anisotropic overburden on azimuthal amplitude analysis in horizontal transverse isotropic media

Azimuthal amplitude variation in fractured media, commonly used to characterize fracture systems, is a function not only of reflection at the target but also of transmission through the overburden. This study investigates the sensitivity of amplitudes to various anisotropic overburden effects in horizontal transverse isotropic (HTI) media. Issues considered here are geometric spreading, transmission coefficients and attenuation due to fluid flow. Their influence on the azimuthal amplitude variation is evaluated quantitatively over a wide model space.
Only the variation of transmission coefficients with azimuth proves to be negligible. Geometric spreading alters the amplitude signature significantly over a relatively narrow range of models, and its influence increases with layer thickness. The most severe effect of an anisotropic fractured overburden is attenuation due to fluid flow between the cracks or the cracks and pores in the surrounding matrix. The relative changes in amplitudes between the symmetry directions due to anisotropic absorption are of the same order of magnitude as the changes in the reflection coefficient. The effect is significant over a very wide range of petrophysical parameters. Thus it leads to considerable problems in the amplitude analysis for almost any case of an overburden that contains cracks and pores. A correct amplitude analysis at the target will not be possible unless the effect of attenuation is removed.     

大斜度井/水平井随钻方位电磁波测井资料实时反演方法

水平井钻进过程中,地层边界的准确预测对地质导向决策至关重要.但水平井测井环境复杂,随钻方位电磁波测井响应不直观,储层性质、地层界面等信息的精确、定量评价困难.本文提出了一种适用于大斜度井/水平井随钻方位电磁波测井资料的电阻率及地层界面信息的实时提取方法.该方法主要基于以下几个技术：（1）引入滑动开窗策略将高维随钻方位电磁波测井资料的反演问题转化为多个连续窗口的一维反演问题;（2）对反演参数进行多初值初始化,结合正则化Levenberg-Marquardt算法保证解的全局最优性;（3）充分利用随钻方位电磁波测井多种探测深度信息,采用分级反演策略提高反演速度.数值和实测资料反演结果表明：降维反演方法可实现对地层界面的实时、准确提取,且适用于多数复杂地层结构和任意井眼轨迹;利用浅探测地质信号可实现随钻方位电磁波测井资料的分级反演,最终反演速度可提升10倍以上.

     

The mechanism of small mining tremors from amplitude inversion

A method of inverting seismic wave amplitudes from mining tremors is proposed to reveal their mechanism. A simple shear-implosive model of a mechanism is proposed. The shapes of theoretical and observed radiation patterns are compared in the least squares norm and the parameters describing the shear-slip together with the ratio of implosive and shear dislocations are optimized. The efficiency of the method, even in the case of a small number of stations, is demonstrated by applying it to weak mining tremors from two mine regions in Czechoslovakia.The weak mining tremors are found to be of a predominantly shear mechanism with a clear coincidence of nodal planes to local tectonics, which could imply their tectonic nature.     

井地磁异常模量联合反演

三维反演是磁测数据定量解释的重要方法,在金属矿勘探中扮演着重要的角色.但是在实际矿区的应用中,传统的磁总场异常反演方法依然存在两个问题：一是地面磁异常反演的深度分辨率较低,深部场源体的成像效果差;二是金属矿中可能包含强剩磁,反演结果可能是完全错误的.尽管前人对上述两个问题分别进行了广泛的研究,但尚未尝试同时解决这两个问题.本文在前人研究的基础上,提出了一种井地磁异常模量联合反演方法,该方法需要的控制参数少,无需加入额外的地质信息,且可用于多场源复杂磁异常的反演,具有较强的适用性.本文方法首先将地面和井中磁异常转化为模量数据,然后利用基于核函数或距离的加权函数将井地模量数据结合起来,使得该方法适用于联合反演.我们利用井地多种异常参量进行反演的模型试验表明,在强剩磁存在时,本文方法的效果优于其他方法,在减少剩磁影响的同时,也改善了深部成像效果,具有良好的应用前景.

     

Exact computation of input‐energy spectra from Fourier amplitude spectra

A theoretically exact relation is presented between the Fourier amplitude spectrum and the elastic input‐energy spectrum. This energy spectrum is computed from the Fourier acceleration amplitude spectrum and the real part of the relative velocity transfer function of the single‐degree‐of‐freedom elastic oscillator. Some accelerometric records collected worldwide are used to prove the solution. It is shown that the relation presented in this work opens the door for a better understanding of the relationship between the input‐energy spectrum and other seismological variables, such as magnitude and focal distance. Copyright © 2003 John Wiley & Sons, Ltd.     

P-wave amplitude spectra of Nevada underground nuclear explosions

Summary Short-period vertical-componentP-wave spectra have been determined for twelve Nevada underground nuclear explosions recorded by the Swedish seismic station network. Selected events cover the magnitude range fromm _b =5.9 tom _b =7.0 and the shot-depth range from 460 m to 1400 m. All explosion spectra show pronounced minima near 1 sec period. Within individual test areas the period at which the minimum occurs increases systematically with increasing shot-depth. This dependence makes it possible to interpret the observed minima in terms of interference betweenP- andpP-waves. One atmospheric explosion from Novaya Zemlya is analyzed to verify the inferred minima intrpretation.The minimum-period shows also a systematic linear dependence upon the magnitude. However, the physical explanation for this dependence is not evident and it may just be a consequence of a magnitude-depth relation.On leave from the Geophysical Institute, Czechosl. Acad. Sci., Prague 4-Spoilov, Boni II.     

Response spectra for nuclear structures on rock sites considering the near-fault directivity effect

Near-fault ground motions, potentially with large amplitude and typical velocity pulses, may significantly impact the performance of a wide range of structures. The current study is aimed at evaluating the safety implications of the near-fault effect on nuclear power plant facilities designed according to the Chinese code. To this end, a set of nearfault ground motions at rock sites with typical forward-directivity effect is examined with special emphasis on several key parameters and response spectra. Spectral comparison of the selected records with the Chinese and other code design spectra was conducted. The bi-normalized response spectra in terms of different corner periods are utilized to derive nuclear design spectra. It is concluded that nuclear design spectra on rock sites derived from typical rupture directivity records are significantly influenced both by the earthquake magnitude and the rupture distance. The nuclear design spectra specified in the code needs to be adjusted to reflect the near-fault directivity effect of large earthquakes.     

Comparison of waveform inversion, part 3: amplitude approach

In the second paper of this three part series, we studied the case of conventional and logarithmic phase‐only approaches to full‐waveform inversion. Here, we concentrate on deriving amplitude‐only approaches for both conventional‐ and logarithmic‐based methods. We define two amplitude‐only objective functions by simply assuming that the phase of the modelled wavefield is equal to that of the observed wavefield. We do this for both the conventional least‐squares approach and the logarithmic approach of Shin and Min. We show that these functions can be optimized using the same reverse‐time propagation algorithm of the full conventional methodology. Although the residuals in this case are not really residual wavefields, they can both be considered and utilized in that sense. In contrast to the case for our phase‐only algorithms, we show through numerical tests that the conventional amplitude‐only inversion is better than the logarithmic method.     

Long period Fourier amplitude spectra of strong motion acceleration

The empirical equations for scaling Fourier amplitude spectra in the frequency band from ~0·1 to 25 Hz can be extrapolated to describe the long period strong motion amplitudes. The results of this extrapolation can agree with (1) the seismological and field estimates of permanent ground displacement (near field), and with (2) the independent estimates of seismic moment and the observed frequencies of far field Fourier spectrum amplitudes.     

Rapid least-squares inversion of apparent resistivity pseudosections by a quasi-Newton method1

A fast inversion technique for the interpretation of data from resistivity tomography surveys has been developed for operation on a microcomputer. This technique is based on the smoothness-constrained least-squares method and it produces a two-dimensional subsurface model from the apparent resistivity pseudosection. In the first iteration, a homogeneous earth model is used as the starting model for which the apparent resistivity partial derivative values can be calculated analytically. For subsequent iterations, a quasi-Newton method is used to estimate the partial derivatives which reduces the computer time and memory space required by about eight and twelve times, respectively, compared to the conventional least-squares method. Tests with a variety of computer models and data from field surveys show that this technique is insensitive to random noise and converges rapidly. This technique takes about one minute to invert a single data set on an 80486DX microcomputer.     

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

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

     

转换波方位各向异性裂缝检测技术研究及应用

HTI裂缝各向异性介质中,转换波随方位角的变化比较复杂,目前还没有解析公式可以表达其变化特征,只能通过物理实验或数值模拟来分析其应用的可能性.数值模拟结果表明,转换波在裂缝各向异性介质中传播时,其R分量和T分量的振幅属性都具有明显的方位各向异性特征,R分量振幅方位各向异性拟合椭圆的长轴方向指示裂缝方位,这与纵波方位各向异性特征相似;根据P波AVAZ方位各向异性分析原理,对转换波R分量振幅方位各向异性曲线进行方位椭圆拟合,寻找椭圆的长轴方向,即裂缝主方位,再由P波AVAZ技术中振幅响应与炮检方向和裂缝走向之间的夹角关系式得到裂缝的发育密度,从而构建转换波方位各向异性AVAZ裂缝检测技术.该技术已用于川西新场气田某区块的裂缝储层预测,取得了较好的应用效果.     

Seismic moment tensor of Pribaikalye earthquakes from the surface-wave amplitude spectra

The method of surface-wave amplitude spectra inversion for the seismic moment tensor (SMT) is implemented and tested in the Pribaikalye region. The SMTs are calculated for 39 events with M _w = 4.4–6.3, which occurred in the region in 2000–2011. Based on the obtained data, the seismotectionic deformations of the crust are estimated in two seismically active areas-the Northern Pribaikalye and northeastern Baikal rift zone. It is found that on a level of moderate-magnitude events, the region is dominated by the regimes of subhorizontal northwestern extension and strike-slip faulting, which reflects the long-term trends in the stress field of the crust in these parts of the rift.