Strain-threshold- and frequency-dependent seismic simulation of nonlinear soils

A one-dimensional equivalent linear method (EQL) is widely used in estimating seismic ground response. For this method, the shear modulus and damping ratio of inelastic soil are supposed to be frequency independent. However, historical earthquake records and laboratory test results indicate that nonlinear soil behavior is frequency-dependent. Several frequency-dependent equivalent linear methods (FDEQL) related to the Fourier amplitude of shear strain time history have been developed to take into account the frequency-dependent soil behavior. Furthermore, the shear strain threshold plays an important role in soil behavior. For shear strains below the elastic shear strain threshold, soil behaves essentially as a linear elastic material. To consider the effect of elastic-shear-strain-threshold- and frequency-dependent soil behavior on wave propagation, the shear-strain-threshold- and frequency-dependent equivalent linear method (TFDEQL) is proposed. A series of analyses is implemented for EQL, FDEQL, and TFDEQL methods. Results show that elastic-shear-strain-threshold- and frequency-dependent soil behavior plays a great influence on the computed site response, especially for the high-frequency band. Also, the effect of elastic-strain-threshold- and frequency-dependent soil behavior on the site response is analyzed from relatively weak to strong input motion, and results show that the effect is more pronounced as input motion goes from weak to strong.     

准噶尔盆地腹部隐蔽油气藏物理模拟及地震响应特征

准噶尔盆地腹部侏罗系以隐蔽油气藏为主,地震资料在目的层分辨率较低、地震反射特征不清晰,制约了对该区的低幅构造、岩性圈闭的识别和评价.本文根据侏罗系的沉积特点,建立河流、（扇）三角洲沉积体系地质模型,制作尖灭体、透镜体、河道砂、浊积体等地质体的三维物理模型,通过实验测试方法,研究地层条件下岩石物理性质的变化特征,系统研究储层地球物理参数变化的总体规律,并结合测井资料,采用波场物理模拟方法,研究波场响应特征与储层参数变化之间的规律,分析各种情况下的波场特征、属性特征.     

团树断层场地地震反应的数值模拟

团树断层场地在多次地震中出现了烈度异常。为了弄清烈度异常的原因，应用显式有限元方法，数值模拟了在脉冲地震波斜入射情形下断层场地的地震反应。结果表明：在基岩山脚附近土层场地地震动放大比基岩山区平均高4倍、比远离基岩山脚土层场地平均高2倍，与历史震害分布有较好的对应。     

可控震源非线性扫描地震响应的数值模拟

讨论了几种非线性扫描信号自相关函数的旁瓣特性.对可控震源非线性扫描的地震响应进行了数值模拟.当震源扫描持续时间是检波器接收时间的一半时,相关运算把反射扫描信号压缩成脉冲信号的效果是显著的.有限度地增加震源扫描持续时间以及与之相适应的接收时间可以使反射扫描信号得到进一步的压缩,但在固定扫描与接收时间的条件下,增加采样点个数对反射扫描信号的压缩并不起作用.采用具有低频相关子波特性的非线性扫描震源信号,将有利于相关噪声的消除.模拟结果还证实了二次扫描在某些特定条件下会出现鞍点效应这样一个事实.     

地震资料中的多尺度断裂与流体响应特征：数值模拟分析

流体在断裂和岩石骨架间的交换被认为是影响岩石弹性参数各向异性的主要原因,理论研究表明断裂尺度同样对弹性参数的各向异性也有影响.为了说明两者对各向异性影响以实现多尺度断裂裂隙的识别,本文在等效介质模型的基础上,应用数值分析的方法研究速度和衰减(1/Q)随多尺度断裂、频率和流体因子变化规律.结果表明介质弹性参数是频率依赖的,并且参数中存在衰减项,而这种频率依赖性与介质物性参数中的断裂尺度及流体性质存在一定的联系;当断裂定向分布时,参数结果显示为各向异性;不同断裂尺度具有不同的波速频散特性,剪切波分裂程度依赖于频率,断裂尺度起着控制作用,高频时对小尺度的敏感,低频段对大尺度敏感.在地震频段Thomsen参数随着频率的增大而减小,随着断裂尺寸的增大而减小.因此地震数据可能区分断裂和微裂隙引起各向异性,从而可测量断裂尺度.     

Numerical simulation of catastrophic flood: the case study of hypothetical failure of the Bielkowo hydro-power plant reservoir

The numerical modeling of flood wave propagation following the hypothetical breaks of the embankments of the Bielkowo hydro-power plant storage reservoir (Kolbudy II Reservoir) on the Radunia River in Poland has been presented. The results of computations were used to estimate the parameters of the flood waves, which are indispensable for the flood zone determination and mapping and then for the flood risk analysis. When estimating the reach and area of the inundation, related to the embankments failures, digital terrain model, and mathematical model of flood wave propagation are necessary. For the numerical simulations of flood, the mathematical model of free surface, two-dimensional unsteady water flow was applied. Four locations of potential breaks of the reservoir embankments were considered. The computed flood zones were presented on the flood hazard maps. The maps have been used by the local authorities and the dam owner to manage the flood risk related to hydro-power plants operations on the Radunia River. This type of research has been done for the first time for the water plant managed by the ENERGA Elektrownie Straszyn.     

