2018年8月5日印度尼西亚地震震前电离层异常特征

基于2018年8月5日印度尼西亚地震震中（116.45°E,8.33°S）附近的澳大利亚达尔文站和位于震中磁力线共轭区域的武汉站的电离层测高仪观测数据,以及相同区域全球卫星导航系统接收机的观测数据,对该地震震前的电离层异常扰动特征进行了分析。结果显示,震前同时观测到了电离层F₂层临界频率（f_oF2）和总电子含量（TEC）时间序列的异常。基于岩石圈-大气层-电离层直流电场耦合模式和美国国家大气研究中心的热层-电离层环流耦合模式（NCAR/HAO TIEGCM）对震前震中及其对应磁力线共轭点区域的电子密度异常进行了模拟,模拟结果表明,在地震异常电场的作用下,地震区域及其对应半球的磁共轭区域的TEC和电离层F₂层峰值电子密度N_mF2发生了明显扰动。      

二维最佳线性数字滤波器的设计原理

针对如何在干扰场的背景上区分出低缓异常,以及在位场的向下延拓一类计算中如何限制因误差的高频放大所导至的解的不稳定性等问题,本文探讨了在“最小二乘”意义下的最佳线性数字滤波器的设计原理,并将它转化为下述数学问题,即在L₂线性赋范函数空间中如何选取最佳滤波函数的问题。在空间域中直接解这个问题是十分复杂和困难的,我们发现在波数域中用变分法中的等周问题的解法直接选取最佳线性滤波器的传输函数（或波数响应）,则在数学方法上既简单又严格。这样选取的最佳线性滤波器的传输函数L（f,k）其表达式也很简单,即L（f,k）=|S_i（f,k）|²/{|S_i（f,k）|²+λ|N_i（f,k）|²}。式中,|S_i（f,k）|²及|N_i（f,k）|²分别代表滤波器输入端讯号和干扰的能谱（或功率谱）,f、k分别代表x、y方向上的波数,λ为大于零的常数。 对上述两类问题以及相关的两种最佳线性滤波器而言,L（f,k）的表达式是相同的,而区别仅在于其参变量λ的选取条件不同而已。 有了最佳线性滤波器的传输函数L（f,k）的理论公式,就可以在最小二乘的意义下分析和评价国内外所发表的解决上述两类问题的各种线性滤波方法,并能指出在不同的讯号与干扰条件下,在理论上线性滤波可能达到的最佳效果,从而为设计二维线性数字滤波器时,提供一个理论上的准则。 对位     

电阻率测深的数字解释

本文主要介绍了应用积分变换的方法和采样定理将视电阻率ρ_s曲线作线性滤波,得出一新的电阻率转换函数T′曲线,然后,以层参数（各层的电阻率和厚度）算得的T用最优化数值方法在DJS-6型电子计算机上与其进行自动拟合,以达到解释电测深曲线的目的。 文中简述了戈什（Ghosh）提出的对ρ_s作线性滤波的原理,介绍了与国外不同的取样间距和滤波系数的确定以及阻尼最小二乘法和变尺度最优化法的计算框图和应用,最后附有实例和简要的讨论。     

随机激发过程对地极移动的作用

以随机微分方程组（1）作为地极移动（钱德勒摆动）的数学模型,求出了（1）的解。此解是二维的正态、平稳、马尔科夫随机过程,它具有遍历性.找到了它的转移概率密度和相关函数.通过对相关函数的研究,如果取ω=5.274,则钱德勘周期重合于强相关周期2π/ω;此外,还发现了另一个周期,即独立性周期π/ω.最后,研究了地极移动模型的最佳预测公式和预测误差.     

2016年熊本MJ7.3地震的工程地震动参数模拟及分布特征分析

