近地表低速带参数反演

阐明利用瑞雷面波反地表速度和厚度的可行性，并利用地震记录中的瑞雷面波数据，反演近地表低速带的厚度和速度，文中给出了实例试算，结果表明，面波资料的近地表参数反演将有利于勘探目标区的物理参数建模。     

对汶川地震小鱼洞地表破裂带的调查分析

2008年5月12日的汶川8.0级地震使龙门山断裂带形成了3条同震地表破裂带,这表明有多条活动断层同时参与地震破裂,其过程复杂,现象丰富。本文对小鱼洞地表破裂带及其与另2条地表破裂带的交汇区域进行了野外调查,并对小鱼洞地表破裂带的活动性质和展布特征进行了分析。小鱼洞地表破裂带位于彭州市小鱼洞镇附近,是汶川8.0级地震形成的一条走向NW的逆冲并具有左旋走滑分量的同震地表变形带。调查结果显示,小鱼洞地表破裂带表现出明显的分段性特征：小鱼洞镇一带的中段,逆冲量和走滑量最大;小鱼洞镇向东南方向延伸的南段,逆冲量和走滑量逐渐变小;小鱼洞镇向西北方向进入山区的北段,则表现为以逆冲为主的活动性质。     

Spatial dispersion of surface waves at the free surface of a general anisotropic elastic medium

A new technique relates the wave velocity of the surface waves in anisotropic elastic medium to its elastic constants. Anisotropic propagation of surface waves is studied in a half-space occupied by a general anisotropic elastic solid. The phase velocity expressions of quasi-waves, in three-dimensional space, are used to derive the secular equation of surface waves. The complex secular equation is resolved, analytically, into real and imaginary parts and is then solved, numerically, for phase velocity along a given phase direction on the surface. The complete procedure is thus analogous to the one used for conventional Rayleigh waves in isotropic medium. A non-linear equation relates the ray direction of the surface waves to its phase direction on the (plane) surface of the medium. The analytical differentiation of secular equation yields the directional derivative of phase velocity. This derivative is used to calculate the wave velocity of surface waves. Spatial variations of phase velocity, wave velocity and ray direction over the free plane surface are plotted for the numerical models of crustal rocks with orthorhombic, monoclinic and triclinic anisotropies.     

河北磁县西部山区最新地表破裂带的发现与1830年磁县7.5级地震的关系

本文介绍了磁县西部山区新近发现的最新地表破裂带。该地表破裂带近东西走向，长约７ｋｍ，由断层崖及地表沟槽组成。在地表粉末状灰岩破碎带中保存大量坚硬灰岩岩块，并见断层破碎带碎石盖在古腐植土之上。断层北盘下降，最新事件的最大垂直落差约７－１０ｍ，以正倾滑为主。个别地点断层显示左旋位移。该地表破裂带沿岔口活动断裂带分布，可认为是１８３０年磁县地震的地表破裂。     

Analysis on effect of surface fault to site ground motion using finite element method

Dynamic contact theory is applied to simulate the sliding of surface fault.Finite element method is used to analyze the effect of surface fault to site ground motions.Calculated results indicate that amplification effect is obvious in the area near surface fault,especially on the site that is in the downside fault.The results show that the effect of surface fault should be considered when important structure is constructed in the site with surface fault.     

浅层断层对场地地震动影响的有限元分析

应用动接触理论模拟断层的滑移问题，用有限元方法计算浅层断层对场地地震动的影响．计算结果表明，在地表断层出露的周围地区，断层下盘地表地震动放大效应较大，断层上盘地表地震动放大作用有明显的降低．这一结果说明，在有浅层断层的场地建设重大工程时需要考虑断层对地震动的影响．      

A Simple Relation to Constrain Groundwater Models Using Surface Deformation

A simple relation between pore pressure change and one-dimensional surface deformation is presented. The relation is for pore pressure change in a confined aquifer that causes surface deformation. It can be applied to groundwater models of any discretization and is computationally efficient. The estimated surface deformation from model results can be compared to observed surface deformation through geodetic techniques such as Differential Interferometric Synthetic Aperture Radar. Model parameters then are constrained using the observed surface deformation. The validity of this relation is shown through constraint of model parameters for surface uplift due to pore pressure increase caused by wastewater disposal injection.     

三峡水库水面蒸发年内滞后效应分析