浮放设备地震滑移反应数值模拟研究

以前人研究建立的浮放设备滑移反应分析理论为基础,给出了浮放设备在水平和竖向地震输入下滑移反应运动微分方程,选取连续型的库仑摩擦力模型,采用Runge-Kutta法求解浮放设备地震滑移反应运动微分方程,可以得到浮放设备在地震作用下的绝对加速度、相对速度和相对位移反应时程。采用编制的计算程序,进行了浮放设备地震滑移反应参数影响研究,结果表明:浮放设备水平滑移反应随着水平地震地面输入或楼层反应输入的增大而增大,随着设备与支撑面之间摩擦系数的增大而减小;当水平向输入加速度峰值大于0.3g时,需要考虑竖向激励的影响。     

Monitoring viscosity changes from time‐lapse seismic attenuation: case study from a heavy oil reservoir

Heating heavy oil reservoirs is a common method for reducing the high viscosity of heavy oil and thus increasing the recovery factor. Monitoring of these viscosity changes in the reservoir is essential for delineating the heated region and controlling production. In this study, we present an approach for estimating viscosity changes in a heavy oil reservoir. The approach consists of three steps: measuring seismic wave attenuation between reflections from above and below the reservoir, constructing time‐lapse Q and Q^?1 factor maps, and interpreting these maps using Kelvin–Voigt and Maxwell viscoelastic models. We use a 4D relative spectrum method to measure changes in attenuation. The method is tested with synthetic seismic data that are noise free and data with additive Gaussian noise to show the robustness and the accuracy of the estimates of the Q‐factor. The results of the application of the method to a field data set exhibit alignment of high attenuation zones along the steam‐injection wells, and indicate that temperature dependent viscosity changes in the heavy oil reservoir can be explained by the Kelvin–Voigt model.     

河流相储层特征及识别:一个老河口油田的实例分析(英文)

河流相储层是一种重要的储层类型,寻找河道砂体是油气勘探的一项重要工作。但是河道变迁比较频繁,叠置比较严重,加之地震资料本身信噪比和分辨率的限制,河流相储层的精细描述至今仍存在一定的困难。本文以胜利油田老河口三维工区为例,首先研究了本区河流相储层的基本特征,发现该类储层剖面上呈强相位、切片上连续性好和频谱上呈低频特征;通过颜色处理和地震精细解释,发现本区河流相储层存在类串珠状特征,为了验证这一观测结果,我们用正演模拟分析了其形成机理,为在类似工区识别此类古河道找到了新的途径;应用多属性融合和RGB显示技术,河道下切特征更加明显,很好地彰显了河道结构特征,提高了河道识别能力;最后,我们研究并开发了多子波检测技术,检测出更多河流相弱反射信息。     

Understanding the fracture role on hydrocarbon accumulation and distribution using seismic data: A case study on a carbonate reservoir from Iran

The South Pars, the largest gas field in the world, is located in the Persian Gulf. Structurally, the field is part of the Qatar–South Pars arch which is a regional anticline considered as a basement-cored structure with long lasting passive folding induced by salt withdrawal. The gas-bearing reservoir belongs to Kangan and Dalan formations dominated by carbonate rocks. The fracture role is still unknown in gas accumulation and distribution in this reservoir. In this paper, the Scattering Index (SI) and the semblance methods based on scattered waves and diffraction signal studies, respectively, were used to delineate the fracture locations. To find the relation between fractures and gas distribution, desired facies containing the gas, were defined and predicted using a method based on Bayesian facies estimation. The analysis and combination of these results suggest that preference of fractures and/or fractured zones are negligible (about 1% of the total volume studied in this paper) and, therefore, it is hard to conceive that they play an important role in this reservoir. Moreover, fractures have no considerable role in gas distribution (less than 30%). It can be concluded from this study that sedimentary processes such as digenetic, primary porosities and secondary porosities are responsible for the gas accumulation and distribution in this reservoir.     

近断层地震下可液化场地土反应的数值模拟

液化条件下场地土大变形是造成工程结构失效的主要原因之一。文中以某高速公路特大桥的可液化岸坡近场场地为研究对象,基于PL-Finn液化本构模型,利用FLAC 3D程序对其在地震作用下的动响应全过程进行了数值模拟分析。结果表明,PL-Finn模型可较好地反映地震过程中孔隙水压、液化区的变化规律,并能较好的预测液化后场地土的变形规律。液化引起的地基土大变形对桥梁桩基的影响需引起重视。     

唐山地震人员伤亡数值模拟与评估

应用作者提出的地震人员伤亡指数和状态函数的概念,进一步给出初始人员伤亡矩阵,并借助于几个概率分布函数确定了砖砌体结构的初始人员伤亡矩阵,建立了地震人员伤亡动态量化评估方法,以此方法对唐山地震进行了人员伤亡数量的模拟与评估。     

