从昌平井体应变、水位对地震波的响应特征求算含水层的Skempton常数

以理想孔隙弹性介质线性孔隙弹性理论为基础, 利用井含水层系统孔隙压力变化与体应变之间的数学表达式. 以2004年12月26日苏门答腊MＳ8.7（根据中国地震台网中心目录）地震波引起的昌平地震台水位和体应变波动资料为例,分别从时间域和频率域分析了水位与体应变之间的相互关系,并计算出地震波通过含水层时引起的水位波动对体应变的响应系数,给出了一种估算Skempton常数B的方法,为求解含水层特性提供了途径.      

Precursory groundwater level changes in the period of activation of the weak intraplate seismic activity on the NE margin of the Bohemian Massif (central Europe) in 2005

tWe analyse continuous measurements of groundwater level in two deep wells VS-3 and V-28 at the experimental hydro-meteorological station situated on the NE margin of the Bohemian Massif, central Europe, characterized by the weak intraplate seismic activity. The aim of our study is to examine the relationships between changes in the groundwater level and earthquake occurrence. Based on the tidal and barometric response of the water level, we estimated selected elastic parameters of the observed aquifers: the shear modulus G, the Skempton ratio B, the drained matrix compressibility β and the undrained compressibility β_u. Using these parameters and assuming the homogeneous poroelastic material, we derived the sensitivity of the wells to the crustal volume strain. During the observation period from November 1998 to December 2005 we detected in the VS-3 well two pre-seismic steps, related to August 10, 2005 (M = 2.4) and October 25, 2005 (M = 3.3) earthquakes. Amplitudes of the recorded precursory changes (+6 cm and +15 cm) are several times higher than the values predicted from the theoretical precursory crustal strain and the strain sensitivity of the well. Therefore, we presume that the observed pre-seismic water level steps can be attributed to heterogeneity of poroelastic material. We consequently propose the hypothesis of the origin of precursory events based on the presumption of a sensitive site, at which the well is situated.     

Study on Volume Strain Inversion from Water Level Change of Well-aquifer Systems

Based on linear poroelasfic and hydrogeology theory, a mathenmatical expression describing the relationship between water level clmnge and aquifer volume strain is put forward. Combined with earth tidal theory, we analyze the response characteristics from well-aquifer water level change to earth tide of volume strain and present a method of volume strain inversion from water level clmnge. Comparing the results of inversion with real observed data, we found that there is a good consistency. This suggests that the method of volume strain inversion from water level clmnge is proper. It will offer a reference for learning about hydrogeology characteristics, volume strain and searching for precursor anomalies.     

由井-含水层系统的水位动态反演含水层体应变

以孔隙弹性理论和水文地质学原理为基础,给出井水位波动与含水层介质体应变变化关系的数学表达式,结合固体潮理论分析了井—含水层系统水位潮汐波动对体应变固体潮的响应特征,用水位潮汐波来反演含水层体应变的变化,将反演结果与体应变实测资料得到的结果进行对比,发现二者有较好的一致性,说明用水位潮汐波来反演含水层体应变的方法是可行的,它为了解含水层的水文地质特性、体应变的变化及探索地震前兆提供了借鉴。     

Analysis of reservoir heterogeneity-induced amplification effect on time-lapse seismic responses of fluid substitution: A physical modelling study

Seismic wave propagation through a fluid-saturated poroelastic layer might be strongly affected by media heterogeneities. Via incorporating controlled laboratory simulation experiments, we extend previous studies of time-lapse seismic effects to evaluate the wave scattering influence of the heterogeneous nature of porous permeable media and the associated amplification effects on 4D seismic response characteristics of reservoir fluid substitution. A physical model consisted of stratified thin layers of shale and porous sandstone reservoir with rock heterogeneities was built based on the geological data of a real hydrocarbon-saturated reservoir in Northeast China. Multi-surveys data of good quality were acquired by filling poroelastic reservoir layers with gas, water and oil in sequence. Experimental observations show that reservoir heterogeneity effect causes significantly magnified abnormal responses to the fluid-saturated media. Specifically, reflection signatures of the gas-filled reservoir are dramatically deviated from those of the liquid fluid-filled reservoir, compared with ones of the homogeneous media. By removing the influences unrelated to reservoir property alterations, 4D seismic estimates of travel-time and frequency-dependent characteristic are reasonably consistent with fluid variations. Nevertheless, strong 4D amplitude difference anomalies might not correspond to the regions where fluid variations occur. We also find that 4D seismic difference attributes are evident between oil- and water-filled models, whereas significant between oil- and gas-filled models. Meanwhile, rock physics modelling results reveal the predicted 4D seismic differences are obviously smaller than those calculated from seismic observations. The results in this paper, therefore, implicate that the effect of a reservoir's heterogeneous nature might be beneficial for hydrocarbons detection as well as monitoring small variations in pore fluids.     

汤池一号井水震波与地震波的对比分析

