维尼迪科夫调和分析对大甸子井水位潮汐因子和相位的气象影响研究

利用维尼迪柯夫潮汐调和分析方法计算内蒙古大甸子井水位原始资料和剔除气压、降雨资料后的潮汐因子和相位。结果显示,两种情况下得到的潮汐因子和相位具有较好的一致性,说明采用维尼迪柯夫潮汐调和分析方法得到的结果基本不受降水、气压等因素的影响。这对于缺少气压和降水等辅助观测的水位观测资料,使用原始资料进行维尼迪科夫潮汐调和分析,可获取比较可靠的潮汐因子和相位参数。     

用褶积滤波处理井水位对固体潮响应的滞后

本文给出了一种用褶积滤波处理井水位对固体潮响应滞后的新方法。鲁07井水位观测资料用这种方法处理比一般回归方法的平均中误差少4.6mm。文中讨论了褶积滤波积分区间长度与水井含水层渗透系数的关系及褶积滤波比别尔采夫滤波的优越性。     

香港地区重力固体潮和海潮负荷特征研究

介绍了在香港地区重力固体潮合作观测成果, 获得了该地区完整的重力固体潮实测模型. 利用全球和近海海潮模型以及岛屿验潮站数据较系统地研究了海潮负荷特征, 反演了全球海潮模型的适定性. 数值结果说明周日频段内的海潮模型要比半日频段内的模型更加稳定, 实施验潮站潮位高变化改正对精密确定重力固体潮相位滞后起重要作用. 文章还研究了重力观测残差和台站背景噪声水平. 本项研究填补了中国地壳运动观测网络在该地区重力固体潮观测空白, 为地表和空间大地测量提供有效参考和服务.     

地震子波处理的二步法反褶积方法研究

针对玛湖斜坡区三块三维地震资料和赛汉塔拉凹陷二块三维地震资料连片处理中的特点,结合地质任务和处理目标要求,提出了地震数据连片处理中的地震子波处理的方法.该方法主要体现了两次反褶积,一次是采用地表一致性反褶积,将不同震源的频带拓宽到一个标准上;再一次采用相位校正反褶积,将不同震源的数据校正到相同相位上.为了保证提取的相位校正反褶积算子稳定,采用叠后地震道提取（主要考虑到叠后地震道信噪比高,算子稳定性强）,然后将该算子应用到叠前地震道,进行相位校正.     

井约束非稳态相位校正方法

在地震勘探资料处理中,子波的零相位化有助于提高地震资料的分辨率、改善叠加剖面的质量.常规的相位校正方法是利用测井合成记录对井旁地震记录进行相位估计,然后对整条剖面进行常相位校正,该方法没有考虑地震子波相位的非稳态性（相位随时间和空间变化）.虽然通过局部相似度方法利用最大方差模准则或包络最大相似度准则可以估计出随时间和空间变化的相位属性,但是由于零相位判别准则本身具有一定的局限性,因此精度有限.针对这一问题,本文在局部地震属性和局部平面波模型下,提出了一种井约束的非稳态相位校正方法,该方法不仅考虑了子波相位的非稳态性,而且充分利用了测井合成地震记录进行相位校正精度高的优点.理论模型和实际资料处理表明,本文方法可以有效实现信号的零相位化,有利于改善叠加效果,提高资料的分辨率.     

临夏台钻孔系统性质的论证