单体边坡地形的地震动力响应及其放大效应的数值分析

本文采用隐式动力有限单元法研究了不同的边坡角度和边坡高度对地形放大效应的影响,并以位移峰值放大系数为衡量地震动放大效应的标准,计算了不同边坡角度和边坡高度条件下的地震响应,在此基础上对模型关键监测点的输出波形以及位移峰值放大系数的变化趋势进行了分析,获得了不同监测点处的地震动时程曲线,揭示了坡角和坡高对单体边坡地震动放大效应的定量作用规律。数值结果表明,相同高度处坡面监测点的水平向位移峰值放大系数大于坡内监测点的,地形放大效应在水平方向具有趋表效应。由于坡面存在入射波和反射波的叠加,因此竖直向位移峰值放大系数的最大值出现在坡体内部。      

叠前反演和地震吸收技术在复杂天然气藏地震预测中的应用

针对某复杂断块天然气目标储层,在岩石物理分析的指导下,综合利用地质、地震、测井等资料,提出了一套面向复杂天然气藏的叠前地震预测技术.首先基于地震岩石物理分析得到的初始横波信息,采用叠前贝叶斯非线性三参数反演得到了井旁控制点处精确纵横波速度和密度信息,然后通过叠前/叠后联合反演技术实现了面向目标的弹性阻抗体反演及含气储层敏感参数直接提取,最后结合小波变换时频谱分析的方法从叠前地震资料中估算地层吸收参数值,提高天然气藏识别精度.实际应用表明,综合各种叠前地震预测技术,可以大大提高对复杂天然气藏的识别精度,降低勘探风险.     

Numerical simulation of reservoir sediment and effects on hydro‐dynamic response of arch dams

Based on the dynamic theory for saturated porous media by Biot (Journal of the Acoustical Society of America 1956; 28 : 168–178), a numerical model is presented to analyse the reflection behaviours of reservoir sediment and compared with those from the visco‐elastic model. It is concluded that the two models give very similar results of reflection coefficient α within the frequency range of interest. Then, using the two models, the change of the reflection coefficients α with various sedimentation parameters and excitation frequencies are studied in detail. The results are further used in the analysis of response functions of hydro‐dynamic pressures on, and structural displacements of the Xiang Hong Dian arch dam, for which some results from a field vibration test are available. It appears that effects of water compressibility with sediment reflection on hydro‐dynamic pressures and structural response are not significant for this specific case. Copyright © 2001 John Wiley & Sons, Ltd.     

地震散射波模拟成像方法在铜陵某矿区的应用

针对深部隐伏的金属矿床,以铜陵某矿区为例,应用散射波正演模拟技术,对金属矿地震数据野外采集中涉及的偏移距、覆盖次数、道间距等参数进行了不同炮记录的模拟对比,辅助指导非均匀地质体金属矿地震勘探野外数据采集的观测系统设计;针对该地区实际采集的地震数据,应用基于等效偏移距的共散射点成像方法,对金属矿地震数据进行了处理,通过与常规反射波法处理结果对比,散射成像效果较好,说明散射波技术在金属矿地震勘探中应用有一定的可行性和适用性.     

地震低频伴影的数值模拟与应用

Strong low-frequency energy beneath a hydrocarbon reservoir is called a seismic low-frequency shadow and can be used as a hydrocarbon indicator （Tarter et al., 1979） bu the physical mechanism of the observed low-frequency shadow is still unclear. To stud） the mechanism, we performed seismic numerical simulation of geological models with a hydrocarbon-bearing zone using the 2-D diffusive-viscous wave equation which car effectively model the characteristics of velocity dispersion and transform the seismic dat~ centered in a target layer slice within a time window to the time-frequency domain by usinl time-frequency signal analysis and sort the frequency gathers to common frequency cubes. Then, we observe the characteristics of the seismic low-frequency shadow in the common frequency cubes. The numerical simulations reveal that the main mechanism of seismic lowfrequency shadows is attributed to high attenuation of the medium to high seismic frequency components caused by absorption in the hydrocarbon-filled reservoir. Results from a practical example of seismic low-frequency shadows show that it is possible to identify the reservoir by the low-frequency shadow with high S/N seismic data.     

2.5维非均匀介质中的地震波数值模拟

目前,完全的3D弹性波数值模拟计算仍然需要庞大的计算资源,而2D弹性波数值模拟,又不能准确地近似3D数值模拟.因此,本文在较高数值精度的一阶应力一速度弹性波动方程的基础上,采用时间错格差分拟谱方法在2D模型介质中计算3D弹性波场,即2.5D数值模拟.并将3D计算结果与2.5D和2D结果进行比较,发现2.5D模拟结果在波场形态和振幅方面比2D模拟结果更接近3D模拟结果,从而证实了2.5D数值模拟比2D数值模拟可以更好地近似3D波场数值模拟,同时证实了时间错格差分拟谱方法是一种高精度、高效率的正演模拟方法.