流体饱和周期性层状孔隙介质模型平行层面传播的纵波频散和衰减研究

孔隙介质中的地震波传播一直是油气地震勘探领域的研究热点和难点问题.该科学难题源自不同尺度的裂隙、孔隙、溶洞与岩石骨架之间的耦合作用,导致地震波场特征复杂.目前相关的研究主要集中于探索孔隙介质中地震波的传播机制及地震响应的特征与变化规律,包括对地震波在复杂孔隙介质中传播,进行比较精确的数学物理描述以及数值实现.地球物理学家们集中于研究垂直于地层层面方向入射的地震波频散和衰减,而忽略了实际地球介质中的地震波是以任意角度(方向)入射并进行传播的普遍性情况.在前人的研究基础上,本文的创新之处在于将纵波的入射方向扩展到平行于流体饱和的周期性层状孔隙介质模型层面方向.针对流体饱和的周期性层状孔隙介质模型,提出了介观波致流(Wave-induced Fluid Flow, WIFF)对流体饱和孔隙层状介质中平行于层面方向入射的纵波频散、衰减及频变各向异性的新模型.利用准静态Biot孔弹性方程推导出了模型的孔隙压力、流体流动速度、平均应力和平均应变等物理量的解析表达式,进而得到流体饱和的周期性层状孔隙介质复纵波模量的精确解析解．然后,利用复纵波模量讨论了纵波速度频散、衰减和频变各向异性特征,讨论了背景...     

基于Biot-Squirt方程的波场模拟

Biot流动和喷射流动是含流体多孔隙介质中流体流动的两种重要力学机制,对地震波和声波的传播均产生重要影响. Dvorkin和Nur提出了同时包含Biot流动和喷射流动力学机制的统一的BISQ(Biot-Squirt)模型,基于这一模型,尽管有关弹性波在多孔隙介质中的衰减和频散问题已被广泛研究,然而,基于BISQ波传播方程的波场数值模拟至今仍未见报道. 本文从同时包含两种力学机制的孔隙弹性波方程出发,利用FCT有限差分法对含流体孔隙各向同性介质中的地震波和声波进行了数值模拟,并与基于Biot流动的Biot理论之模拟结果进行比较. 数值模拟结果表明:同时包含Biot流动和喷射流动影响的地震波和声波速度比仅包含Biot流动作用的地震波和声波速度慢,慢P波的衰减比根据Biot理论模拟的慢P波衰减更强.     

裂缝诱导双相HTI介质地震波场错格伪谱法模拟与波场特征分析

裂缝诱导的双相具有水平对称轴的横向各向同性(HTI)介质模型是由一组平行排列的垂直裂缝嵌入到统计各向同性的流体饱和多孔隙岩石中而组成的,它综合考虑了裂缝型储层岩石的各向异性和孔隙性.高精度的地震波场数值模拟技术是研究该介质中地震波传播规律的主要方法.本文结合错格伪谱法和时间分裂法,求解描述该介质中地震波传播的一阶速度-应力方程.模拟了单层和双层模型中的地震波场,并对其进行了特征分析.研究结果表明:错格伪谱法能有效消除标准网格伪谱法波场模拟结果中出现的数值伪影现象,与时间分裂法结合能够获得稳定的、高精度的模拟结果;裂缝诱导双相HTI介质中的地震波场兼具裂缝各向异性介质和双相介质中传播的地震波的波场特征.     

随机弹性介质中地震波散射衰减分析(英文)