水面蒸发在季节变化上相比净辐射等气象要素的相位延迟反映了水体储热对水面蒸发的影响,量化这一对气象要素波动的滞后响应对理解和估算深水水库(湖泊)蒸发非常重要。三峡水库等河道型深水水库的水位和面积具有显著周期性变动,使得其水面蒸发的响应模式更为复杂,而目前对其认识非常薄弱。本文利用2013年8月—2020年7月三峡水库巴东站水面蒸发场和陆面蒸发场的蒸发和气象观测数据,分析了水面蒸发的季节变化及其年内滞后效应。结果表明:水面蒸发场蒸发量在年内的8和12月呈现双峰变化,与只有8月单峰值的陆面蒸发场蒸发量显著不同。水面蒸发场蒸发量相对净辐射、平均气温和水面温度分别存在4、3和2个月的滞后,而陆面蒸发场蒸发量相对滞后时间均在1个月以内;水面与陆面蒸发场相比,水温、蒸发量和水面与大气饱和差之间的滞后时间分别为1、3和4个月,而平均气温和净辐射之间不存在滞后。本文揭示出三峡水库巴东段水面蒸发年内滞后效应主要受到水库水温引起的水面与大气饱和差在季节上滞后的影响,需通过深入分析水温的时空变化来明确整个三峡水库的水面蒸发特征。     

一种适合沙漠地区的静校正方法

静校正问题是我国西部复杂地形和复杂近地表地区地震勘探迫切需要解决的一个关键问题,也是一个十分复杂的问题.而复杂近地表地区的静校正问题往往归结为近地表层速度的求解问题,地表越复杂,近地表速度的求解就越困难.塔中地区表层地质条件变化剧烈,潜水面以上覆盖着大小不一的疏松沙体,沙体高度从几米到近百米,地震波在地表的传播速度随沙丘厚度而变化,因此,研究该区表层速度变化规律,成为该区静校正的主要问题.本文在总结前人所研究的沙漠地区静校正方法的基础上,提出了一种新的静校正方法,通过多项式拟合来求取表面速度,实际应用见到了良好的效果.     

起伏海面对反射地震偏移成像结果的影响

实际海面受风浪的影响是个随机起伏的复杂曲面.一个随机起伏的自由表面对地震波形成相当复杂的散射,进而影响波场在靶区的二次照明与成像结果;此外,起伏海面具有小尺度随机起伏的特性,难以用贴体网格等处理大尺度起伏地表的常规方法和技术对其进行逼近和处理.鉴于此,利用不等距有限差分法实现小尺度起伏海面的波场自由表面边界条件,以此进行正演、波场照明分析;利用逆时偏移（Reverse Time Migration,RTM）进行起伏海面下的成像并对成像结果进行分析.不同模型的试算结果表明：起伏海面引起的复杂散射使靶区地下照明不均匀、成像界面发生弯曲和畸变,从而降低成像结果的分辨率和信噪比.     

Utilization of satellite data in land surface hydrology: sensitivity and assimilation

This paper investigates the sensitivity of potential evapotranspiration to input meteorological variables, i.e. surface air temperature and surface vapor pressure. The sensitivity studies have been carried out for a wide range of land surface variables such as wind speed, leaf area index and surface temperatures. Errors in the surface air temperature and surface vapor pressure result in errors of different signs in the computed potential evapotranspiration. This result has implications for use of estimated values from satellite data or analysis of surface air temperature and surface vapor pressure in large‐scale hydrological modeling. The comparison of cumulative potential evapotranspiration estimates using ground observations and satellite observations over Manhattan, Kansas for a period of several months shows a variable difference between the two estimates. The use of satellite estimates of surface skin temperature in hydrological modeling to update the soil moisture using a physical adjustment concept is studied in detail, including the extent of changes in soil moisture resulting from the assimilation of surface skin temperature. The soil moisture of the 1 cm surface layer was adjusted by 0·9 mm over a 10‐day period as a result of a 3 K difference between the predicted and the observed surface temperature. This is a considerable amount given the fact that the top layer can hold only 5 mm of moisture. Copyright © 2001 John Wiley & Sons, Ltd.     

北极夏季冰面上近地层特征及热量收支问题

应用1999年8月19～24日我国第1次北极考察资料对北极夏季近地层气象要素梯度特征和冰面热量收支问题进行了讨论.北极夏季近地层相对于冰面的相对湿度很大,经常接近于饱和状态,冰温具有明显的日变化,夜间冰面附近冰温梯度较大,白天很小.在冰温垂直分布中20cm冰层温度有时可出现极大值.冰面热量收支计算表明,夜间冰面辐射冷却损失热量主要由冰层向冰面热量输送来补偿.白天冰层的热量主要来源于穿透短波辐射,冰面和以下冰层间的热量交换很低.白天冰面冰雪融化热量是不可忽视的.     

岩石和断层泥摩擦特性的现场大尺度试件实验研究

