The influence of mesoscopic flow on the P-wave attenuation and dispersion in a porous media permeated by aligned fractures

In sedimentary rocks attenuation/dispersion is dominated by fluid-rock interactions. Wave-induced fluid flow in the pores causes energy loss through several mechanisms, and as a result attenuation is strongly frequency dependent. However, the fluid motion process governing the frequency dependent attenuation and velocity remains unclear. We propose a new approach to obtain the analytical expressions of pore pressure, relative fluxes distribution and frame displacement within the double-layer porous media based on quasi-static poroelastic theory. The dispersion equation for a P-wave propagating in a porous medium permeated by aligned fractures is given by considering fractures as thin and highly compliant layers. The influence of mesoscopic fluid flow on phase velocity dispersion and attenuation is discussed under the condition of varying fracture weakness. In this model conversion of the compression wave energy into Biot slow wave diffusion at the facture surface can result in apparent attenuation and dispersion within the usual seismic frequency band. The magnitude of velocity dispersion and attenuation of P-wave increases with increasing fracture weakness, and the relaxation peak and maximum attenuation shift towards lower frequency. Because of its periodic structure, the fractured porous media can be considered as a phononic crystal with several pass and stop bands in the high frequency band. Therefore, the velocity and attenuation of the P-wave show an oscillatory behavior with increasing frequency when resonance occurs. The evolutions of the pore pressure and the relative fluxes as a function of frequency are presented, giving more physical insight into the behavior of P-wave velocity dispersion and the attenuation of fractured porous medium due to the wave-induced mesoscopic flow. We show that the specific behavior of attenuation as function of frequency is mainly controlled by the energy dissipated per wave cycle in the background layer.     

Formation of systems of incompact bands parallel to the compression axis in the unconsolidated sedimentary rocks: A model

Uniaxial compression of poorly lithified rocks leads to the formation of thin incompact layers (or bands, in the two-dimensional case) parallel to the compression axis, which are characterized by increased porosity. The standard model of the formation of such bands, as well as deformation bands of other types, associates them with the narrow zones of localization of plastic deformations. In the case of decompaction, it is assumed that transverse tensile deformations are localized within the band, which cause the band to dilate. Here, the formation of a band of localized deformations is treated as a loss-of-stability phenomenon. Based on observations, we propose a fundamentally different model of incompact bands formation, according to which the microdefects in sediment packing (pores) rather than the deformations are localized in the narrow zones. The localization of pores, which are initially randomly distributed in the medium, occurs as a result of their migration through the geomaterial. The migration and subsequent localization of pores are driven by a common mechanism, namely, a trend of a system to lower its total energy (small variations in total energy are equal to the increment of free energy minus the work of external forces). Migration of a single pore in a granular sedimentary rock is caused by the force f driving the defect. This force was introduced by J. Eshelby (1951; 1970). An important feature of our model is that the formation of an incompact band here does not have a sense of a loss of stability. Quite the contrary, the formation of incompact bands is treated as a gradual process spread over time. In this context, the origination of incompact band systems directly follows from our model itself, without any a priori assumptions postulating the existence of such systems and without any special tuning of the model parameters. Moreover, based on the proposed model, we can predict the incompact bands to always occur in the form of systems rather than as individual structures. A single incompact band may only be formed when the force resisting the pore motion, f _c , is absent.     

应用远震有限频率层析成像反演首都圈上地幔速度结构

本文选用首都圈数字地震台网2003年9月~2005年12月记录的300多个远震事件的波形资料,采用分频带多道互相关方法得到三个不同频段的P波相对走时数据共18499个,计算了每个频段的走时灵敏度核,应用有限频率层析成像反演得到首都圈地区的上地幔三维P波速度结构模型.利用检测板估计了反演结果的分辨率,并与射线层析成像方法的结果进行了比较,说明了反演结果的可靠性.研究结果表明,各构造单元具有明显不同的速度结构特征,其差异可到150 km深:燕山隆起区表现高速;太行山隆起区整体以低速为主并存在小范围高速块体;华北盆地、渤海湾下浅层上地幔中存在大范围的强低速异常,其顶面在50~70 km,可视为软流圈顶面的埋深,这一结果说明华北盆地、渤海湾下岩石圈明显减薄;张家口-蓬莱断裂带是上地幔浅部速度结构的变异带,也是岩石圈减薄的边界带,区内大部分强震都发生在该构造带上,由此看来该带上强震的发生不仅与地壳结构的不均匀性有关,还可能有较深的构造背景.     

