消除地震前兆观测中干扰的一种DAI数字滤波器

本文介绍了一种DAI数字带通滤波器,这种滤波器能根据数据处理需要选择滤波参数。参数确定后,可计算出滤波器的频率响应,根据频率响应的好坏确定是否重新选择滤波参数,直到满意为止。滤波的误差可以在滤波器设计中给定,适当增加褶积的项数即滤波器的长度,可以使滤波的误差很小。本文用这种方法处理了地下水位观测资料,消去固体潮汐的影响。其效果比用潮汐理论值校正的方法消去潮汐影响更好。本文还用频移的方法将这种带通滤波器变成低通滤波器,用来消去视电阻率的年变化。直接用带通滤波提取水氡观测值中周期在半个月到两年的变化成分。这种滤波方法在以上应用中均得到了较好的效果。     

基于褶积原理和改进广义S变换的震动波衰减补偿试验研究

为了恢复震动波能量在传播过程中产生的衰减损耗,提出基于褶积原理求取品质因子Q的方法与改进广义S变换相结合的反Q滤波法。通过震动波衰减补偿模型试验,对试验数据进行改进广义S变换的时频特性分析,得出了信号的能量分布情况以及时间频率对应关系;采用基于褶积原理求取品质因子的方法,得到时变Q值;对试验数据进行反Q滤波处理,使震动波能量得到了补偿。结果表明本文提出的反Q滤波法提高了对震动波能量衰减补偿的效果,拓宽了地震资料的频带,提高了地震资料分辨率,有利于进行高分辨率地震勘探、深部信号增强和油气藏预测工作的开展。     

新丰江短水准资料中的降雨影响

短水准资料中的降雨干扰相当可观.本文详细讨论了这种干扰的影响机制和影响特点,尝试用褶积运算方法建立单位降雨量影响模式和降雨总影响模型,并提出运用褶积滤波方法消除测值中降雨影响的具体方法.作者以新丰江测量站的资料为例进行了整套计算。结果表明,此模式能真实地再现资料受降雨影响的情况.用此法消除降雨影响后,测值恢复了正常的年变化形态,提高了资料的分析使用价值.     

混合法计算GB体波理论地震图

高斯射线束（GB）方法是一种用于计算不均匀介质中波场的高频近似方法。本文在详细讨论了几种用高斯射线束叠加计算理论地震图的方法--频谱法、褶积法和波包法之后,提出了适用范围更广泛的混合方法（褶积-波包法）,并给出了一个便于数值计算的褶积公式。混合法在计算GB理论地震图时既用褶积法又用波包法,可以得到较高的计算速度与精度。最后,就二维情况分别用褶积法、波包法和混合法进行了一些实际计算和比较。     

聊古-1井氮氦比的动态灰箱分析

聊古-1井位于聊考断裂带的北段.多年来的观测分析表明,影响该井氮氦比变化的干扰最主要的是气象因素降雨和井水流量.其中流量的影响存在着“记忆”滞后效应.用一般相关分析难以将这种影响排除干净.本文从系统论、控制论出发,对有记忆影响的流量因子采用褶积滤波方法进行处理,以最小二乘法为建立模型的判据,用动态灰箱分析法拟合并分析了观测系统的变化.在分析中用了4个状态变量:长趋势成分、年周期变化成分、井水流量影响部分和随机变化部分.     

聊古-1井氮氦比区间概率信息量与地震的关系

采用了非线性的褶积滤波结合多元回归的方法排除聊古１井氮氦比观测值中的干扰影响,然后计算了其观测值区间概率的信息量．分析了信息量与１９８１年宁晋ＭＳ５．８地震和１９８３年菏泽ＭＳ５．９地震的关系．结果显示,在这两次地震前２个月左右,氮氦比信息量都明显地出现升高异常,地震后恢复到正常值     

时频域时变反褶积方法

常规反褶积方法具有很多局限性,往往要求地震子波是最小相位并且是平稳的,反射系数序列为白噪等等.本文将常规褶积模型扩展成时变褶积模型,通过S变换对地震记录进行谱分解,估计出震源子波和衰减因子,求出时变子波,在S域进行时变反褶积,再反S变换到时间域输出结果.此方法完全打破了常规反褶积方法的局限性,数值试验和实际资料处理均证明了此方法的有效性.     

用褶积滤波及多元回归方法处理聊古—1井氮氦比观测资料

聊古-1井位于聊城——兰考断裂带北段。自1981年以来对该井氮氦比每日进行观测,其主要干扰因素为降雨和井水流量。其中降雨对氮氢比观测的影响不是简单的线性关系,而是带有滞后记忆效应的非线性关系。为此,我们采用了非线性的褶积滤波的方法,来处理降雨对氮氦比的影响。为了消除井水流量变化对氮氦比的影响,同时采用了多元回归的方法,对氮氦比观测值进行综合处理。在以上处理的基础上,分析了1981年宁晋地震和1983年菏泽地震前,氮氦比的变化特征。两次地震前的2个月左右,氮氦比值都出现了明显地高值异常,且均在异常转折、恢复的过程发震。所以这一规律可做为该井氮氦比地震短期前兆异常指标。     

聊古－1井氮氦比观测值的区间概率与地震的关系