本文选取2016年4月16日日本熊本县M_J7.3 （M_W7.0）地震近场区域内K-net地震台网的47个强震台站所记录的加速度数据,运用经验格林函数法模拟分析此次地震主要的工程地震参数并给出了各工程参数的空间分布。通过对比分析得到结论如下:① 地震动时程的基本频谱的模拟结果较好,尤其是1—15 Hz的高频段内;② 在震源距小于50 km范围内,峰值加速度和阿里亚斯强度的观测值与模拟值拟合较好,二者在近场区域均以椭圆形向周围扩散衰减,且阿里亚斯强度的模拟值整体大于观测值;③ 卓越周期拟合整体较好,但在局部区域存在较大差异,模拟结果难以表征场地环境的复杂性。      

多尺度地球化学勘查聚焦找矿靶区——以滇东Pt地球化学省勘查研究为例

铂矿资源一直是中国紧缺而急需的矿产资源.以滇东Pt地球化学勘查为例,在极低密度泛滥平原沉积物勘查数据圈定的滇东Pt地球化学省的基础上,首先按不同的勘查尺度对研究区内国家地球化学填图计划采集的水系沉积物样品进行重组形成组合样品,然后应用化学预富集光谱法（C-OES）分析组合样的Pt含量,并应用变差函数定量描述不同勘查尺度数据圈定的Pt异常的变化性及控制其变化性的地质因素,最后基于不同勘查尺度Pt数据分别圈定区域和局部Pt异常块体,按地球化学块体方法估算各异常块体内的Pt资源潜力.结论表明：（1）从区域Pt异常块体到局部Pt异常块体,控制因素由一级裂控-岩控向二级岩控-裂控转化,表明变差函数有效揭示了控制不同勘查尺度Pt异常块体变化性的主要地质因素.（2）存在两类Pt异常：一类是在两个观测尺度下都具有Pt异常显示,且异常含量随观测尺度的变小（取样密度的增加）而显著增强;另一类是随勘查尺度的变小Pt异常弱化甚至消失.表明多尺度地球化学聚焦找矿具有逐步排除非矿异常,筛选矿异常的功能.（3）小尺度勘查数据圈出的局部Pt异常块体的Pt平均含量是大尺度勘查数据圈出的区域Pt异常块体的2倍,局部Pt异常块体Pt资源量占区域Pt异常块体的60%,而前者面积仅占后者的35%.表明Pt资源量主要集中在局部异常区域,聚焦找矿具有高效浓缩靶区的勘查效果.     

大地电磁资料精细处理和二维反演解释技术研究（四）——阻抗张量分解的多测点-多频点统计成像分析

基于大地电磁阻抗张量分解技术,本文提出了两种电性主轴方位的统计描述图像：随频率变化的统计分布成像（频率分布云图）和随测点序列变化的统计分布成像（测点分布云图）.这两种图像与传统的统计玫瑰图一起,较全面地描述了最佳主轴的分布特征.在进行构造维性分析过程中,通过定义二维有效因子e_2d,来压制一维结构和三维结构、突出纯二维结构的影响.e_2d被用于电性主轴的统计加权,有效地起到了滤波的作用;同时,统计成像中还考虑了数据质量的影响.为了得到稳定、高质量的区域阻抗张量数据,提出并实现了共主轴的多测点-多频点阻抗张量分解新算法.最终,完成了以上各项处理手段的可视化实现.本文通过两个理论模型和一个实测算例,以共轭阻抗法（CCZ法）为基础,展示了这一新技术的有效性.     

2013年芦山MS7.0地震前南北地震带数字化水位、 水温高频信息异常特征及效能分析

数字化观测资料中含有丰富的高频异常信息, 这些高频信息中可能隐含着与地震有关的异常信号。 引入概率密度分布法对2013年4月20日芦山M_S7.0地震发生前南北地震带水位、 水温分钟值观测数据高频成分进行分析, 结果表明, 芦山M_S7.0地震前南北地震带共有17个测点、 19组数据出现高频异常信息, 其中水位5组、 水温14组, 高频异常多集中在四川南部构造带, 且随着时间进程发展, 异常幅值逐渐减小且异常出现向震中区迁移的特征。 通过应用Molchan图表法对提取的高频信息进行统计检验, 结果显示, 概率密度分布法提取的高频异常对地震具有较好的预测效果。 通过对震源区及附近区域大地构造演化、 前兆异常发展及地壳运动速度场等特征分析, 结果发现, 计算得到的地下流体高频异常空间分布特征与南北地震带地壳运动场具有很好的一致性, 表明利用概率密度分布法提取的高频异常与区域构造活动存在一定相关性, 同时也为前人关于对异常形成机理的研究提供了旁证。     