地震波衰减一直是许多学科研究的热点,因为可以反映介质的特性。导致地震波衰减的因素很多,如:传播过程中由于能量扩散导致的几何衰减,固体岩石内部晶粒间相对滑移导致的摩擦衰减,岩石结构不均匀引起的地震波散射衰减。本文主要从统计的观点出发,通过多次数值模拟的方法研究纵波散射在随机弹性介质中所引发的衰减。首先用随机理论建立了二维空间随机弹性介质模型,然后用错格伪谱法的数值方法模拟了波在随机介质中的传播,再通过波场中虚拟检波器的记录,用谱比法估计了弹性波在随机介质中的散射衰减。不同非均匀程度随机弹性介质中的数值结果表明:介质不均匀程度越高,散射衰减越大;在散射体尺寸小于波长的前提下,不同散射体尺寸的计算结果说明:散射体尺寸越大,弹性波衰减越明显。最后提出了一种不均匀孔隙介质中流体流动衰减的方法。通过对随机孔隙介质中地震波的总衰减和散射衰减分别进行了计算,并定量得出了随机孔隙介质中流体流动衰减,结果表明:在实际地震频段下,当介质不均匀尺度101米量级时,散射衰减比流体流动衰减要大,散射衰减是地震波在实际不均匀岩石孔隙介质中衰减的主要原因。     

方位各向异性粘弹性介质波场数值模拟

当地震信号通过复杂地球介质时，地层除了表现为各向异性，还表现为内在的粘弹性特征.因此，为准确描述地震波在地球介质中的传播特征，理想的地球介质模型应该能够模拟岩石的各向异性特征和衰减特征.本文给出了各向异性粘弹性介质模型的波动方程及其差分格式，并利用有限差分法实现了地震波波场数值模拟.结果表明了该介质模型中地震波场特征与各向异性主轴方位和介质的粘滞性参数之间的关系.     

各向同性介质中射线理论与有限差分地震波场模拟方法比较研究

地震波场数值模拟不仅是研究复杂地区地震资料采集、处理和解释的有效辅助手段,而且是研究地球深部精细构造和地球深部探测的有效工具.射线理论和波动方程理论是地震波场数值模拟的理论基础.射线理论主要刻画地震波在介质中传播的走时场、地震射线等运动学属性;波动方程理论通过求解波动方程来描述地震波在介质中传播的弹性动力学响应(能量衰减、相位特征、偏振属性、以及全波形等).基于波动方程理论的波场数值模拟由于能够引入丰富的波场信息,使得人们对不同介质中地震波的传播过程有了较全面的了解.本文以二维层状均匀介质模型为例,通过射线追踪法和交错网格有限差分法模拟得到的波场快照图、单炮地震记录剖面、合成理论地震图的分析比较,不但对地震波在各向同性层状均匀介质中的传播规律和特点有了深刻的理解和认识,同时又可以相互验证两种不同方法的正确性和有效性.     

含夹层孔隙介质分界面上反透射系数特征研究

研究地震波在含黏性流体孔隙介质中的传播特性,对于识别油气储层、刻画流体边界、判识油水界面等有重要作用.通常利用地震波在孔隙介质分界面上的反透射系数特征,来实现这些方面的研究.顶、底层孔隙介质之间存在含黏性流体孔隙介质夹层时,夹层内岩石骨架、孔隙因子、孔隙流体等参数会对反透射系数产生影响.本文推导出含夹层孔隙介质分界面上的反透射系数方程,针对方程建立不同孔隙度的含夹层孔隙介质模型;基于模型实现孔隙介质分界面上的反透射系数特征研究,分别分析不同夹层厚度、不同频率情况下,反透射系数随入射角的变化规律.研究表明,快慢纵波及横波的反透射系数都会受到孔隙夹层的影响,受影响程度主要取决于夹层厚度和地质背景等因素;慢纵波反射和透射系数数值极小,且后者随夹层厚度的增加而迅速衰减;频率通过反透射系数方程本身和地震波速度两个因素直接或间接的影响地震波反透射特征,其值的增加会引起快纵波和横波反透射系数不同程度的衰减.     

含混合裂隙、孔隙介质的纵波衰减规律研究

