《地球物理学报》稿约

正《地球物理学报》创刊于1948年,是中国科学院地质与地球物理研究所、中国地球物理学会联合主办的有关地球物理科学的综合性学术刊物.主要刊载固体地球物理、应用地球物理、地磁和空间物理、大气和海洋地球物理,以及与地球物理密切相关的交叉学科研究成果的高质量论文.中文版为月刊,在国内外发行,并由美国地球物理学家联合会(简称AGU)同步发行英文网络版.作者和读者对象主要为从事地球物理学、地球科学及其他相关学科的国内外科技工作者和大专院校师生.     

一阶弹性波方程数值模拟中的混合吸收边界条件

Liu和Sen（2010和2012）在地震波场数值模拟中提出一种混合吸收边界条件. 该方法具有计算量小、容易实现及吸收效果好等优点. 但现有的混合吸收边界条件是针对二阶位移-应力方程设计的,存在稳定性问题. 本文首先推导了两种速度-应力单程波方程：二阶Higdon单程波方程和一阶Higdon单程波方程. 进而提出基于一阶弹性波方程的混合吸收边界条件方法. 在内部区域和边界之间引入一个过渡区域,通过单程波与双程波方程平滑过渡来消除人工边界反射. 为了改善混合吸收边界条件的吸收效果和稳定性,我们采用了能同时吸收纵、横波反射的一阶单程波方程和与变量位置有关的加权系数. 为了验证混合吸收边界条件的有效性,将其与常规分裂完全匹配层（PML）方法进行了比较. 数值模拟结果表明,与PML边界条件相比,混合吸收边界条件在耗用更小计算时间和存储量的前提下,可以获得更好的吸收效果. 另外,本文提出的两种混合吸收边界条件中,混合一阶Higdon吸收边界条件具有更好的稳定性.     

基于声学全波形反演的油气藏地震成像方法

岩性油气藏在我国天然气勘探开发中占有非常重要的位置,其分布区域的成像是合理布设井位,提高钻井成功率的关键之一.本文首先基于地下介质的声学近似和波场回传理论,利用频率域单程声波方程延拓计算地震波场,进行全波形反演,获得地层密度和体积模量的定量成像,并依据油气藏物性特征和流体饱和多孔介质岩石物理模型,简要讨论了孔隙度和饱和度与密度及体积模量的关系,明确了地震油气藏成像新概念.在此基础上,定义了基于流体体积模量和孔隙度的成像函数,进行油气藏成像.理论模型计算表明该方法是可行的.通过对西部地区某气田二维地震数据处理,实现了致密砂岩气藏成像,钻井结果证实了气藏区域成像位置的准确性和方法的有效性.     

基于单能窄束伽马射线算法计算高分辨率屏蔽自然伽马探测器的三维测井响应

本文应用单能窄束伽马射线理论建立高分辨率屏蔽自然伽马探测器三维测井响应的数值模拟算法.首先利用单能窄束原理,用稳态扩散方程描述伽马光子密度的空间分布,根据扩散方程基本解、放射源空间分布、探测器位置和屏蔽情况,将屏蔽探测器上的总伽马通量表示成放射性地层中的有效探测区域上的体积分和晶体表面上有效接收面的面积分形式.并根据伽马射线传播路径和探测器屏蔽情况,给出有效探测区域的解析表达式,通过数值积分法计算任意复杂情况下探测器上的自然伽马通量,获得伽马测井响应的3D数值模拟算法.并通过数值模拟结果与模型井数据的对比验证了该算法的有效性.最后通过3D数值模拟算法系统研究考察晶体形状（圆柱形晶体、方形晶体）、晶体长度、仪器在井轴中的位置（居中屏蔽探测器、贴井壁屏蔽探测器）、以及测速等不同情况下自然伽马测井响应,设计出新型高分辨率自然伽马测井仪器.     

地震波初至走时的计算方法综述

在地震波场中,初至波到时信息由于初至震相可追踪、易识别性,在地震学领域占有重要的位置,广泛地应用于叠前偏移、叠前速度分析、地震走时层析成像及地震定位等.本文主要介绍了四类具有代表性的计算初至波走时的方法:(1)基于高频近似射线理论方法,如最短路径方法(SPM),及修正后的最短路径方法(MSPM);(2)基于程函方程的数值解方法,如有限差分方法(FD)、快速推进法(FMM)和快速扫描法(FSM);(3)基于惠更斯原理的波前构建法(WFC);(4)基于频率域波动方程数值解法(FWQ).最短路径方法计算精度较高,稳定性较好,但其需要采用更多的网格节点,因此计算效率低;程函方程数值解法无需计算射线路径,具有计算效率高、稳定性较好、易于实现等优势,但其计算精度较低,可以通过引入高阶差分格式得到提高;波前构建法计算精度高,稳定性好,但其需要在射线网格和规则网格之间做网格转换,因此计算效率较低;频率域波动方程方法能适应任意复杂介质,但其计算精度和计算效率较低.     

相关检测技术在低频交流地电阻率观测中的应用