储层地震物理模拟的震源指向性及介质非均质性影响分析

随着地震物理模拟的研究目标由地质构造发展为不同类型的油气储层,明确震源指向性及介质非均质性对实验结果的影响,对提高物理模拟数据的可信性及准确性具有重要意义.利用激光测振新方法对震源指向性进行了实验测试,得到了更加精确的指向性结果,形成了模拟数据振幅补偿方法.制作了具有不同尺度特征的非均质样品,在50 kHz至1000 kHz超声频段下对纵波速度及能量进行了测试,明确了不同的波长/介质尺度比(λ/a)范围内纵波速度及干扰波与有效波能量比的变化特征.结果表明,在地震物理模拟超声频段下射线速度向散射速度的转折点出现在λ/a≈2处;散射速度向等效介质速度的转折点出现在λ/a≈100处.当满足λ/a>200条件时,干扰波与有效波能量比小于5%时,可认为满足储层模拟的等效介质条件.     

2013年芦山地震震源区地壳介质地震波速变化的特征分析

为探究芦山M7.0级地震后5年多来,震源区龙门山断裂带西南段介质波速的变化规律,本文基于2012年4月至2018年4月共6年的连续波形数据,运用移动窗互谱与频域偏振等分析方法,结合背景噪声源的特性,对不同深度范围内的相对波速变化以及震后的恢复过程与机制进行了研究.获得的主要认识包括：（1）年尺度而言,震源区周期为1~20 s的背景噪声场相对稳定,但成分复杂、2~10 s频带内至少存在2个能量相对稳定的噪声源;不同周期噪声的能量,在月变与季节性上的变化特征差异明显.（2）获得了长时间尺度、不同频带内介质相对波速的背景变化水平,1~2 s、2~4 s的波动幅度（约为±0.04%）与季节性变化规律强于4~10 s、10~20 s的,结合与降雨量相关的地下水位模型能很好地解释其变化规律.（3）震源区的同震波速降低现象清晰,降幅约为0.08%~0.1%;空间上,波速下降最为显著的区域主要集中在龙门山断裂带两侧约70 km范围内,其中四川盆地一侧平均约为0.1%,略高于青藏高原（0.08%）一侧;在断裂带内的降速不显著.对不同子频带进行测量的结果显示,震后除10~20 s外,其余3个子频带的相对波速在震后较短时间内（约20天左右）均出现较大幅度的波速降低现象,其中4~10 s的平均降速最大（约为0.08%）,分析认为主震及大量余震的松弛效应是引起介质波速下降的主要原因.（4）震后大约1年左右,波速变化基本恢复到震前水平,且至2018年4月前未观察到大幅的波速变化现象,总体上各频带内的结果均沿零线小幅波动.     

Wave dispersion in high‐rise buildings due to soil–structure interaction

Nonparametric techniques for estimation of wave dispersion in buildings by seismic interferometry are applied to a simple model of a soil–structure interaction (SSI) system with coupled horizontal and rocking response. The system consists of a viscously damped shear beam, representing a building, on a rigid foundation embedded in a half‐space. The analysis shows that (i) wave propagation through the system is dispersive. The dispersion is characterized by lower phase velocity (softening) in the band containing the fundamental system mode of vibration, and little change in the higher frequency bands, relative to the building shear wave velocity. This mirrors its well‐known effect on the frequencies of vibration, i.e. reduction for the fundamental mode and no significant change for the higher modes of vibration, in agreement with the duality of the wave and vibrational nature of structural response. Nevertheless, the phase velocity identified from broader band impulse response functions is very close to the superstructure shear wave velocity, as found by an earlier study of the same model. The analysis reveals that (ii) the reason for this apparent paradox is that the latter estimates are biased towards the higher values, representative of the higher frequencies in the band, where the response is less affected by SSI. It is also discussed that (iii) bending flexibility and soil flexibility produce similar effects on the phase velocities and frequencies of vibration of a building. Copyright © 2014 John Wiley & Sons, Ltd.     