地下多孔介质中的孔隙类型复杂多样,既有硬孔又有扁平的软孔.针对复杂孔隙介质,假设多孔介质中同时含有球型硬孔和两种不同产状的裂隙(硬币型、尖灭型裂隙),当孔隙介质承载载荷时,考虑两种不同类型的裂隙对于孔隙流体压力的影响,建立起Biot理论框架下饱和流体情况含混合裂隙、孔隙介质的弹性波动方程,并进一步求取了饱和流体情况下仅由裂隙引起流体流动时的含混合裂隙、孔隙介质的体积模量和剪切模量,随后,在此基础上讨论了含混合裂隙、孔隙介质在封闭条件下地震波衰减和频散的高低频极限表达式.最后计算了给定模型的地震波衰减和频散,发现地震波衰减曲线呈现"多峰"现象,速度曲线为"多频段"频散.针对该模型分析讨论了渗透率参数、裂隙纵横比参数以及流体黏滞性参数对于地震波衰减和频散的影响,表明三个参数均为频率控制参数.     

频率域有限差分数值模拟地震波衰减和震电效应特征

地震波在地下含流体孔隙介质中传播时,会引起中观尺度的"局域流",进而产生地震波震电效应.基于Biot(1941)固结理论的准静态方程,在频率域中采用空间有限差分方法,正演模拟虚岩石物理岩样的地震波衰减和震电效应.与时间域虚岩石物理方法相比,该方法既可以直接求取任一频率下的地震波衰减和电势,便于应用于实际岩样的预测分析,也避免了讨论岩样外表面施加的力源函数表达式及时间剖分稳定性条件等问题.首先利用周期性层状介质模型验证了本文所描述方法的有效性,并进一步求取分析了周期性层状介质两种不同特征单元的渗流电流密度及电势,数值模拟结果表明由中观尺度"局域流"引起的震电效应电势振幅数量级在实验室测量范围之内,随后,分析研究了四种不同高渗介质占比值的地震衰减及震电效应特征.最后,将本文提出的震电效应数值计算方法推广至二维,并求取了二维斑块饱和模型的地震波衰减、速度频散、电势的振幅和相位角数值结果.     

双相介质中地震波的频率-空间域数值模拟

本文利用优化的25点频率-空间域有限差分算法对基于BISQ模型双相各向同性介质中的地震波进行了数值模拟.通过与经典的Biot模型理论模拟结果进行对比,分析了Biot流动（宏观流体流动）和Squirt流动（微观流体流动）耦合作用对地震波在孔隙介质中传播特性的影响.数值模拟在地震频段进行,结果显示：在理想相界和黏滞相界情况下,Squirt流动机制都比Biot流动机制产生了更大的速度频散和能量衰减.其中,在Biot流动和Squirt流动耦合作用下的快P波的速度和振幅小于仅考虑Biot流动影响下快P波速度和振幅,而且慢P波的衰减也更加强烈.本文还研究了地震波在双层双相各向同性介质分界面处的反射和透射特征,双相介质中波的反射与透射现象类似于单相介质的情况.模拟结果表明,利用优化25点频率-空间域有限差分法模拟双相孔隙介质中的地震波场是可行的,这为开展双相孔隙介质全波形反演问题的研究提供了可能.     

Pattern analysis of dune‐field parameters

As a basic form of pattern analysis, the parameters of dune spacing, defect density, crest orientation and crest length are measured from remote images and treated statistically for dunes at White Sands in New Mexico, the Algodones in California, the Agneitir in Mauritania, and the Namib in Namibia. Statistical populations are identified from frequency plots of dune spacing and crest length, field‐scale calculations of defect density, and rose diagrams of crest orientation. Single populations characterize simple dune fields (White Sands), whereas multiple populations characterize compound/complex dunes (Algodones, Namib), and complex dune fields (Agneitir). As time increases, dune fields show an increase in dune spacing and crest length, a decrease in defect density, more tightly clustered crest orientation, and a reduction in the variance associated with measurements of these parameters. The results are consistent with models of dune fields as self‐organizing complex systems in which a characteristic pattern emerges as a function of constructional time. Because pattern evolution is a function of time, it may be possible to use pattern analysis to augment current methods of age determination. Statistically defined populations can be used in geomorphic backstripping to unstack generations of simple patterns that give rise to complex patterns, and to reconstruct each generation in terms of construction time and palaeo‐wind regime. Copyright © 2006 John Wiley & Sons, Ltd.     