本文介绍了计算聊古－１并氮氦比观测值区间概率的两种方法。采用了非线性的褶积滤波结合多元回归的方法排除氮氟比观测值中的干扰因素。并分析了它与１９８１年宁晋Ｍｓ５．８级地震和１９８３年菏泽Ｍｓ５．９级地震之间的关系。在这两震前２个月，其两种区间概率都明显地出现异常，而且都在曲线的极值处发震，震后又恢复到正常值。     

褶积与反褶积的一种计算方法

应用FFT方法计算实型线性褶积与反褶积问题时,通常需要将线性褶积转化为k=1的循环褶积,因而所得到的计算结果是近似的并且产生较为严重的边界效应问题.本文提出了一种新方法,它通过将线性褶积转化为k=±-1的循环型褶积,然后引入FFT进行计算;在保证不增加计算量的情况下,所得到的实型褶积运算是精确的、反褶积运算中边界效应可以得到较好的改善.     

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

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

苏05井数字化水位雨荷效应的处理与改正

根据苏05井水位受降雨影响具有滞后“记忆”效应的特点,利用褶积滤波方法,对苏05井水位的雨荷效应进行处理和改正。处理结果表明,井水位与固体潮理论值的相关系数由处理前的0．24l上升为改正后的0．65l,达到未受降雨影响的水平,因此对于受降雨影响严重的井孔水位进行降雨量的改正是有效而又极其必要的。     

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

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

深井水位的固体潮效应

本文从体应变固体潮对深井水位影响的偏微分方程出发,考虑到含水层和井孔之间相互渗流的边界条件,用叠加原理、冲量定理和分离变量法等方法得出了方程的解.通过对这个解中水井含水层参数给予一些可能的值进行数值计算,讨论了水井固体潮系数和位相滞后与水井含水层参数间的关系,较好地解释了井水位对固体潮响应的位相滞后现象.计算表明,井孔的半径、含水层的孔隙度及固体骨架的体压缩系数愈大,含水的导水系数愈小,则水井水位的固体潮系数愈小,而水位对固体潮响应的位相滞后愈大.井水对长周期的潮汐响应比对短周期的更好.      

Effect of earth tide on deep well water level

On the basis of the partial differential equations of the water level in a deep well that is affected by the bulk strain tide, the boundary conditions of permeating each other between aquifer and well are considered. The solutions of the equations have been obtained by the superposition principle, inpulse theorem and separation of the variations. Some suitable values are given to the parameters of well aquifer in the solutions. By means of the numerical calculation the relations among the well tide coefficients and phase lag phenomenon of response of the well level to earth tides is well interpreted. The parametric studies and the obtained graphs show that the more the radius of casing in a well, the porosity of an aquifer and the bulk compressibility for solid matrix of an aquifer, the less the transmissibility of an aquifer, then the less the tidal coefficient of well level, while the more the phase lag for well level to respond to the tide is. The response of the well level to the tide of longer period is better than the shorter period events. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 66–75, 1991.     

地下水潮汐现象的物理机制和统一数学方程

用新的分层承压含水层模式 ,不但考虑含水层的力学压缩性质 ,而且考虑含水层的渗流特性 ,并结合扰动信息源的频率特性 ,分别研究扰动源地球固体潮、大气潮和地表负荷潮对承压井水位和流量的影响机理 ,给出相应的偏微分方程。从方程的解释或数值解讨论扰动源与承压井含水层的力学压缩参数、渗流特性参数及与频率特性频数的关系 ,进而给出承压井水位和流量对地球固体潮、大气潮和地表负荷潮汐响应的统一数学方程及其潮汐响应函数 ,并揭示了上述几类潮汐扰动信息源对承压井水位和流量影响的物理机理     

中国地下水潮汐的观测研究和分析

本文简单介绍了中国地下水潮汐的观测研究情况,介绍了对承压水位潮汐机理研究的理论结果,给出了固体潮体应变、气压和海洋潮汐对承压水头影响的偏微分方程,从这些方程的简化解讨论了它们与含水层参数间的关系.并从全国选取了33口有代表性的井孔,对这些井孔的水位和气压观测资料采取反复调和分析的方法,求出了它们的水位固体潮体应变各分波的潮汐系数和位相滞后,也求出了水位气压各分波的系数和位相滞后.      

鲁07井水位对固体潮滞后的处理

本文给出了一种计算考虑滞后影响的井水位固体潮系数的新方法，导出了计算固体潮系数和滞后时间的公式，给出了考虑滞后影响的井水位固体潮校正公式。以鲁07井为例，用该法求出其井水位固体潮系数为16.5mm/10^-8，水位对固体潮响应的滞后时间为0.48小时，与以胶的回归分析方法相比，校正后的水位中误差明显减小，校正效果比较好。     

菏泽5＼9级地震前豫01井水位对气压和固体潮响应的变化

本文给出了一种处理井水位的新方法，它考虑了井水位对气压和固体潮响应的滞后影响。利用泰勒和杜得森展开法，导出了计算井水位气压系数，固体潮系数以及它们的滞后时间的公式，并给出了考虑滞后影响的井水位气压与固体潮改正公式。我们利用这一新方法处理了豫０１井１９８３年的观测资料，计算出每天的井水位固体潮，气压系数和它的滞后时间。该４种参数在１９８３年１１月７日菏泽Ｍｓ５．９级地震时大幅度变化。而且在地震前，后     

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

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