Biot流动和喷射流动耦合作用下波传播的FCT紧致差分模拟

本文从含流体多孔隙介质中同时包含Biot流动和喷射流动两种力学机制的BISQ(Biot-Squirt)方程出发,利用FCT(Flux-Corrected Transport)紧致差分方法数值模拟了在Biot流和喷射流共同作用下的波在含流体多孔隙各向同性介质中的传播.通过与仅受Biot流动作用下的波场结果对比,我们研究了...     

Ground roll wave suppression based on wavelet frequency division and radial trace transform

Ground roll waves interfere with seismic data. The suppression of ground roll waves based on the division of wavelet frequencies considers the low-frequency characteristics of ground roll waves. However, this method will not be effective when the ground roll wave and the effective signal have the same frequency bands because of overlapping. The radial trace transform (RTT) considers the apparent velocity difference between the effective signal and the ground roll wave to suppress the latter, but affects the low-frequency components of the former. This study proposes a ground roll wave suppression method by combining the wavelet frequency division and the RTT based on the difference between the ground roll wave velocity and the effective signal and their energy difference in the wavelet domain, thus making full use of the advantages of both methods. First, we decompose the seismic data into different frequency bands through wavelet transform. Second, the RTT and low-cut filtering are applied to the low-frequency band, where the ground roll waves are appearing. Third, we reconstruct the seismic record without ground roll waves by using the inverse RTT and the remaining frequency bands. The proposed method not only improves the ground roll wave suppression, but also protects the signal integrity. The numerical simulation and real seismic data processing results suggest that the proposed method has a strong ability to denoise while preserving the amplitude.     

饱和储层砂岩中流体粘度引起的超声波速度频散实验研究(英文)

在实验室对5块储层砂岩进行了模拟地层压力条件下的超声波速度测试。砂岩样品采自WXS凹陷的W地层,覆盖了从低到高的孔隙度和渗透率范围。实验选用了卤水和4种不同密度油作为孔隙流体,结合温度变化,实现了对流体粘度引致的速度频散研究。对实验结果的分析表明:(1)对于高孔隙度和渗透率的样品,无论是哪种流体饱和,观察到的超声波速度测试值和零频率Gassmann预测值的差异较小(约2-3%),基本上可以用Biot模型解释;对于中等孔隙度和渗透率的样品,低粘度流体(<约3mP?S)的频散效应也可以用Biot模型得到合理解释;(2)对于低、中孔隙度和渗透率样品,当流体粘度增加时,喷射流机制起主导作用,导致严重的速度频散(可达8%)。对储层砂岩的微裂隙纵横比进行了估计并用于喷射流特征频率的计算,当高于该特征频率时,Gassmann理论的假设条件受到破坏,实验室测得的高频速度不能直接用于地震低频条件下的W地层砂岩的Gassmann流体替换研究。     

基于格林函数理论的波场预测和鬼波压制方法

作为一种特殊的噪声,鬼波对一次波的波形及频带宽度产生极大的影响,鬼波压制是提高海上地震资料分辨率及保真度的重要因素.以格林公式为基础,详细论述了基于格林函数理论的鬼波压制方法,在不需要地下介质信息的条件下,进行地震数据驱动鬼波压制,并根据"Double Dirichlet"(双狄利克雷)边界条件,预测压力波场和垂直速度波场.建立了基于格林函数理论鬼波压制的处理流程,数值模拟和实际资料处理结果表明,基于格林函数理论鬼波压制方法在很好地去除鬼波的同时极大地拓宽了地震资料的频带,尤其提升了低频端能量,有利于后续资料的处理解释.     

五株微囊藻毒素降解菌的分离与鉴定