裂隙弱度参数的波场特征

本文将Schoenberg线性滑移裂隙介质理论应用到HTI(LSD)介质的弹性波波场数值模拟之中.主要针对描述裂隙状态的2个裂隙参数--法向弱度EN和切向弱度ET对波场的影响进行了有限元模拟.模拟结果揭示了EN主要影响快波而ET主要影响慢波的波场特征,通过EN和ET变化对应的裂隙含水程度说明了EN对裂隙中的流体敏感而E...     

On the parameters of the normal earth's gravity field

Summary The parameters of the normal Earth's gravity field, determined from contemporary satellite data are discussed and compared with the parameters of the normal Helmert system (1901 to 1909).Dedicated to Academician Alois Zátopek on His 65th Birthday     

利用多种地震学参数研究中国大陆地壳应变场

地壳应变场研究是地球动力学研究的一个基础性领域, 以往文献的探索大多集中在地壳形变的测量和研究上, 包括大面积水准测量和GPS观测等。 文中利用多种地震学参数研究中国大陆地壳应变场, 利用哈佛CMT目录和现代中国地震目录得到了中国大陆不同区域面波震级与标量地震矩的统计关系, 进而采用中国大陆历史地震资料研究中国大陆地壳应变场, 给出地震最大剪切应变率场的分布。 研究结果显示, 青藏高原及其周边地区是最大剪切应变率的高值区, 帕米尔和阿萨姆地区的应变率值最大。 与GPS得到的应变场结果进行比较, 两者的结果存在一致趋势。 利用NEIC宽频带地震辐射能量目录研究了中国大陆的地震视应变分布, 显示地震的视应变分布与地壳应变水平存在一定相关关系, 且发生在青藏高原周边地区地震的视应变水平较高。 这些研究结果为地球动力学研究的进一步深入探讨提供了新的科学资料。     

Research of pattern dynamics parameters of crustal deformation field in seismogenic process

ＲｅｓｅａｒｃｈｏｆｐａｔｔｅｒｎｄｙｎａｍｉｃｓｐａｒａｍｅｔｅｒｓｏｆｃｒｕｓｔａｌｄｅｆｏｒｍａｔｉｏｎｆｉｅｌｄｉｎｓｅｉｓｍｏｇｅｎｉｃｐｒｏｃｅｓｓＳｈｕｏ－ＹｕＺＨＯＵ（周硕愚），ＹｕｎＷＵ（吴云），Ｒｕｏ－ＢｏＷＡＮＧ（王若柏）ａｎｄ...     

Correlation between the near-Earth solar wind parameters and the source surface magnetic field

The solar magnetic field B _s at the Earth’s projection onto the solar-wind source surface has been calculated for each day over a long time interval (1976–2004). These data have been compared with the daily mean solar wind (SW) velocities and various components of the interplanetary magnetic field (IMF) near the Earth. The statistical analysis has revealed a rather close relationship between the solar-wind parameters near the Sun and near the Earth in the periods without significant sporadic solar and interplanetary disturbances. Empirical numerical models have been proposed for calculating the solar-wind velocity, IMF intensity, and IMF longitudinal and B _z components from the solar magnetic data. In all these models, the B _s value plays the main role. It is shown that, under quiet or weakly disturbed conditions, the variations in the geomagnetic activity index Ap can be forecasted for 3–5 days ahead on the basis of solar magnetic observations. Such a forecast proves to be more reliable than the forecasts based on the traditional methods.     

Magnetopause magnetic barrier parameters in dependence on the solar wind magnetic field orientation