水震波是地震波引起的井孔水位的振荡现象.利用同一地质单元内的水震波与地震面波资料进行对比研究,分析出水震波与地震波数据上的异同.同时运用地震波与水震波的数据异同,探求两者数据在应用方面的联系,以更好理解两者的对比关系,本课题估算出水位对含水层应力应变的响应"格值",通过对庐江地震台等井的井水位观测资料分析,结合白山数字...     

Reaearch on Skempton's coefficient B based on the observation of groundwater of Changping station

Based on isotropic linear poroelastic theory and under the undrained condition, we summarize three equations connecting the Skempton's coefficient B with the groundwater level. After analysis, we propose a method to calculate the Skempton's coefficient B according to the relationship between water level and tidal strain. With this method we can get the value of B without the earthquake occurrence, which can provide the high frequency waves for research. Besides, we can also get the in-suit Skempton's coefficient B without the experiment of rock physics. In addition, we analyze the observed data of Changping station recorded in groundwater monitoring network (abv., GMN) before and after the Wenchuan M_S8.0 with this method, and find out there's a slight change of the value of B after the seismic waves passed by, which implies that the propagation of seismic waves may have brought some variations to the poroelastic medium of the well.

     

Hydrologic response of underground reservoirs to seismic vibrations

The response of the water level in wells to seismic impacts from remote earthquakes and explosions are analyzed. It is shown that in most cases the magnitude of the postseismic change in the water level scales as a square root of the amplitude of the deformation wave. The intensity of persistent changes averages a dynamical deformation of 1–5 cm/μstrain. Noticeable deviations from the mentioned range are possible, depending on the particular structural features of the layer.     

Soil pressure and pore pressure for seismic design of tunnels revisited: considering water-saturated,poroelastic half-space

This paper describes a systematic study on the fundamental features of seismic soil pressure on underground tunnels, in terms of its magnitude and distribution, and further identifi es the dominant factors that signifi cantly infl uence the seismic soil pressure. A tunnel embedded in water-saturated poroelastic half-space is considered, with a large variety of model and excitation parameters. The primary features of both the total soil pressure and the pore pressure are investigated. Taking a circular tunnel as an example, the results are presented using a fi nite element-indirect boundary element(FE-IBE) method, which can account for dynamic soil-tunnel interaction and solid frame-pore water coupling. The effects of tunnel stiffness, tunnel buried depth and input motions on the seismic soil pressure and pore pressure are also examined. It is shown that the most crucial factors that dominate the magnitude and distribution of the soil pressure are the tunnel stiffness and dynamic soil-tunnel interaction. Moreover, the solid frame-pore water coupling has a prominent infl uence on the magnitude of the pore pressure. The fi ndings are benefi cial to obtain insight into the seismic soil pressure on underground tunnels, thus facilitating more accurate estimation of the seismic soil pressure.     

Influencing factors of seismic signals generated by un-tuned large volume airgun array in a land reservoir

Un-tuned large volume airgun array in a water reservoir is recently proposed as a new way to generate seismic waves on land. It can be used to explore the earth velocity structure and its temporal variations as well. However, the characteristics of seismic signals (especially far-field signals) from an airgun array in a reservoir and its affecting factors (firing pressure, airgun towing depth, water level of the reservoir, etc.) has not been adequately studied. We analyzed the seismic data collected from field experiments at Binchuan Transmitting Seismic Station in 2011 and 2013 and found that (1) The similarity of seismic signals decrease with distance, which is most likely induced by the decay of signal amplitude and signal to noise ratio (SNR); (2) The amplitudes of far-field airgun signals are almost linearly proportional to the firing pressure; (3) The towing depth of airgun has less effects on the far-field signals; (4) The amplitudes of far-field airgun signals are proportional to the water level of the reservoir.     

福州地震台TJ-2型体积式应变仪安装及注水试验

在安装福州地震台TJ-2型体积式应变仪过程中,发现存在钻孔选点与观测室分开考虑等问题,提出一些建议。井孔注水试验结果揭示,该井孔探头顶部可能存在一个与探头近距离的含水层干扰。初步对气压、水位、降雨影响进行分析。     

Co-Seismic Groundwater Level Changes Induced by the May 12, 2008 Wenchuan Earthquake in the Near Field

The large scales of co-seismic water level changes in mainland China were observed in response to the tragic 2008 Ms 8.0 Wenchuan earthquake. To better understand the mechanism of these hydrogeological phenomena, groundwater-level data at 17 confined wells, with an epicentral distance of <500 km, were collected. We compare the static strain predicted by dislocation theory with the volumetric strain calculated by the tide effect of the groundwater based on poroelastic theory. The results show that the sign of the co-seismic groundwater level change is consistent with the sign predicted by dislocation theory. Additionally, the magnitude of the strain calculated by the two methods is also concordant in half of the wells. In the rest of the wells, the strains inversed from the groundwater level are one or two orders of magnitude larger than the fault dislocation model. These wells mostly have an epicenter distance larger than 300 km; therefore, the dynamic stress induced by the seismic wave may be responsible for the co-seismic water level changes in these wells. According to these results, we roughly estimate that the effect range of the static stress is approximately 300 km for the Wenchuan earthquake, and the dynamic stresses dominate beyond this epicenter distance. In addition, geological and hydrogeological conditions and other mechanisms may be responsible for these changes.     