陕西地磁短周期变化异常解释

1996年9月从陕西南部石泉至甘肃宁县布设地磁短周期变化观测临时测线．通过测点地磁转换函数,分析了该区地磁短周期变化特征．根据子波变换进行多尺度分析的理论,提取了该区地磁短周期变化异常的区域性源的影响成分．并以此为观测数据,反演研究了测线地下相应于区域性源的二维导电率结构．在反演研究中,将Zhdanov等提出的基于观测电磁场和导电率模型理论计算场的剩余场总能量最小,并且引入迁移剩余场概念的反演方法,引入磁变测深数据的反演中．反演结果给出了测线地下100 km深度范围内的导电率结构,指出高导电率区位于武功、乾陵一带地下;并且讨论了该反演方法的主要特点,即能较好地反映异常导电率区域的深度和导电率分布,且反演中能允许大数量的反演模型参数．      

论中国地震台网测得的面波震级

文中使用了一个新的量规函数σ₃₀₀（Δ）,它是基于面波传播理论,考虑了面波的几何扩散、非弹性吸收、频散衰减以及面波优势周期随震中距的变化导出的,取全球平均Q=300。在震中距Δ=20°-160°的范围内与IASPEI推荐的面波震级的量规函数σ_IASPEI（Δ）一致,又推广了σ_IASPEI（Δ）的震中距使用范围,从1.0°-179°,以及可使用的面波周期范围（7.5-30s）,因此它更适合我国由中周期宽频带SK仪组成的台网情况。若采用此量规函数σ₃₀₀（Δ）和进行方位角改正,则可使我国台网测出的M_s与国际地震中心（ISC）测出的一致,震级误差降到0.19级,此误差与苏联台网的误差一样。     

二度体重、磁异常频谱与三度体重、磁异常频谱关系的理论问题

设F（x-ξ,y,z,η,ζ,α,…）为一在x方向无限薄的层的异常,若物体的异常可表示为: T（x,y）=integral from n=a to （a+2c）（F（x-ξ,y,z,η,ζ,α,…））dξ T（y）=integral from n=-∞ to ∞（F（x-ξ,y,z,η,ζ,α,…））dξ 则有 式中T（v）和T（u,v）分别为T（y）和T（x,y）的频谱,2c是物体在x方向的宽度,（x,y,z）和（ξ,η,ζ）分别为观测点和物体上任一点的坐标。 如果2c=dξ,则物体即为薄层本身;如果2c=dξ,薄层的面积无限缩小,则物体即为一质点;如果2c=dξ,薄层只有长度而无宽度,则物体即为一曲线。 对于形状复杂的物体,可以将物体分成若干简单而又符合所述条件的物体元,而有式中2c_m为第m个物体元在x方向的宽度,T_m（v）为第m个物体元当其沿x方向的长度为无限长时的异常的频谱。     

重、磁异常的走向与其振幅谱走向的关系

在文中利用坐标变换的方法,证明了重、磁异常的走向与其振幅谱的走向是垂直的。     

优化滤波方法及其在中国大陆布格重力异常数据处理中的应用

在优选延拓法的理论基础上,研究提出基于格林等效层概念和维纳滤波器的优化滤波法,用于对重力异常数据进行去噪和分离.与传统向上延拓法和优选延拓法相比,优化滤波法分离异常与延拓高度无关,不需要已知延拓高度,具有一定的优势.理论重力模型数据的去噪和异常分离试验表明优化滤波法有效,异常分离效果优于传统向上延拓法和带通滤波法.利用优化滤波法对中国大陆重力异常数据去噪和异常分离,得到有效的布格重力异常和区域重力异常.以中国大陆深地震探测推断的莫霍面深度信息为约束,对区域重力异常数据进行密度界面约束反演,得到中国大陆莫霍面深度分布.本文方法为中国大陆深部探测和区域构造研究提供一定的技术支撑.     