对中国临夏台2013年至2014年两年的井水位和四分量钻孔应变资料进行了预处理以消除趋势与突跳.根据文献(Means,1982;Young and Budynas,2005),相互正交的两条测线的应变观测值之和等于面应变.文献(刘序俨等,1988)证明近地表的面应变的2/3等于体应变,因此,可由4分量钻孔应变观测值得到钻孔体应变,然后根据体应变与井水位观测资料,从时域和频域对该钻孔系统的性质进行了论证.结果表明,在时域,体应变与井水位高度负相关.钻孔系统的灵敏度为—0.1620mm/10-9.把两年中的某两个月份的两者时间坐标轴和纵轴比例尺放大,发现井水位曲线的峰/谷与体应变观测曲线的谷/峰一一对应,两者的相位滞后非常小.在频域内,本文采用Venedikov调和分析方法分别取得了井水位与体应变9个月的半日波与全日波数个波群的逐月潮汐因子与相位滞后,然后作简单计算,得到了钻孔系统对上述波群的灵敏度与相位滞后.结果表明9个月中大多数波群的灵敏度不但十分相近,且非常接近由时域得到的周年频率分量的灵敏度,但相位滞后误差较大,本文对此进行了分析,认为由反正切得到的相位滞后受计算误差影响较大,应以时域经审视所得的相位滞后接近于零为准.通过时域与频域的分析,表明井水位对体应变的响应基本是线性时不变的,论证了临夏台钻孔系统基本上满足了叠加性、齐次性与时不变性,基本上为一线性时不变系统.     

井水位的“记忆”滞后效应

观测资料表明,井水位对信息响应存在的“记忆”滞后现象,它与一般的位相滞后不同,在鲁２９井现场试验也证明了井水位对井吕空气压力变化的响应存在的“记忆”滞后现象。利用水平层状承压含水层模式,从理论上解释了井水位对井口空气压力变化响应“记忆”滞后现象,认为这种现象与水井含水层的导水系数有关,含水层导不系数越小,这种现象越明显,用一般多元回归方法无法较好地扣除井水位中“记忆”滞后影响,作给出了一种可以扣     

Quantitative correction of influence of rain on well water level

In this paper a new method of the quantitative correction of influence of the rain on the water level in a well is given. Using the faltung filtering and polynary regression, considering the effect of the rain on the well water level with the lag “memory”, the correction of the influence of the rain on the well water level is made. The result of correction of the water level in the well Lu-08 by the rain shows that the correction result of this method is better. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, 202–207, 1993. This study was supported by the Chinese Joint Seismological Science Foundation. Thank The Seismic Office of Victory Oil Field for their help in gathering data of the well Lu-08.     

福建省井水位与体应变观测一致性分析

本文从时域和频域分别对福州台和永安台的体应变和井水位资料观测的一致性进行了分析,并对其力学机制进行了深入探究,在此基础上阐明了二者观测一致性的物理含义.时域内的分析结果显示:福州台体应变与井水位的相关系数为0.7939,钻孔系统的灵敏度为0.0558 mm/10-9,井水位与体应变观测曲线十分相似;永安台体应变与井水位的相关系数为-0.7543,钻孔系统的灵敏度为-0.0005 mm/10-9,井水位与体应变曲线呈镜像对称,这表明两台站的井水位与体应变观测高度相关,且二者的相位滞后近似为零.在频域内,分别对福州台2012年3—5月及永安台2014年3—5月的井水位和体应变资料进行维尼迪柯夫调和分析,分别获取半日波、全日波数个波群的逐月潮汐因子和相位滞后,并计算得到了钻孔系统的灵敏度和井水位相对于体应变的相位滞后.结果表明,3个月中大多数波群的灵敏度十分相近,而且与时域内所得年周期低频分量的灵敏度较接近,但相位滞后误差较大,这可能与采用反正切公式计算相位滞后受计算误差影响有关,即非线性的反正切公式对计算误差具有放大累积效应,因此相位滞后应以时域分析所得接近于零的结果为准.综上,时域和频域的分析结果表明,福州台和永安台的井水位对体应变的响应基本上是线性时不变的,二者观测一致性的根本原因在于两台站的承压含水层井水位观测系统基本上是满足叠加性、齐次性和时不变的物理可实现因果系统.      

基于最高频相位法和空间滤波法的二维模型的大地电磁静态校正

静态效应一直是影响大地电磁测深法精确性的主要原因之一.在相位校正法的基础上提出了一种更适用于电性变化较为平缓的地质情况的静态校正方法—最高频相位法,其核心是用需要校正测点两侧受静态效应影响较小测点的最高频视电阻率的算术平均值代替相位法递推公式中每个频点前一个频点的视电阻率值,以消除相位法中的误差积累.以二维模型的正反演...     

基于自适应非稳态相位校正的时频域多尺度全波形反演