为了研究岩石与断层泥的摩擦特性,在远安断裂带附近的灰岩体中开挖试洞、制作试体,进行大尺度试件的现场实验研究.实验用直剪法进行.试体的尺度为506030cm3.试体与岩体间的摩擦面积为3000cm2.摩擦面分为天然断层面、软弱夹层面和节理裂隙面三种.测得的摩擦系数平均值是:（1）天然断层面————0.48;（2）软弱夹层面————0.41;（3）节理裂隙面————0.84.与室内小试件实验结果进行对比,发现现场断层泥的摩擦系数与室内实验结果相近.现场岩石的摩擦系数与室内粗糙岩石的实验结果相近.由于摩擦面积的增大,摩擦面上各点介质屈服的时间有明显的差别.部分摩擦面上介质的较早屈服可作为预报整体失稳的前兆.      

含流体孔隙介质中面波的传播特性及应用

基于单相介质中地震波理论的高频面波法已广泛应用于求取浅地表S波的速度.然而水文地质条件表明,普遍的浅地表地球介质富含孔隙.孔隙中充填的流体会显著地影响面波在浅地表的传播,进而造成频散和衰减的变化.本文研究了地震勘探频段内针对含流体孔隙介质边界条件的面波的传播特性.孔隙流体在自由表面存在完全疏通、完全闭合以及部分疏通的情况.孔隙单一流体饱和时,任何流体边界条件下存在R1模式波,与弹性介质中的Rayleigh波类似,相速度稍小于S波并在地震记录中显示强振幅.由于介质的内在衰减,R1在均匀半空间中也存在频散,相速度和衰减在不同流体边界下存在差异.Biot固流耦合系数(孔隙流体黏滞度与骨架渗透率之比)控制频散的特征频率,高耦合系数会在地震勘探频带内明显消除这种差异.介质的迂曲度等其他物性参数对不同流体边界下的R1波的影响也有不同的敏感度.完全闭合和部分疏通流体边界下存在R2模式波,相速度略低于慢P波.在多数条件下,如慢P波在时频响应中难以观察到.但是在耦合系数较低时会显现,一定条件下甚至会以非物理波形式接收R1波的辐射,显示强振幅.浅表风化层低速带存在,震源激发时的运动会显著影响面波的传播.对于接收点径向运动会造成面波的Doppler频移,横向运动会造成面波的时频畸变.孔隙存在多相流体时,中观尺度下不均匀斑块饱和能很好地解释体波在地震频带内的衰减.快P波受到斑块饱和显著影响,R1波与快P波有更明显关联,与完全饱和模型中不同,也更易于等效模型建立.频散特征频率受孔隙空间不同流体成分比例变化的控制,为面波方法探测浅地表流体分布与迁移提供可能性.通常情况孔隙介质频散特征频率较高,标准线性黏弹性固体可以在相对低频的地震勘探频带内等效表征孔隙介质中R1波的传播特征,特别在时域,可在面波成像反演建模中应用.     

三维地震与地面微地震联合校正方法

由于地面微地震监测台站布设在地表,会受到地表起伏、低降速带厚度和速度变化的影响,降低了微地震事件的识别准确度和定位精度,限制了地面微地震监测技术在复杂地表地区的应用.因此,将三维地震勘探技术的思路引入到地面微地震监测中,提出了三维地震与地面微地震联合校正方法,将油气勘探和开发技术更加紧密地结合在一起.根据三维地震数据和低降速带测量数据,通过约束层析反演方法建立精确的近地表速度模型,将地面微地震台站从起伏地表校正到高速层中的平滑基准面上,有效消除复杂近地表的影响.其次,根据射孔数据和声波测井速度信息,通过非线性反演方法建立最优速度模型,由于已经消除复杂近地表的影响,在进行速度模型优化时不需要考虑近地表的影响,因而建立的速度模型更加准确.最后,在精确速度模型的基础上,通过互相关方法求取剩余静校正量,进一步消除了复杂近地表和速度模型近似误差的影响.三维地震与地面微地震联合校正方法采用逐步校正的思路,能够有效消除复杂近地表的影响,提高微地震数据的品质和速度模型的精确度,保证了微地震事件的定位精度,具有良好的应用前景.     

Directional statics in common reflection surface stack for rugged surface topography