有毒蓝藻产生的蓝藻毒素导致地表水的污染已是全球普遍的现象,最常见的毒素是能导致肝癌的肝毒素微囊藻毒素(MC),迄今为止尚无适用的MC驱除技术,微生物降解MC可能是一条有效的途径．本研究以铜绿微囊藻(Microcystis aeruginosa PCC 7806)所产生的MC为底物,从多个形成微囊藻水华的水体中分离获得能显著降解MC的细菌18株,相对降解能力从18．4％到48．2％．对其中5株MC相对降解能力较强的菌,通过形态特征和生理生化特征分析,鉴定出5株菌分别为:不动细菌属(Acinetobacter)、肠杆菌属(Enterobacter)、微杆菌属(Microbacterium)、芽孢杆菌属(Bacillus)和弗拉特氏菌属(Frateuria),均不同于已报道的属种．     

远场地震波辐射能计算以及对近断层地面运动的约束

从圆盘断层模型出发,根据地震波能量表象定理推导出了滑移弱化过程中远场S-波辐射能量表达式,并同已有的动力学模型作了比较.结果表明,得到的模型能量值或视应力的取值强烈地依赖于断层上的动态、静态应力降和破裂传播速度,而破裂速度则对应了断层带模型中断层破坏过程所耗散的能量.动摩擦应力上调和应力下调力学机制在能量求解中得到了充...     

2.5D modelling, inversion and angle migration in anisotropic elastic media

2.5D modelling approximates 3D wave propagation in the dip‐direction of a 2D geological model. Attention is restricted to raypaths for waves propagating in a plane. In this way, fast inversion or migration can be performed. For velocity analysis, this reduction of the problem is particularly useful. We review 2.5D modelling for Born volume scattering and Born–Helmholtz surface scattering. The amplitudes are corrected for 3D wave propagation, taking into account both in‐plane and out‐of‐plane geometrical spreading. We also derive some new inversion/migration results. An AVA‐compensated migration routine is presented that is simplified compared with earlier results. This formula can be used to create common‐image gathers for use in velocity analysis by studying the residual moveout. We also give a migration formula for the energy‐flux‐normalized plane‐wave reflection coefficient that models large contrast in the medium parameters not treated by the Born and the Born–Helmholtz equation results. All results are derived using the generalized Radon transform (GRT) directly in the natural coordinate system characterized by scattering angle and migration dip. Consequently, no Jacobians are needed in their calculation. Inversion and migration in an orthorhombic medium or a transversely isotropic (TI) medium with tilted symmetry axis are the lowest symmetries for practical purposes (symmetry axis is in the plane). We give an analysis, using derived methods, of the parameters for these two types of media used in velocity analysis, inversion and migration. The kinematics of the two media involve the same parameters, hence there is no distinction when carrying out velocity analysis. The in‐plane scattering coefficient, used in the inversion and migration, also depends on the same parameters for both media. The out‐of‐plane geometrical spreading, necessary for amplitude‐preserving computations, for the TI medium is dependent on the same parameters that govern in‐plane kinematics. For orthorhombic media, information on additional parameters is required that is not needed for in‐plane kinematics and the scattering coefficients. Resolution analysis of the scattering coefficient suggests that direct inversion by GRT yields unreliable parameter estimates. A more practical approach to inversion is amplitude‐preserving migration followed by AVA analysis. SYMBOLS AND NOTATION A list of symbols and notation is given in Appendix D .     

Scaling exponents of the velocity structure functions in the interplanetary medium

We analyze the scaling exponents of the velocity structure functions, obtained from the velocity fluctuations measured in the interplanetary space plasma. Using the expression for the energy transfer rate which seems the most relevant in describing the evolution of the pseudo-energy densities in the interplanetary medium, we introduce an energy cascade model derived from a simple fragmentation process, which takes into account the intermittency effect. In the absence and in the presence of the large-scale magnetic field decorrelation effect the model reduces to the fluid and the hydromagnetic p-model, respectively. We show that the scaling exponents of the q-th power of the velocity structure functions, as obtained by the model in the absence of the decorrelation effect, furnishes the best-fit to the data analyzed from the Voyager 2 velocity field measurements at 8.5 AU. Our results allow us to hypothesize a new kind of scale-similarity for magnetohydrodynamic turbulence when the decorrelation effect is at work, related to the fourth-order velocity structure function.     