Some theories predict the magnetosheath magnetic field strength will decrease and the density increase just outside the dayside magnetopause as the interplanetary magnetic field turns southward. Two studies have recently reported results which confirm these expectations. In contrast, we briefly review our own theoretical predictions which indicate that precisely the opposite effect is expected. We survey new and previously reported magnetosheath observations and demonstrate that they are consistent with the predictions of our model. The conflicting results indicate a need for further theoretical and observational work.     

Synchronization trends and rhythms of multifractal parameters of the field of low-frequency microseisms

The field of low-frequency microseisms is investigated with the use of data from 83 stations of the F-net broadband network in Japan over the period from the beginning of 1997 through June 2008. Vertical components with a sampling step of 1 s are used for analysis, as well as signals with a sampling step of 1 min obtained from the initial data by averaging and thinning. Long-period regularities of changes in the singularity spectrum support width Δα and the generalized Hurst exponent α* for the field of low-frequency microseisms were revealed by estimating multifractal singularity spectra in consecutive time windows 30 min long for 1-s data and 24 hour long for 1-min data. The average value of the parameter α* for 1-s data significantly decreased before the Hokkaido earthquake of September 25, 2003 (M = 8.3), and was not restored subsequently to its previous level. Prior to September 2003, 1-min data on α* variations experienced strong annual changes, which completely ceased afterwards. Both these effects are interpreted as an increase in the degree of synchronization of microseismic noise on Japan’s islands after the September 25, 2003, earthquake. This hypothesis is also supported by estimates of the measures of correlation and spectral coherence between variations in the average values of Δα and α* calculated for 1-min data inside five spatial clusters of stations from consecutive time fragments two months long. Based on the well-known statement of the theory of catastrophes that synchronization is one of the flags of an approaching catastrophe, it was suggested that the Hokkaido event could be a foreshock of an even stronger earthquake nucleating in the region of Japan’s islands.     

Identification of modal parameters of a 600-m-high skyscraper from field vibration tests

Accurate prediction of the dynamic responses of a high-rise building subjected to dynamic loads such as earthquake and wind excitations requires the information of its structural dynamic properties such as modal parameters including natural frequencies and damping ratios. This paper presents the identification results of the modal parameters based on field vibration tests on a 600-m high skyscraper. A set of tests, including ambient vibration test (AVT) and free vibration test (FVT), were conducted on the skyscraper to identify its modal parameters. Firstly, this paper presents and discusses the modal parameters of the skyscraper assessed by several identification methods applied to the AVT measurements. These methods include the wavelet transform (WT) method, the stochastic subspace identification (SSI) method, and the random decrement technique (RDT). Secondly, an active mass damper (AMD) system with total mass 1000 tons equipped into the skyscraper was used to excite the building for estimation of the modal parameters by FVT. Thirdly, this paper presents observations on the structural dynamic behavior of the skyscraper with the operation of the AMD system during a typhoon event. The field measurement results show that the AMD system functioned efficiently for suppression of the wind-induced vibrations of the skyscraper during the typhoon. This paper aims to further understand the structural dynamic properties of super-tall buildings and provide useful information for structural design and vibration control of future skyscrapers.     

Influence of solar wind parameters on ground magnetic field near equator during intense storms

Using a simple time-lagged correlation technique, present study aims to identify the solar wind (SW) parameter, which is better associated with the ground magnetic field variations of shorter time duration near equator, during intense geomagnetic storms. It is found that out of all SW parameters, successively occurring enhancements in the SW dynamic pressure have substantial influence on the horizontal component of magnetic field at ground. Present analysis reveals a time lag of ~30–45 min between the SW pressure changes seen at L1 location and ground magnetic field variations, and hence providing a good approximation of an averaged propagation time during entire storm interval; the time lag varies with solar wind velocity. Separate study during day and nighttime suggests that the SW dynamic pressure enhancements recorded by the dayside outer magnetospheric satellite have impact on the ground horizontal magnetic field measurements near equator, irrespective of day or nighttime.