Seismic data acquired along rugged topographic surfaces present well‐known problems in seismic imaging. In conventional seismic data processing, datum statics are approximated by the surface consistence assumption, which states that all seismic rays travel vertically in the top layer. Hence, the datum static for each single trace is constant. In case this assumption does not apply, non‐constant statics are required. The common reflection surface (CRS) stack for rugged surface topography provides the capability to deal with this non‐vertical static issue. It handles the surface elevation as a coordinate component and treats the elevation variation in the sense of directional datuming. In this paper I apply the CRS stack method to a synthetic data set that simulates the acquisition along an irregular surface topography. After the CRS stack, by means of the wavefield attributes, a simple algorithm for redatuming the CRS stack section to an arbitrarily chosen planar surface is performed. The redatumed section simulates a stack section whose acquisition surface is the chosen planar surface.     

Application of a sigma coordinate sea model to the calculation of wind-induced currents

A three-dimensional numerical sea model is formulated in terms of sigma coordinates in the vertical. The vertical grid spacing in the model is arbitrary and can be refined to give enhanced resolution in high shear regions (e.g., close to the sea surface in wind-driven flows, and/or across the thermocline in stratified flows). A method of accurately determining surface currents and indicating how fine a grid is required in the surface layer is described.The problem of determining a suitable formulation of vertical eddy viscosity to use in a model of wind-induced flow in a tidal sea is considered in detail. A formulation in which surface eddy viscosity depends upon the roughness of the sea surface and the transfer of momentum to depth by surface waves appears reasonable. Below the surface layer turbulence is related to the current at depth.Idealized calculations are performed to demonstrate the accuracy and stability of the sigma coordinate model. Results of these calculations indicate that the formulation of eddy viscosity developed in this paper can explain the high surface shears reported in lake measurements of wind-induced surface currents, and the lack of shear under strong wind conditions in the open sea (GORDON, 1982, Journal of Geophysical Research, 87, 1939–1951).Surface current to surface wind ratio are also computed.     

一个大气边界层非线性三维数值模式及其在太湖流域的应用

本模式用以K半经验理论为基础的闭合方程组求解,引进了地形作用、动量、热量和水汽的侧向扩散和垂直湍流交换以及凝结加热等物理过程,其中湍流交换系数K是地表粗糙度、几何高度、大气稳定度及风速切变等因子的函数。本模式还考虑了山地、平原和水面在地形高度、地表粗糙度、辐射、蒸发及热交换等方面的差异,在下垫面上建立热平衡方程与边界层内控制方程耦合。将模式应用于太湖流域,计算得到的该流域边界层内温、压、湿、风的分布特征与实际情况相似,模式具有一定的实用性。     

3D surface‐wave estimation and separation using a closed‐loop approach

Surface waves in seismic data are often dominant in a land or shallow‐water environment. Separating them from primaries is of great importance either for removing them as noise for reservoir imaging and characterization or for extracting them as signal for near‐surface characterization. However, their complex properties make the surface‐wave separation significantly challenging in seismic processing. To address the challenges, we propose a method of three‐dimensional surface‐wave estimation and separation using an iterative closed‐loop approach. The closed loop contains a relatively simple forward model of surface waves and adaptive subtraction of the forward‐modelled surface waves from the observed surface waves, making it possible to evaluate the residual between them. In this approach, the surface‐wave model is parameterized by the frequency‐dependent slowness and source properties for each surface‐wave mode. The optimal parameters are estimated in such a way that the residual is minimized and, consequently, this approach solves the inverse problem. Through real data examples, we demonstrate that the proposed method successfully estimates the surface waves and separates them out from the seismic data. In addition, it is demonstrated that our method can also be applied to undersampled, irregularly sampled, and blended seismic data.     

Model‐based attenuation for scattered dispersive waves

Coherent noise in land seismic data primarily consists in source‐generated surface‐wave modes. The component that is traditionally considered most relevant is the so‐called ground roll, consisting in surface‐wave modes propagating directly from sources to receivers. In many geological situations, near?surface heterogeneities and discontinuities, as well as topography irregularities, diffract the surface waves and generate secondary events, which can heavily contaminate records. The diffracted and converted surface waves are often called scattered noise and can be a severe problem particularly in areas with shallow or outcropping hard lithological formations. Conventional noise attenuation techniques are not effective with scattering: they can usually address the tails but not the apices of the scattered events. Large source and receiver arrays can attenuate scattering but only in exchange for a compromise to signal fidelity and resolution. We present a model?based technique for the scattering attenuation, based on the estimation of surface‐wave properties and on the prediction of surface waves with a complex path involving diffractions. The properties are estimated first, to produce surface?consistent volumes of the propagation properties. Then, for all gathers to filter, we integrate the contributions of all possible diffractors, building a scattering model. The estimated scattered wavefield is then subtracted from the data. The method can work in different domains and copes with aliased surface waves. The benefits of the method are demonstrated with synthetic and real data.