重、磁场三维傅里叶变换的若干问题

本文利用极限变换的概念,推导了质点的重力位及偶极子的磁位的三维傅里叶变换。在这个基础上,推导了任意形状物体的重、磁场的三维傅里叶变换的一般表达式。根据这些式子,讨论了文[1]中由三维谱求二维谱存在的问题,即对三维谱而言,不存在u²+v²+w²=0。     

分形滤波(S-A)技术在物探数据异常分解中的应用

分形滤波(S-A)技术能够有效地揭示地球物理场的空间相关性,对背景场以及异常场进行分解。应用分形滤波(S-A)技术分解内蒙古某银矿区地面高精度磁测异常,结果表明:对研究区高精度磁异常进行了有效分解,获得了研究区背景场以及异常场特征,并剖析所获得的异常场,分解异常8处,并据此结合地质情况提出研究区的成矿有利区。     

斜压大气中热带气旋运动特征的动力分析

从绝热、无摩擦、无环境场作用条件下的原始方程组出发,讨论了热带气旋的动力平衡特征,并分别导出了正、斜压大气中综合描述热带气旋中心强度变化的物理量(dI/dt)及其动力平衡关系式,进而分析了正、斜压大气中的热带气旋强度变化及移动的特征.结果表明:热带气旋的定常运动具有"相当地转风平衡”的特征,其等"压能”线和等位涡线均与流线重合;"β涡旋对”及由高(低)层水平风正(负)垂直切变产生的"切变涡旋对”之间的"通风气流”均使热带气旋的强度减弱,高(低)层涡度的负(正)垂直平流则增强热带气旋的强度;热带气旋的中心有向极、向高(低)层涡度垂直平流的梯(升)度方向移动的趋势,垂直运动及高、低层水平风速垂直切变的非对称结构对热带气旋的移动有影响;"γ效应”则是"通风气流”牵引热带气旋移动的动力机制之一.     

Applications of Savitzky-Golay Filter for Seismic Random Noise Reduction

This article utilizes Savitzky–Golay (SG) filter to eliminate seismic random noise. This is a novel method for seismic random noise reduction in which SG filter adopts piecewise weighted polynomial via leastsquares estimation. Therefore, effective smoothing is achieved in extracting the original signal from noise environment while retaining the shape of the signal as close as possible to the original one. Although there are lots of classical methods such as Wiener filtering and wavelet denoising applied to eliminate seismic random noise, the SG filter outperforms them in approximating the true signal. SG filter will obtain a good tradeoff in waveform smoothing and valid signal preservation under suitable conditions. These are the appropriate window size and the polynomial degree. Through examples from synthetic seismic signals and field seismic data, we demonstrate the good performance of SG filter by comparing it with the Wiener filtering and wavelet denoising methods.     

台湾海峡地震定点形变的异常分析

结合定点形变前兆异常的常用判定标准,用多种数学方法处理了厦门台与漳州台1991～1995年水平摆地倾斜资料,发现在9,16海峡地震之前,两个台站的资料都表现出不同程度的年变规律畸变、γ因子异常、一阶差分法异常、滤波法异常等多种地形变前兆指标异常。     

华北地区重力场变化特征与孕震模型研究

本文采用基于贝叶斯原理的重力平差方法对华北地区2009—2017年的16期陆地流动重力观测资料进行了统一处理,结合研究区内2009年以来4次M≥4.5 地震活动,分析区域重力场动态演化特征。基于华北地区均衡重力异常、布格重力异常、沉积层重力异常和莫霍面重力异常,研究华北地区三维密度结构特征,分析区域重力场变化与深部密度结构之间的关系。结果表明:研究区内近年来M≥4.5地震均位于区域差分重力场和累积重力场“零值”等值线附近、异常变化剧烈的梯度带及梯度带转弯部位,以及深部结构中理论均衡厚度与实际地壳厚度存在差异的非均衡区和三维密度结构的高低密度过渡带内。