PAHs contamination and bacterial communities in mangrove surface sediments of the Jiulong River Estuary, China

Sixteen sediment samples collected from eight transects in a mangrove swamp of the Jiulong River Estuary, Fujian, China were investigated for their content of polycyclic aromatic hydrocarbons (PAHs) and the biodegradation potential of the indigenous microorganisms. The bacterial community structures in the mangrove sediments and in enrichment cultures were also investigated. The results showed that the total PAHs concentration of mangrove sediments ranged from 280 to 1074 ng g(-1) dry weight, that the PAHs composition pattern in the mangrove sediments was dominated by high molecular weight PAH components (4-6 rings), and that Benzo[ghi]perylene and Indeno[1,2,3-cd]pyrene were the most dominant at different stations. Abundant PAH-degrading bacteria were found in all the stations, the values of phenanthrene-degrading bacteria ranged from 5.85 x 10(4) to 7.80 x 10(5) CFU g(-1) dry weight, fluoranthene-degrading bacteria ranged from 5.25 x 10(4) to 5.79 x 10(5) CFU g(-1) dry weight, pyrene-degrading bacteria ranged from 3.10 x 10(4) to 6.97 x 10(5) CFU g(-1) dry weight and the benzo(a)pyrene-degrading bacteria ranged from 5.25 x 10(4) to 7.26 x 10(5) CFU g(-1) dry weight. DGGE analysis of PCR-amplified 16S rDNA gene fragments confirmed that there was a remarkable shift in the composition of the bacterial community due to the addition of the different model PAH compound phenanthrene (three ring PAH), fluoranthene(four ring PAH), pyrene(four ring PAH) and benzo(a)pyrene(five ring PAH) during enrichment batch culture. Eleven strains were obtained with different morphology and different degradation ability. The presence of common bands for microbial species in the cultures and in the native mangrove sediment DNA indicated that these strains could be potential in situ PAH-degraders.     

Periodic modulation of Pc3 and Pc4 pulsations in the polar cap by interplanetary and atmospheric processes

The variations in the daily average energy of geomagnetic pulsations and noise in the Pc3 (20–60 mHz) and Pc4 (10–19 mHz) frequency bands in the polar cap have been studied based on the data from P5 Antarctic station (corrected geomagnetic latitude ?87°) from November 1998 to November 1999. The daily average pulsation energy has been calculated using the method for detecting the wave packets, the spectral amplitude of which is higher than the threshold level, from the dynamic spectrum. A spectral analysis of the energy of pulsations and noise in the Pc3 and Pc4 bands, performed using the maximal entropy method, has revealed periodicities of 18 days in the local winter and 26, 13, and 7–9 days during the local summer. The simultaneous and coherent variations with periods of 26, 13, and 7–9 days in the solar wind velocity and IMF orientation indicate that the variations in the Pc3–4 wave energy in the polar cap at a sunlit ionosphere are mainly controlled by the parameters of the interplanetary medium. The variations in the Pc3–4 wave energy with a period of 18 days are observed only during the local winter and are supposedly related to the variations in the ionospheric conductivity modulated by planetary waves.     

Dynamics of planktonic bacteria and heterotrophic microflagellates in the Parker Estuary, northern Massachusetts

The temporal and spatial abundance and activities of members of the planktonic microbial food web in a salt marsh estuary were analysed over the course of a year. Analysis of the data was focused on the factors that control the density and productivity of the heterotrophic bacteria. Temperature exercises a long-term influence on substrate sources, but short-term control is a function of the interaction of grazing and substrate supply. A model has been developed to give a framework for understanding how the microbial food web may work under steady-state assumptions. The model assumptions and predictions appear to be supported by the field data. Equations were generated for estimating substrate flux through the bacteria and estimating bacterial density as a result of grazing and substrate input. The results give some insights into the question of whether the heterotrophic bacteria function as a link to higher organisms or as an energy sink.     

Nonlinear isotropic elastic reduced Cosserat continuum as a possible model for geomedium and geomaterials. Spherical prestressed state in the semilinear material