本文提出非稳态相位校正时频域目标函数,通过缩小观测数据与模拟数据在波形相位上的差异来缓解全波形反演过程中对应波形匹配错位的问题（周波跳跃）.同时引入自适应相位校正因子,可以根据观测数据与模拟数据的差异来调整相位校正量的大小.在构建非稳态相位校正时频域全波形反演目标函数的基础上,利用链式法则详细推导了对应的伴随震源,并从理论上证明了该方法的可行性与优越性.数值测试过程中结合了低通滤波多尺度反演策略,进一步缓解全波形反演过程中的强非线性问题.缺失低频分量测试结果表明,利用自适应非稳态相位校正时频域多尺度全波形反演方法结合常规全波形反演方法在缺失7 Hz以下低频分量的地震数据中仍然能够得到高精度的反演结果.震源不准确测试结果表明,即使震源子波相位差异较大,利用非稳态相位校正方法仍然能够一定程度上缓解周波跳跃现象.测试结果综合证明了本文提出的方法在构建初始速度建模,缓解周波跳跃等方面具有一定的优势.     

确定精密重力仪相位滞后方法的研究

从系统辨识的原理出发,给出了确定精密重力仪相位滞后的一种方法,即记录精密重力仪在阶跃信号作用下的输出数据,利用这些数据对仪器的传递函数进行辨识,从而确定重力仪对各潮汐波产生的相位滞后角度.用这种方法对ET-21地潮重力仪进行了相位滞后的确定,并给出这台仪器对15个主要潮汐波的相位滞后角度.     

鲁03井水位滞后于气压的处理

本文给出了一种处理井水位对气压影响滞后的方法,同时给出了计算公式。以鲁03井为例,用这种方法求出其井水位气压系数为6.41mm/hpa,水位对气压响应的滞后时间为1.52小时。该方法排除气压影响的效果较好,水位的气压校正中误差明显减小。     

PRINCIPLES AND APPLICATION OF MAXIMUM KURTOSIS PHASE ESTIMATION1

Methods of minimum entropy deconvolution (MED) try to take advantage of the non-Gaussian distribution of primary reflectivities in the design of deconvolution operators. Of these, Wiggins’(1978) original method performs as well as any in practice. However, we present examples to show that it does not provide a reliable means of deconvolving seismic data: its operators are not stable and, instead of whitening the data, they often band-pass filter it severely. The method could be more appropriately called maximum kurtosis deconvolution since the varimax norm it employs is really an estimate of kurtosis. Its poor performance is explained in terms of the relation between the kurtosis of a noisy band-limited seismic trace and the kurtosis of the underlying reflectivity sequence, and between the estimation errors in a maximum kurtosis operator and the data and design parameters. The scheme put forward by Fourmann in 1984, whereby the data are corrected by the phase rotation that maximizes their kurtosis, is a more practical method. This preserves the main attraction of MED, its potential for phase control, and leaves trace whitening and noise control to proven conventional methods. The correction can be determined without actually applying a whole series of phase shifts to the data. The application of the method is illustrated by means of practical and synthetic examples, and summarized by rules derived from theory. In particular, the signal-dominated bandwidth must exceed a threshold for the method to work at all and estimation of the phase correction requires a considerable amount of data. Kurtosis can estimate phase better than other norms that are misleadingly declared to be more efficient by theory based on full-band, noise-free data.     

Understanding the effect of elliptical polarization in surface nuclear magnetic resonance method

Currently, surface nuclear magnetic resonance (SNMR) method is the only geophysical method that detects groundwater directly. In this paper, we investigate the effect of elliptical polarization in the perpendicular excitation magnetic field. The effect of elliptical polarization is clearly visible in our ellipticity calculation and it can cause strong distortion to the excitation field in the presence of high subsurface conductivities. By examining the co-rotating and counter-rotating components of the field, we show that elliptical polarization affects transmitting and receiving processes differently and that a clear phase lag exists between transmitter loop and receiver loop. Finally, we derive the response function of coincident loops and calculate proton tip angles, the kernel function and SNMR response curves of a 1D aquifer model. Based on the simulations, we conclude that the elliptical polarization and phase lag can significantly affect SNMR response and it is essential to include elliptical polarization in SNMR modeling and data interpretation.     