Application of Biot's poroelasticity to some soil dynamics problems in civil engineering

This review type of paper shows how the poroelastodynamic theory of Biot can be applied to some soil dynamics problems encountered in transportation engineering, which have been solved by the present authors. These problems involve rigid walls retaining poroelastic soil and subjected to harmonic seismic waves and moving loads on poroelastic soil. Both classes of problems involve a soil layer over bedrock, are of the plane strain type and are solved analytically by two methods: a direct (almost exact and exact for the above two classes of problems) method and an approximate method. The effects of shear modulus, porosity, permeability and hysteretic damping of the soil medium as well as the seismic frequency for retaining walls and velocity for moving loads on the dynamic response are numerically evaluated in order to assess their relative importance on that response.     

宽城地震台体应变仪观测数据干扰因素及其映震能力分析

本文通过分析宽城地震台体应变仪自运行以来所有的观测数据,对宽城地震台体应变仪运行良好的区段（2005年-2008年）的动态特征进行分析研究,在此基础上对影响宽城地震台体应变仪的地下水位、温度、气压等相关因素进行了相关性分析,得出了宽城地震台体应变仪数据主要受到地下水位和气压的影响,温度则对宽城地震台体应变仪影响较小。通过观测数据拟合分析来看,宽城地震台体应变仪映震能力较好,体应变仪在震情监测预报以及地震研究工作起到重要的作用,对地震监测预报工作的长远发展做了较大的贡献。     

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.     

Numerical simulation of frequency-dependent seismic response and gas reservoir delineation in turbidites: A case study from China

To simulate the frequency-dependent response of turbidite reservoirs in the JZ Area, near the Bohai Sea, China, we apply a diffusive and viscous wave equation (DVWE), which takes into account diffusive and viscous attenuation and velocity dispersion in fluid-bearing poroelastic media. We use a seismic data-driven geological model building approach to produce physical parameter sections, which are then used to numerically synthesize the frequency-dependent seismic response in the DVWE-based simulation. The DVWE-based synthetic section shows the characteristic reflection and geometry of turbidites and delineates the phase delay, instantaneous dominant frequency decrease and magnitude attenuation related to the gas-bearing reservoir. The common frequency sections obtained by instantaneous spectral decomposition of the synthetic section show that a low-frequency shadow (LFS) lies immediately beneath the reservoir. Next, following the implications of the numerical simulation, we then apply LFS and fluid mobility to the data volume. Both hydrocarbon indicators clearly delineate the bright gas reservoir and its spatial distribution. The workflow and methodologies can be expected to be applicable to other frequency-dependent hydrocarbon indicators.     

青岛地震台钻孔体应变异常变化因素分析

分时段研究降雨、钻孔水位、工程施工等因素对青岛地震台体应变的影响，研究结果表明:钻孔水位反映区域地下水主要由降雨补给，降雨下渗改变地下水状态，并影响体应变观测；青岛体应变与钻孔水位的相关性及其周期特性均在2015年前后出现变化，可能与2015—2016年降雨偏少及工程施工有关；2018年1月17日台站钻孔施工很可能改变了体应变周围岩石孔隙环境，并扩大孔隙压力作用系数，可解释2018年3月以来水位小幅上升引起的体应变显著压性上升现象，受钻孔施工的影响，体应变与钻孔水位相关性的滞后特征明显减弱。     

应变能积累在地震安全性评价中的应用探讨

理论及实验研究显示,震源区的位置和最大震级不仅可能与构造和历史有关,更重要的是与应力状态包括其方向、大小、增加速率、集中位置等以及岩体自身性质,包括其强度、杨氏模量、应变大小及增长速率、裂缝(断层)大小和数量等有关,同时与作用方式有关.综合考虑各种因素,探讨了利用应变能积累确定强震位置和震级的方法:首先对利用形变空间特征变化判断强震震源区位置的方法进行了讨论;接着着重提出估算应变能的初步方法:一是由地表裂缝最大错距和裂缝(断层)长度比估算地震释放能量;二是由新构造运动分区块体体积、年平均应变速率及杨氏模量估算块体应变能积累;进而依据能量估算地震震级;最后以实例论述了应变能积累过程及特点对地震安全性评价的重要作用.     

模糊数学方法在分析日本东海地区地震前兆中的应用

本文把模糊模式识别的直接方法用于检测日本东海地区地震活动性、体积应变、地下水位、地下水氡含量、地倾斜和潮位中的前兆变化。此方法可消除诸如降雨和气温变化等干扰,使我们能更清晰地识别出地震前兆。此外,还把模糊相似方法和模糊聚类方法应用于识别地震空区及检验地震活动随时间变化的相似性。上述方法在日本其他地区同样有效,即将编入日本气象研究所板内地震预报计划中的计算机系统内。