目前我国的地震地电阻率观测主要采用单一极距四极对称装置的直流观测测方法,该观测越来越多的受到城市建设中出现的各种电磁干扰的影响,特别是地铁、轻轨等的干扰,使得观测数据严重超差.为此我国的地震工作者提出了交流地电阻率观测方法,该方法由于采用固定频率信号进行观测,可以抑制人文干扰.本文对低频交流地电阻率观测中所采用的相关检测技术原理及在观测系统研制中的应用进行了阐述,并采用数值模拟方法对相关检测技术及观测系统参数选择进行了详细分析,研究了在不同参数条件下的相关检测误差及影响因素,给出了观测系统主要性能指标的基本要求和设计参数,为观测系统研制提供基础和依据.     

合肥地震台地电阻率与气温、降水、地下水位关系的初步研究

经对地电阻率ρs观测值的分析研究,发现观测数据因受本身台址电性条件的限制,观测深度较浅,受降雨等气象因素的影响明显.运用合肥地电阻率北南、北西向2007-2013年观测资料月均值与温度、降水、地下水位进行相关分析,采用多元线性回归方法对地电阻率进行校正,以反映地电阻率的真实变化,达到提取地震异常的目的.     

数据填补及能量闭合对半干旱区生态系统年净碳交换的影响

基于半干旱区2种不同下垫面(草地和旱作农田)2005和2008年涡动相关法取得的通量资料,分析了数据填补、能量收支闭合率以及摩擦风速(u*)阈值等对生态系统年净碳交换的影响.通过加入4种不同长度的人工空缺(空缺长度从0.5 h~12 d),比较了平均日变化法(MDV)、边缘分布抽样法(MDS)和非线性回归法等6种填补方法的填补效果.结果表明,MDS的整体表现最好,特别是对长空缺的填补效果优于其他方法,估算的年NEE偏差在5 g C m-2 a-1以内.非线性回归法估算的夜间NEE具有较大的正偏差,表明非线性回归法估算的夜间生态系统呼吸偏高.4种非线性回归法估算的年NEE偏差在8.0~30.8g C m-2 a-1.由于在半干旱区土壤含水量是生态系统碳交换的重要限制因子,非线性回归法中综合考虑土壤温度和土壤含水量影响的Non_linear3和Non_linear4表现较好.MDV对白天NEE空缺的填补优于夜间,估算的年NEE偏差在-2.6~-13.4 g C m-2 a-1.总体上,数据填补的精确度受下垫面类型、空缺长度以及空缺出现时间(白天、晚上)影响.2个观测站点的能量收支闭合率在80%左右.能量收支闭合率受湍流强度影响显著;当夜间摩擦风速较低时,湍流混合不充分,能量收支闭合率也较低.生态系统在某个风向的累积通量印痕较大时,有效能量和湍流通量源区的不匹配造成这一风向上的能量收支闭合率也较低.通过假设能量收支不闭合全部由感热通量和潜热通量的低估引起,评估了能量闭合订正对生态系统CO2通量的影响.结果表明经过订正后的草地、农田的年净碳交换量平均增加近10 g C m-2 a-1.此外当u*阈值从0.1增加到0.2 m s-1,年净碳交换平均增加37.5 g C m-2 a-1,这表明u*阈值的设定对生态系统的年净碳交换影响较大.     

辅助微分方程完全匹配层在声波方程数值模拟中的应用

在地震波数值模拟中,为提高算法精度,需要使用高阶时间更新格式,而普通的非分裂完全匹配层（PML）吸收边界局限于低阶时间格式。辅助微分方程完全匹配层（ADE PML）是一种可以适应任意阶时间格式的非分裂完全匹配层技术,且可以直接应用复频移拉伸算子以提高PML在高角度入射时的效果。作者将ADE PML应用于声波方程四阶Runge Kutta时间格式的数值模拟中,对其吸收效能进行了检验。数值模拟表明,复频移ADE PML在高角度入射时表现优于非复频移ADE PML。另外,不同辅助变量更新格式的吸收效果存在微小差异,显格式下计算结果与解析解吻合较好。长时间能量衰减计算表明ADE PML可以稳定至2 × 10⁵时间步。     

西北干旱区荒漠生态系统通量贡献区模型研究

涡度相关仪通量值所代表的通量贡献区范围,对于通量观测塔的选址、仪器安装高度的确定以及通量观测数据的质量控制等具有重要的指导意义。利用通量贡献区模型对位于古尔班通古特沙漠试验场通量观测资料的空间代表性进行初步分析。结果表明：该荒漠区在大气稳定条件下90%的通量贡献区最远可以达到686.40 m,通量贡献函数最大点的位置在162.50 m;大气稳定时各风向的通量贡献区范围在生长末期均达到最大,生长初期和中期的源区变化因受到各风向风速和植被下垫面的影响而有差异;大气不稳定时不同生长时期各风向通量贡献区没有固定变化规律;通量源区大约有58.71%的信息来自于荒漠区通量观测塔西南至西北方,整个生长季生长末期通量贡献最多,所占比例为40.16%。由FSAM模型测得的通量贡献区范围可以较准确地反映荒漠生态系统下垫面的通量信息。