进入21世纪以来,全球范围内强烈地震频频发生,势已进入到一个新的地震活动时期.我国大陆内部、周边和台湾岛内一系列强烈地震(MS≥6.5)的发生,特别是2008年5月12日汶川MS8.0大地震的突发与其震前的异常平静表明:地表与上地壳均未见有明显的或确切能构成短、临强烈地震预测的浅表层活动过程,大地震由深部壳、幔物质运动,即深层动力过程所致.基...
利用GRACE卫星重力资料,计算了中国大陆及周边的卫星重力时变场和地表密度变化分布,获取了具有代表性的点位区域的每月重力变化时间序列.同时获得了WUSH、LHAS、KUNM、LUZH站相对于区域参考框架的GPS位移时间序列.卫星重力观测结果显示喜马拉雅弧形带的重力在2004年苏门答腊Mw9.3地震后快速下降, 2006~2008年尤为明显,西域地块西北边界带上震后...
2001年11月14日昆仑山口7.8级地震、2008年3月21日于田7.1级和5月12日的汶川7.9级地震,全部发生在青藏高原中部,构成新的地震活动组.昆仑山口和汶川地震分别位于巴颜喀拉活动地块的北部和东部边界,于田地震发生在该地块的西端.GPS的观测资料分析表明该块体整体向东运动.发生在块体不同部位的昆仑山口、于田和汶川3次大地震震源机制类型体现了巴...
为分析2008年5月12日四川汶川MS8.0级地震对周边断层地震活动的影响,本文首先基于Burgers体黏滞松弛模型计算汶川MS8.0级地震引起的库仑应力动态演化,分析认为2008年汶川MS8.0级地震在周边断层上引起的库仑应力显著增加的主要有四个断层段,分别为鲜水河断裂道孚-康定段、东昆仑断裂东段玛曲段、青川断裂和龙门山断裂南段.而且震后4年内黏滞松弛...
选取分布在北川-映秀中央断裂两侧断层距120 km以内的40个强震动台站的记录,对汶川地震近断层地震动竖向和水平向加速度峰值、速度峰值、竖向和水平向加速度反应谱及谱比值进行了统计分析.研究表明:(1)地震动加速度峰值有显著的上盘效应,经验衰减模型的结果表明,在距地表破裂3~60 km的范围内,龙门山发震断层上盘一侧竖向与水平向的加速度峰...
本文选用首都圈数字地震台网2003年9月~2005年12月记录的300多个远震事件的波形资料,采用分频带多道互相关方法得到三个不同频段的P波相对走时数据共18499个,计算了每个频段的走时灵敏度核,应用有限频率层析成像反演得到首都圈地区的上地幔三维P波速度结构模型.利用检测板估计了反演结果的分辨率,并与射线层析成像方法的结果进行了比较,说明了...
本文利用华北遥测地震台网和首都圈数字地震台网112个台站记录到的1993~2004年发生在首都圈地区3983次地震的P波绝对到时资料和相对到时资料,采用双差地震层析成像方法联合反演了京津唐地区地壳三维P波速度结构和震源参数.京津唐地区的三维P波速度结构图像在浅层上很好地反映了地表地质、地形的特征.在平原和凹陷的盆地处呈现P波低速速度异常,...
为了解决核爆地震自动识别中最优特征子集的选择问题,根据克隆选择原理,提出了一种过滤与封装相结合的特征选择方法.该方法融合了封装式与过滤式特征选择方法的优点,利用局部化的类别可分性判据来处理核爆地震样本的多峰分布问题,通过设定独立的记忆抗体能够保证最终结果是搜索过的最佳特征组合,并且可以处理设定和不设定最优特征子集维数两种情...
采用有限元方法模拟了俯冲带耦合作用对巽他弧及其邻区的影响.根据模拟结果,对比GPS、地震和地质学观测数据,定量分析了苏门答腊及其周边地区的应变强度和主应变方向的分布特征,据此探讨了该区构造特征、地震发生模式与耦合面积之间的关系.模型由具有黏弹性性质的岩石圈和软流圈上地幔组成,其中岩石圈包括了大陆岩石圈和大洋岩石圈以及俯冲至上...
反射地震方法一般用于地下地质结构的探测.本文采用南海东北部的一条地震测线记录,通过对水体反射信号的处理,获取了该测线上海水的声速剖面,并获得了海水水体结构的直观图像.在反射地震剖面上,可以看出南海东北部上层海水结构较为稳定,没有明显的水体侵入痕迹,但中、深层海水中内波较为发育,这与物理海洋学上的研究结果是一致的.本文通过设计合...
鲜水河断裂是青藏高原东南缘的一条北西向大型左旋走滑断裂,其南东段逐渐向南偏转,并与近南北向的安宁河断裂相接,在两个断裂相接处西侧耸立着海拔7556 m高的贡嘎山.磷灰石裂变径迹(AFT)测试可知,贡嘎山及其邻区12个样品的年龄分布在0.2±0.1 Ma~2.7±0.7 Ma之间,平均径迹长度在13.64~15.19 μm之间,表明贡嘎山及其邻区第四纪时期一直处于快速...
年龄-高程法是地质热年代学中利用垂向分布的样品高程和年龄关系计算山体剥露速率的技术,其主要优势在于避免了未知古地温梯度的约束,被广泛应用于造山带研究中.但在实际应用中,由于样品采集剖面无法完全满足垂直的要求,会引入封闭温度面受地形起伏影响所带来的计算误差.本文利用稳态地温场数值模拟,在理论三角函数地形条件下,对不同采样剖面和...
以塔里木前陆盆地地质背景为模型,采用三维任意差分精细积分方法并行算法实现了库车地区三维正演模拟.三维任意差分精细积分方法通过时间域采用局部积分半解析方法求得波动方程的递推算子,与常规的差分法相比计算精度有较大提高;文中计算稳定性根据实际算例进行分析,采用稳定因子约束,得到较好的稳定性;边界条件采用改进的自适应吸收边界,并通过...
本文针对时间推移地震本身包含不同时期的两次或者两次以上的勘探反问题,构造了一种快速有效的反演方法--多尺度全变分法.通过引入全变分正则化来代替传统的Tikhonov正则化,针对待反演参数不连续的情况,提高了算法精度.为了提高计算效率,引入了多尺度方法,从而构造了多尺度全变分方法.在数值模拟中,针对一个时间推移地震反演问题对多尺度-Tikhon...
本文基于拟微分算子理论和李代数积分法,根据程函方程和波场坐标变换,提出一种新的适于横向变速介质Kirchhoff叠前深度偏移的地震波走时算法.该算法与Kirchhoff叠前时间偏移所用李代数时间积分表达相比,差异在于增加了波数一次项,且二次项的系数在求积时亦需进行修正.针对单平方根算子象征、李代数积分、指数映射和走时多项式的求解而言,皆需对...
反射地震勘探中的偏移成像技术是获取地下介质构造形态最有效的手段之一.在叠前深度域偏移方法中,目前工业界采用的方法包括基于射线理论的波动方程积分解法和基于波动理论的微分波动方程单程波解法,这两类方法难以处理地震波横向速度变化剧烈的高陡倾角构造成像问题.近年来勘探地震学研究领域发展起来的叠前逆时偏移采用了双程波求解微分波动方...
相速度和偏振方向是研究地震波传播规律和描述介质特性的重要参数,在理论研究和实际应用中有重要作用.本文假定倾斜横向各向同性(TTI)介质对称轴位于观测坐标系XOZ面内,在此观测坐标系下直接推导了TTI介质弹性波相速度和偏振方向的解析表达式,再进一步利用Thomsen弱各向异性理论,推导了弱各向异性近似条件下弹性波相速度以及qP波和qSV波偏振方...
在前人研究基础上,对Groom-Bailey(GB)张量分解畸变因子和区域阻抗的求解方法进行了改进.首先,通过Swift旋转与GB分解的扭变和剪切矩阵的求逆变换,利用变换后区域阻抗主对角元素为0的条件获得关于扭变因子和剪切因子的超定方程组,采用模拟退火全局优化算法进行求解.其次,由得到的扭变因子和剪切因子,结合Swift旋转确定的走向角和区域阻抗元素...
本文采用电磁场的磁矢量位和电标量势,将Maxwell方程组化为位势的类似于Helmholtz型方程,并引入罚项及稳定化方法克服了电磁三维有限元算法中的伪解及数值不稳定性;采用人工边界把计算区域局域化,将均匀半空间中水平电偶极子源产生的位势值作为人工边界上的第一类边界条件以表示源的作用,减少了实际的计算区域.理论模型和复杂模型的计算结果均表...
为了探索我国海域油气和水合物等高阻目标体CSEM勘探的可行性和方法技术,本文研究了在海水中水平电性源激励下有限水深海洋地电模型的频率域电磁响应,为进一步的1D和3D仿真计算奠定了理论基础.在推导电磁响应公式时,首先给出了各层介质的Lorentz势,然后根据Coulomb势与Lorentz势的关系,得到了各层介质的Coulomb势.各层介质中的电磁场均可以由Lor...
根据孔隙介质的动电耦合理论设计了一套岩芯渗透率测量系统.实验采用交流锁相放大技术,在低频12~42 Hz范围内完成了砂岩岩样流动电势和电渗实验,得到了流动电势系数KS和电渗压力系数KE,进而计算出岩样动电渗透率,对于中、高渗透率岩样,测量得出的动电渗透率与常规气测渗透率差异较小,两者具有很好的相关性.实验表明,动电测量可作为岩样渗透率测...
阐述了熵平均法对地面γ能谱测量局部干扰的消除的基本原理和方法,应用此方法对浅覆盖区地质填图的γ能谱测量剖面数据及应用γ能谱测量确定古城墙位置的平面测量数据进行了处理,处理结果表明该方法对于消除地表植被及局部土壤元素富集随机干扰、突出弱异常具有明显的效果,使地质填图中进行岩性划分界线明显,考古中古城墙表现的弱异常突出,边界清晰...
随钻地震参考信号处理的关键是从众多噪音成分中提取较弱的钻头源信号.钻头随机激励源信号由于其频率范围宽,持续时间短,易受钻柱和其他机器噪音的影响.倒谱分析方法是一种非线性滤波技术,该技术不但可以将时域卷积信号转变为倒时域信号的相加,而且通过窗函数的选择,还可消除结构混响,提取源信号.本文根据钻柱模拟实验数据,利用了倒谱滤波技术进...
在黄海和东海采样测定了水体中的镭同位素分布,用平流扩散模型描述镭同位素分布,最小二乘方法计算了垂直涡动扩散系数和上升流或下降流流速.结果给出北黄海中部、南黄海中部、浙江沿岸和台湾北部海域存在上升流,流速分别为0.46×10-3cm·s-1、0.17×10-3~1.39×10-3cm·s-1、2.02×10-3~3.04×10-3 cm·s-1和1.06×10-3~2.51×10-3 cm·s-1.北黄海中部和...
加权近似解析离散化(WNAD) 方法是近年发展的一种在粗网格步长条件下能有效压制数值频散的数值模拟技术. 在地震勘探的实际应用中, 不是所有情况都适合使用空间大网格步长. 为适应波场模拟的实际需要, 本文给出了求解波动方程的非一致网格上的WNAD算法. 这种方法在低速区、介质复杂区域使用细网格, 在其他区域采用粗网格计算. 在网格过渡区域,...
本文基于声波方程的哈密尔顿系统,构造了一种新的保辛数值格式,简称NSPRK方法.该方法在时间上采用二阶辛可分Runge-Kutta方法,空间上采用近似解析离散算子进行离散逼近.针对本文发展的新方法,我们给出了NSPRK方法在一维和二维情况下的稳定性条件、一维数值频散关系以及二维数值误差,并在计算效率方面与传统辛格式和四阶LWC方法进行了比较.最后,...
支持向量机(Support Vector Machine: SVM)一直作为机器学习方法在统计学习理论基础上被研究和发展,本文从信号与系统的角度出发,证明了平移不变核最小二乘支持向量机(Least Squares SVM: LS-SVM)是一个线性时不变系统.以Ricker子波核为例,探讨了不同参数对最小二乘支持向量回归(Least Squares Support Vector Regression: LS-SVR)滤波...
2010年7月12~15日,第二届“反演问题计算方法及其应用国际会议”(The Second International Workshop on Computational Inverse Problems and Applications)在中国北京召开.会议由中国科学院地质与地球物理研究所油气资源研究重点实验室主办.会议主题侧重于不适定反演问题、正则化理论及最新进展、最优化理论方法、非标准的正则化方法、...
《地球物理学报》创刊于1948年,是中国科学院地质与地球物理研究所、中国地球物理学会联合主办的有关地球物理科学的综合性学术刊物.主要刊载固体地球物理、应用地球物理、地磁和空间物理、大气和海洋地球物理,以及与地球物理密切相关的交叉学科研究成果的高质量论文.中文版为月刊,