基于稀疏反演三维表面多次波压制方法

三维表面多次波压制是海洋地震资料预处理中的重要研究课题,基于波动理论的三维表面多次波压制方法(3DSRME)是数据驱动的方法,理论上来说,可有效压制复杂构造地震数据表面多次波.但该方法因对原始地震数据采集要求高而很难在实际资料处理中广泛应用.本文基于贡献道集的概念,将稀疏反演方法引入到表面多次波压制中,应用稀疏反演代替横测线积分求和,无需对横测线进行大规模重建,进而完成三维表面多次波预测,这样可有效解决实际三维地震数据横测线方向稀疏的问题.基于纵测线多次波积分道集为抛物线的假设,为保证预测后三维表面多次波和全三维数据预测的多次波在运动学和动力学特征上基本一致,文中对预测数据实施基于稳相原理的相位校正.理论模型和实际数据的测试结果表明,本文基于稀疏反演三维表面多次波压制方法可在横测线稀疏的情况下,有效压制三维复杂介质地震资料中的表面多次波,从而更好地提高海洋地震资料的信噪比,为高分辨率地震成像提供可靠的预处理数据保障.     

黄石台ORBES-81型伸缩仪观测结果的调和分析及海潮改正

本文对黄石台ORBES-81型石英伸缩仪的应变固体潮观测资料作了分析计算,其调和分析结果为:对O₁波,振幅因子0.4219,相位滞后6°36＇;对M₂波,振幅因子0.4256,相位滞后-4°42＇.同时,讨论了海潮引起的负荷潮汐应变的改正问题.模拟计算发现负荷潮汐应变与体潮应变大致为同一量级,文中列出了海潮改正后的结果,并作了讨论.最后,对黄石台伸缩仪的架设方位作了评价.     

地震资料处理中小波函数的选取研究

本文给出了常见地震子波的一个模拟公式，可以很好地模拟零相位及混合相位子波，在一定意义上也可以近似模拟最大相位及最小相位子波．模拟出的子波加上适当的修正项后满足允许条件，可用作小波函数．与Ｍｏｒｌｅｔ小波类似，在实际应用中这些修正项在一定条件下可以略去，文中对Ｍｏｒｌｅｔ小波作了改造，使其能更好地适应于地震资料处理．研究了反射波能量及噪声等干扰波在时间－尺度域的分布特征与所选基本小波的关系．提出用地震子波（或与地震子波相近的函数）作为基本小波，对地震资料进行去噪及分频解释的方法．最后用实例证明方法的有效性．     

A two‐phase automatic static correction method

Statics are an effective approach to correct for complex velocity variations in the near surface, but so far, to a large extent, a general and robust automatic static correction method is still lacking. In this paper, we propose a novel two‐phase automatic static correction method, which is capable of handling both primary wave statics (PP statics) and converted‐wave statics (S‐wave statics). Our method is purely data driven, and it aims at maximizing stacking power in the target zone of the stack image. Low‐frequency components of the data are analysed first using an advanced genetic algorithm to estimate seed statics and the time structure for an event of interest, and then the original full‐band data are further aligned via the back‐and‐forth coordinate descent method using the seed statics as initial values and the time structure for event alignment guidance. We apply our new method to two field datasets, i.e., one for 2D PP static correction and the other for 3D S‐wave static correction.     

深井水位对固体潮和气压的响应

本文以体应变固体潮对深井水位影响的偏微分方程为基础，考虑到含水层与井孔之间相互渗流的边界条件，用叠加原理、冲量定理的分离变量等方法得出了方程的解。把水井含水层的参量给予一些可能的值，通过数值计算讨论了水井固体潮系数、位相滞后和含水层参数之间的关系，解释了井水位对固体潮响应的位相滞后现象。结果表明，井孔的半径、含水层的孔隙度及固体骨架的体压缩系数愈大，含水层的导水系数愈小，则水井水位的固体潮系数愈小