唐裕
- 作品数：16 被引量：51H指数：5
- 供职机构：吉林大学地球探测科学与技术学院更多>>
- 发文基金：江苏省自然科学基金国家高技术研究发展计划安徽省国土资源厅科技项目更多>>
- 相关领域：天文地球更多>>
- 基于全局光滑约束的地磁测深C-响应估计被引量：2
- 2019年
- 地磁测深C-响应,包含着地球内部结构的导电信息,因此获得高质量的C-响应估计对于揭示地球内部准确的电性结构至关重要,为此本文提出一种计算地磁测深C-响应的新思路.不同于以往单个频点估计办法,新方法基于相邻频率C-响应连续光滑的特性同时估计所有频率的C-响应.首先,根据求取C-响应的Z/H方法,由观测的Hr和Hθ构造估计C-响应的线性方程组;为了增强方程组求解的稳定性,克服噪声影响,要求估计的C-响应满足光滑性条件,从而将无约束的C-响应估计转化为阻尼最小二乘理论的C-响应估计优化问题.合成数据和实际台站数据的测试结果表明,光滑约束优化技术能比传统的逐频求取方法得到更加合理和连续的响应曲线;虽然正则化参数的选择对C-响应估计的结果有明显影响,但通过L-曲线和V-曲线确定的最优正则化参数可以在保证数据真实性的前提下获得最佳的C-响应估计.基于全局光滑约束获得的C-响应为提高反演结果的可靠性奠定基础,为利用更多台站的C-响应获得高分辨力的地球深部电性结构提供了技术支持.
- 张艳辉翁爱华李世文李建平杨悦唐裕
- 关键词：地磁测深地球内部
- 全球地磁感应测深数据三维反演被引量：4
- 2019年
- 全球地磁感应测深能获得地幔转换带及下地幔上部的导电结构.但目前稀疏的地磁台站分布及部分台站的观测数据稳定性较差,影响了三维反演对地下电性结构的分辨力和反演可靠性.为此,区别于传统的L2-范数反演方法,本文提出并实现了基于L1-范数的地磁测深响应三维反演技术.在反演中,利用L1-范数度量数据预测误差,降低"飞点"数据的影响,将相关系数较小的C-响应估计也纳入反演数据中.三维正演模拟采用球坐标系下的交错网格有限差分法,反演采用有限内存拟牛顿法.文中利用指数概率密度分布函数构造非高斯噪声的合成数据,并采用棋盘模型对反演方法的可靠性进行了验证.之后,我们将本文提出的三维反演方法用于全球129个地磁观测台站的C-响应数据反演,结果表明在地幔转换带深部,中国东北地区为高导电异常,南欧和北非则均为高阻;夏威夷在900km以下为高导;菲律宾海及以东地区在转换带表现为明显的高阻,这些结果与前人研究结果一致.由于采用了更多的台站数据,我们的反演结果还发现一些新的异常:南美洲南端,转换带表现为明显的高导;澳大利亚东南部,地幔转换带深部,也存在一个明显的高导异常,这些异常分布和地震层析成像的低速区一致.因此,L1-范数三维反演能够充分利用全球C-响应数据信息,提高地磁测深对地球深部电性结构的分辨能力,更好的研究全球地幔电性结构.
- 李世文翁爱华张艳辉李建平杨悦唐裕邹宗霖李春成
- 关键词：地磁测深三维反演地幔结构
- 流体诱发松原地震的大地电磁图像
- 松原地区位于松辽盆地的腹地.自2006年以来,地震活动频发,地震震中集中分布在松原市西南的查干花震区和北侧的宁江震区.地震多以震群形式出现,出现多次5级以上的地震.地震机制解表明,部分地震具有剪切分量.因此,部分地震学家...
- 翁爱华唐裕李亚彬杨悦
- 关键词：发震机制流体作用
- 中国异构地磁数据C-响应估计
- 2019年
- 中国分布着世界上较为密集的地磁台站,为利用地磁测深获得转换带及下地幔上部高分辨率电性结构奠定了基础,但相应的C-响应计算还需要进行系统研究。文中针对中国地磁数据可能存在的噪声干扰、记录长度长短不一、数据记录间断以及数据观测种类不同(绝对和相对)等异构特点,采用如下技术来提高C-响应曲线的稳定性:1)基于BIRRP软件包获得地磁各分量时间序列的频率响应,使用整体光滑技术来压制数据噪声;2)对于相对观测台站的数据,采用近似估计法对水平分量进行分解,进而估计C-响应;3)讨论非连续数据和短记录数据对C-响应估计的影响;4)对于沿海台站,可使用比值法消除海洋效应对C-响应函数的影响。使用以上技术手段对中国100余个地磁台站的数据进行处理后,最终得到42个台站的稳定C-响应函数,其中周期为1.3~113.7d的响应个数为24个;周期为1.3~42.6d的响应个数为18个。采用文中的技术手段能够很好地处理异构数据,获得更多、更稳定的C-响应函数,为中国高分辨率地磁测深反演研究提供了更为丰富的基础数据。
- 张艳辉李世文翁爱华张素琴杨悦李建平唐裕
- 关键词：地磁测深海洋效应
- 松原地震与流体作用联系的大地电磁证据被引量：8
- 2021年
- 松原地震区是目前东北地区危险性最大的地震区.特别是自2006年以来,地震以震群形式持续发生,最大震级达到了5.8级.关于松原地震根本原因,目前尚存在较大争议.文章利用松原地区面积性大地电磁探测结果,试图从电性结构角度揭示其原因.对观测到的全阻抗张量数据采用三维反演技术,获得了地震区深部三维电阻率结构.结合地震定位结果,电阻率模型显示,纵贯松原地震区存在北东向的隐伏断层,推测为扶余-肇东断裂.研究还发现,在地震区岩石圈深部存在明显的低阻异常,该异常向上穿过高阻岩石圈地幔后滞留在地壳底部.岩石物理分析表明,地壳深部的低阻异常由含水熔融物质组成.该含水熔融物质在温度和临界压力共同作用下在下地壳发生岩浆侵位释放出流体;它们在弱化扶余-肇东断裂强度的同时,还增加断层的孔隙压力,降低断层的抗剪强度;区域性压应力的剪切作用(区域性近东西向压应力沿北东向断层的走向分量)和近地表可能的扰动共同作用引起断裂失稳,诱发松原地震.来自深部持续不断的流体补给,导致松原地震以震群形式出现.文章认为,松原地震是板块俯冲、既有断裂、区域应力环境和深部流体共同作用结果,并可能由近地表扰动触发.松原地震根本上属于深部流体(水)作用诱发的地震.
- 唐裕翁爱华杨悦李世文牛建军张艳辉李亚彬李建平
- 关键词：大地电磁测深三维反演流体诱发地震
- 基于可控源电磁法阻抗信息的有限内存拟牛顿法三维反演被引量：7
- 2019年
- 本文研究了利用阻抗信息进行可控源电磁勘探有限内存拟牛顿法三维反演的技术。首先用理论模型来验证有限内存拟牛顿法反演的准确性和可行性。观测参数为复阻抗Zxy分量,采用交错网格有限差分方法计算模型响应,反演采用有限内存拟牛顿法。数值模拟结果表明:有限内存拟牛顿法反演迭代速度较快,每4min迭代一次,拟合差由146.00下降到1.78,收敛稳定;异常体的位置与理论模型吻合较好,有效地验证了有限内存拟牛顿法可控源三维反演的正确性。为了进一步验证该方法的实用性,将其应用到隐伏钼矿可控源电磁勘探工作中。工区的反演结果显示:在工区北西段深部存在高阻异常,其上为低阻异常。截取过钻孔的3号测线发现,-500～-100m的位置表现为低阻,东侧存在向上涌起的高阻。推测此低阻为矿化蚀变带,延伸较深。该异常与钻探资料揭示的钼矿脉一致,证明了反演结果的准确性。因此,利用可控源观测得到的阻抗信息进行有限内存拟牛顿法三维反演,可以获得可靠的三维电阻率分布。
- 杨悦翁爱华张艳辉李世文李建平唐裕
- 关键词：可控源电磁法三维反演
- 地磁测深的三维L1-范数反演
- 反演中,通常假设数据和模型参数具有高斯概率分布特征,通过求解L2-范数的目标函数得到一个光滑模型并使其响应与反演数据相拟合.但观测数据中很难避免坏点的出现,地下介质也会出现电性结构的突变,此时将不再满足高斯概率分布的假设...
- 李世文翁爱华李建平张艳辉唐裕杨悦
- 关键词：地磁测深三维反演
- CSAMT一维全频视电阻率反演被引量：2
- 2017年
- 本文采用全频率视电阻率反演,使用2^(-3)~2^(13)Hz频率,并且不做近区和过渡区校正,直接用卡尼亚视电阻率作为反演参数。正演采用虚界面法计算有限长导线在水平层状介质的电磁场分量,反演采用有限内存拟牛顿法。数值模拟采用三层模型和实测数据进行反演,得到了符合地下电性变化的反演模型结果,证实了CSAMT一维全频视电阻率反演具有可实施性。
- 苌云张艳辉翁爱华李斯睿李建平唐裕
- 关键词：CSAMT反演
- 流体诱发松原地震的大地电磁图像
- 2019年
- 松原地区位于松辽盆地的腹地。自2006年以来,地震活动频发,地震震中集中分布在松原市西南的查干花震区和北侧的宁江震区。地震多以震群形式出现,出现多次5级以上的地震。地震机制解表明,部分地震具有剪切分量。因此,部分地震学家认为松原地震为构造地震。然而,活动断层调查没有发现与地震相关的活动断层。
- 翁爱华唐裕李亚彬杨悦
- 关键词：地震震中大地电磁图像活动断层松辽盆地
- 基于球坐标系下有限差分的地磁测深三维正演被引量：5
- 2018年
- 为了计算全球尺度电磁感应的响应,本文介绍地磁测深频率域三维正演。正演算法采用球坐标系下的交错网格有限差分方法,从Maxwell方程的积分形式出发,采用PARDISO对离散后的方程组求解,避免了迭代求解的散度校正。为了验证本文结果的正确性和精度,与前人的有限元和有限差分方法进行了对比,一维层状模型的三维交错网格有限差分数值结果和解析解相对误差小于5%,双半球模型的计算结果与前人的计算结果完全吻合。三维"棋盘模型"计算表明磁场分量对异常体的大小和位置具有很好的分辨能力。
- 李建平翁爱华李世文李大俊李斯睿杨悦唐裕张艳辉
- 关键词：地磁测深球坐标三维正演有限差分