We suggest a nonlinear elastic reduced Cosserat continuum as a possible model for geomedium and geomaterials and also for a soil or rock with heterogeneities. If a medium has a block structure, or if it contains heterogeneities that may have their proper rotational dynamics, the presence of rotational degrees of freedom may influence wave propagation and stability of the medium. In the reduced Cosserat model, translations and rotations are independent, the medium resists to the rotation of each ??particle?? relatively to the background continuum, but it does not resist to the gradient of rotation. We consider a nonlinear spherical stress state in an isotropic elastic reduced Cosserat material. We write down small deviations from this nonlinear equilibrium. They coincide in form with the equations of the linear elastic isotropic reduced Cosserat continuum. Depending on the level of the stress and on the type of elastic energy, the equilibrium can be stable or unstable. In the domain of stability, shear?Crotational wave demonstrates dispersion. There is a resonant frequency corresponding to the independent rotational oscillations. The bulk plane shear?Crotational wave has a forbidden band of frequencies. In this zone, the shear?Crotational wave localises near heterogeneities or external sources. We show that for a semilinear medium (medium with large deformations but linear dependence of stress on the strain tensor), strong compression leads to the material instability caused by shear perturbations, and strong tension for some class of parameters yields in instability caused by rotational perturbations.     

周期性层状含孔隙、裂隙介质模型纵波衰减特征

地震波在含孔隙、裂隙斑块饱和介质传播过程中会诱发多个尺度孔隙流体流动而产生衰减和速度频散.在含有宏观尺度“Biot流”和介观尺度“局域流”衰减诱导机制的周期性层状孔隙介质模型基础上,引入了微观尺度硬币型和尖灭型裂隙“喷射流”的影响,构建了周期性层状含孔隙、裂隙介质模型.利用双解耦弹性波动方程的方法数值计算了该模型地震频带的纵波衰减和速度频散并与周期性层状孔隙介质模型做了对比研究.分析了该模型在不同裂隙参数（裂隙密度、裂隙纵横比）及裂隙体积含量下的纵波衰减和频散特征,裂隙密度越高对于纵波衰减和频散的影响越大,裂隙纵横比越小,由裂隙引起的纵波衰减部分向高频段移动,裂隙体积含量越少,纵波衰减先降低后小幅增加再降低,频散速度增加,并逐渐接近于周期性层状孔隙介质模型的纵波衰减和频散速度曲线.最后研究了周期性层状含孔隙、裂隙介质模型有效平面波模量的高低频极限以及流固相对位移在该模型中的分布特征.     

A new approach to develop the Raman carbonaceous material geothermometer for low‐grade metamorphism using peak width

The Raman spectra of carbonaceous material (CM) from 19 metasediment samples collected from six widely separated areas of Southwest Japan and metamorphosed at temperatures from 165 to 655°C show systematic changes with metamorphic temperature that can be classified into four types: low‐grade CM (c. 150–280°C), medium‐grade CM (c. 280–400°C), high‐grade CM (c. 400–650°C), and well‐crystallized graphite (> c. 650°C). The Raman spectra of low‐grade CM exhibit features typical of amorphous carbon, in which several disordered bands (D‐band) appear in the first‐order region. In the Raman spectra of medium‐grade CM, the graphite band (G‐band) can be recognized and several abrupt changes occur in the trends for several band parameters. The observed changes indicate that CM starts to transform from amorphous carbon to crystallized graphite at around 280°C, and this transformation continues until 400°C. The G‐band becomes the most prominent peak at high‐grade CM suggesting that the CM structure is close to that of well‐crystallized graphite. In the highest temperature sample of 655°C, the Raman spectra of CM show a strong G‐band with almost no recognizable D‐band, implying the CM grain is well‐crystallized graphite. In the Raman spectra of low‐ to medium‐grade CM, comparisons of several band parameters with the known metamorphic temperature show inverse correlations between metamorphic temperature and the full width at half maximum (FWHM) of the D1‐ and D2‐bands. These correlations are calibrated as new Raman CM geothermometers, applicable in the range of c. 150–400°C. Details of the methodology for peak decomposition of Raman spectra from the low to medium temperature range are also discussed with the aim of establishing a robust and user‐friendly geothermometer.