CN115343751A - 一种活动地块边界带断层震级确定方法及系统 - Google Patents

一种活动地块边界带断层震级确定方法及系统 Download PDF

Info

Publication number
CN115343751A
CN115343751A CN202210932909.5A CN202210932909A CN115343751A CN 115343751 A CN115343751 A CN 115343751A CN 202210932909 A CN202210932909 A CN 202210932909A CN 115343751 A CN115343751 A CN 115343751A
Authority
CN
China
Prior art keywords
fault
earthquake
time
magnitude
displacement
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202210932909.5A
Other languages
English (en)
Inventor
刘琦
邵志刚
冯蔚
刘晓霞
王芃
熊仁伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
INSTITUTE OF EARTHQUAKE SCIENCE CHINA EARTHQUAKE ADMINISTRATION
Original Assignee
INSTITUTE OF EARTHQUAKE SCIENCE CHINA EARTHQUAKE ADMINISTRATION
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by INSTITUTE OF EARTHQUAKE SCIENCE CHINA EARTHQUAKE ADMINISTRATION filed Critical INSTITUTE OF EARTHQUAKE SCIENCE CHINA EARTHQUAKE ADMINISTRATION
Priority to CN202210932909.5A priority Critical patent/CN115343751A/zh
Publication of CN115343751A publication Critical patent/CN115343751A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G01V1/01—
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/22—Transmitting seismic signals to recording or processing apparatus
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/282—Application of seismic models, synthetic seismograms
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy

Abstract

本发明适用于震级预估技术领域，尤其涉及一种活动地块边界带断层震级确定方法及系统，所述方法包括：通过反演给出我国大陆活动地块边界带主要断层段现今断层位移加载速率；根据断层段末次强震发生时间的记录以及区域内地震目录完备的起始时间获得各断层段末次强震离逝时间。本发明根据断层位移加载速率、强震离逝时间等的差异，计算获得不同断层段的震间累积位移，用于基于位错可预测的强震复发模型，利用震级与位错经验公式得到震级预测结果，相关结果可以直接服务于防震减灾工作，为地震重点监视防御区等的确定和更新提供重要依据，指导相关部门和机构针对重点区域提前采取相关措施，有效降低未来强震发生可能造成的灾害损失。

Description

一种活动地块边界带断层震级确定方法及系统
技术领域
本发明属于震级预估技术领域，尤其涉及一种活动地块边界带断层震级确定方法及系统。
背景技术
地震灾害是人类面临的最严重的自然灾害之一，地震多、强度大、分布广、灾害重是我国的基本国情，而随着我国城市化的快速推进和国民经济的发展，地震灾害损失将有可能呈指数上升之势，因此未来严峻的地震灾害形势和社会经济持续稳定发展的需求都对我国地震预测工作提出了更高的要求。
关于地震预测，最理想的情况是给出明确的“时、空、强”三要素，而且最好能以定量的方式表述预测结果。然而基于目前的观测水平和理论发展，很难同时给出上述几种要素的准确预测结果。因此，较为可行的做法通常是在适当放宽单个要素要求的基础上，尝试同时给出其它要素的预测结果，这样就使得地震预测问题的难度大大降低，具备了更强的可操作性。2008年汶川8.0级地震、2011年日本9.0级地震等都造成了远超预期的人员伤亡和财产损失，究其原因主要在于震前对于未来震级的错误估计。因此开展相对准确的地点+震级预测，即预测某地未来可能会发生多大震级的地震具有重要的科学意义和实用价值。
我国大陆主要构造变形和浅源强震多发生在活动地块边界带上，有历史记载以来，我国大陆所有8级以上、超过85％的7级以上浅源地震都发生在活动地块边界带上。因此，以我国大陆活动地块边界带主要活动断层为主，基于统一的模型和方法框架，给出相对完备的震级预测结果，对于我国大陆整体的防震减灾工作具有重要意义。当前我国大陆活动地块边界带上不同断层的基础资料丰富程度存在显著差异，因此需要发明一套可兼容不同资料丰富程度的震级确定方法。
发明内容
本发明实施例的目的在于提供一种活动地块边界带断层震级确定方法，旨在解决当前我国大陆活动地块边界带上不同断层的基础资料丰富程度存在显著差异所导致的兼容性问题，基于统一的模型和方法框架，给出相对完备的震级预测结果，为我国大陆整体防震减灾工作提供支撑。
本发明实施例是这样实现的，一种活动地块边界带断层震级确定方法，所述方法包括：
通过反演给出我国大陆活动地块边界带主要断层段现今断层位移加载速率；
根据断层段末次强震发生时间的记录以及区域内地震目录完备的起始时间获得各断层段末次强震离逝时间；
根据各断层段的断层位移加载速率和末次强震离逝时间，计算各断层段累积位移量；
基于位错可预测的强震复发模型，依据震级与位错的经验公式确定主要断层段的发震震级。
优选的，所述通过反演给出我国大陆活动地块边界带主要断层段现今断层位移加载速率的步骤中：
基于各断层段几何产状参数和闭锁深度参数构建用于反演计算的断层模型；
将大地测量观测的地壳运动速率作为反演约束进行引入，并且将速度误差超过第一预设值以及距离断层距离小于第二预设值的点位数据舍弃；
基于弹性回跳和块体运动理论及断层连接元模型，反演计算断层位移加载速率，主要公式为：
x＝[ATC-1A+BTD-1B]-1ATC-1y，
其中，x为断层错动分量构成的向量，y是大地测量观测的东向和北向速度分量构成的速度场向量，A为各个断层错动分量对某一测站速度分量的贡献构成的传播矩阵，B为包含了结点处位移连续性条件和拉张分量限制条件的断层系统信息的传播矩阵，C为各个测站速度分量的方差-协方差构成的矩阵，D为各结点位移连续性误差的方差矩阵。
优选的，所述根据断层段末次强震发生时间的记录以及区域内地震目录完备的起始时间获得各断层段末次强震离逝时间的步骤中，对于存在末次强震发生时间记录的断层段，直接使用记录计算获得末次强震离逝时间。
优选的，所述根据断层段末次强震发生时间的记录以及区域内地震目录完备的起始时间获得各断层段末次强震离逝时间的步骤中，对于缺乏末次强震发生时间记录的断层段，进行末次强震离逝时间期望值的估计。
优选的，进行末次强震离逝时间期望值的估计的过程，具体包括：
依据区域研究结果确定各断层段地震目录完备的起始时间；
根据公式
Figure BDA0003782402270000031
估计末次强震离逝时间的期望值，其中ξ(τ|τ≥TH)为末次强震离逝时间的期望值，TH为地震目录完备的起始时间，F为累积概率密度函数，分别采用对数正态和布朗过程时间两种模型进行计算；
以计算结果平均值作为最终结果。
优选的，各断层段累积位移量由公式S＝v*T计算，其中S为各断层段累积位移量，v为各断层段的断层位移加载速率，T为各断层段末次强震离逝时间。
优选的，用于基于位错可预测的强震复发模型确定震级，该模型假设每次地震会将应力彻底释放至稳定值，强震离逝时间越长，断层加载速率越高，则累积位移量越大，预期发生的地震震级就越高，其计算结果是时间依赖的。
优选的，所述震级与位错的经验公式为：
M＝a+b log10(AD)，
其中，M为预测震级，AD为断层段累积位移量；a和b为线性回归系数，a取6.93±0.05，b取0.82±0.1。
本发明实施例的另一目的在于提供一种活动地块边界带断层震级确定系统，所述系统包括：
断层位移加载速率反演模块，用于通过反演给出我国大陆活动地块边界带主要断层段现今断层位移加载速率；
末次强震离逝时间确定模块，用于根据断层段末次强震发生时间的记录以及区域内地震目录完备的起始时间获得各断层段末次强震离逝时间；
累积位移量计算模块，用于根据各断层段的断层位移加载速率和末次强震离逝时间，计算各断层段累积位移量；
震级确定模块，用于基于位错可预测的强震复发模型，依据震级与位错的经验公式确定主要断层段的发震震级。
优选的，所述断层位移加载速率反演模块包括：
断层模型构建单元，用于基于各断层段几何产状参数和闭锁深度参数构建用于反演计算的断层模型。
观测资料约束单元，用于将大地测量观测的地壳运动速率作为反演约束进行引入，并且将速度误差超过第一预设值以及距离断层距离小于第二预设值的点位数据舍弃。
反演计算单元，用于基于弹性回跳和块体运动理论及断层连接元模型，反演计算断层位移加载速率，主要公式为：
x＝[ATC-1A+BTD-1B]-1ATC-1y，
其中，x为断层错动分量构成的向量，y是大地测量观测的东向和北向速度分量构成的速度场向量，A为各个断层错动分量对某一测站速度分量的贡献构成的传播矩阵，B为包含了结点处位移连续性条件和拉张分量限制条件的断层系统信息的传播矩阵，C为各个测站速度分量的方差-协方差构成的矩阵，D为各结点位移连续性误差的方差矩阵。
优选的，所述末次强震离逝时间确定模块包括：
直接记录单元，用于对存在末次强震发生时间记录的断层段，直接使用记录计算获得末次强震离逝时间。
期望估计单元，用于对缺乏末次强震发生时间记录的断层段，进行末次强震离逝时间期望值的估计。
优选的，所述期望估计单元包括：
目录完备时间子单元，用于依据区域研究结果确定各断层段地震目录完备的起始时间。
期望值单独求解子单元，用于根据公式
Figure BDA0003782402270000051
估计末次强震离逝时间的期望值，其中ξ(τ|τ≥TH)为末次强震离逝时间的期望值，TH为地震目录完备的起始时间，F为累积概率密度函数，分别采用对数正态和布朗过程时间两种模型进行计算。
期望值平均子单元，用于将基于对数正态和布朗时间过程两种模型的强震离逝时间期望值结果进行平均，以该平均值作为最终结果。
与现有技术相比，本发明的有益效果是：
本发明实施例提供的一种活动地块边界带断层震级确定方法，对我国大陆活动地块边界带主要活动断层，根据断层位移加载速率、强震离逝时间等的差异，计算获得不同断层段的震间累积位移，并利用经验公式得到震级预测结果，克服了不同断层段基础资料丰富程度显著差异所导致的兼容性问题。相关结果可以直接服务于防震减灾工作，为地震重点监视防御区等的确定和更新提供重要依据，指导相关部门和机构针对重点区域提前采取相关措施，以提升建筑抗震、应急准备与救援等各方面能力，有效降低未来强震发生可能造成的灾害损失。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例。
图1为本发明实施例提供的方法的流程图；
图2为本发明实施例提供的方法中断层位移加载速率反演的流程图；
图3为本发明实施例提供的方法中末次强震离逝时间确定的流程图；
图4为本发明实施例提供的方法中末次强震离逝时间期望值估计的流程图；
图5为本发明实施例提供的系统的架构图；
图6为本发明实施例提供的系统中断层位移加载速率反演模块的架构图；
图7为本发明实施例提供的系统中末次强震离逝时间确定模块的架构图；
图8为本发明实施例提供的系统中期望估计单元的架构图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
可以理解，本申请所使用的术语“第一”、“第二”等可在本文中用于描述各种元件，但除非特别说明，这些元件不受这些术语限制。这些术语仅用于将第一个元件与另一个元件区分。举例来说，在不脱离本申请的范围的情况下，可以将第一xx脚本称为第二xx脚本，且类似地，可将第二xx脚本称为第一xx脚本。
图1示出了本发明实施例提供的方法的流程图。
具体的，一种活动地块边界带断层震级确定方法，所述方法具体包括以下步骤：
步骤S100，通过反演给出我国大陆活动地块边界带主要断层段现今断层位移加载速率。
具体的，图2示出了本发明实施例提供的方法中断层位移加载速率反演的流程图。
其中，在本发明提供的优选实施方式中，所述断层位移加载速率反演具体包括以下步骤：
步骤S1001，基于各断层段几何产状参数和闭锁深度参数构建用于反演计算的断层模型。
步骤S1002，将大地测量观测的地壳运动速率作为反演约束进行引入，并且将速度误差超过第一预设值以及距离断层距离小于第二预设值的点位数据舍弃。
步骤S1003，基于弹性回跳和块体运动理论及断层连接元模型，反演计算断层位移加载速率，主要公式为：
x＝[ATC-1A+BTD-1B]-1ATC-1y，
其中，x为断层错动分量构成的向量，y是大地测量观测的东向和北向速度分量构成的速度场向量，A为各个断层错动分量对某一测站速度分量的贡献构成的传播矩阵，B为包含了结点处位移连续性条件和拉张分量限制条件的断层系统信息的传播矩阵，C为各个测站速度分量的方差-协方差构成的矩阵，D为各结点位移连续性误差的方差矩阵。
进一步的，所述一种活动地块边界带断层震级确定方法还包括以下步骤：
步骤S200，根据断层段末次强震发生时间的记录以及区域内地震目录完备的起始时间获得各断层段末次强震离逝时间。
具体的，图3示出了本发明实施例提供的方法中末次强震离逝时间确定的流程图。
其中，在本发明提供的优选实施方式中，所述末次强震离逝时间确定具体包括以下步骤：
步骤S2001，对存在末次强震发生时间记录的断层段，直接使用记录计算获得末次强震离逝时间。
步骤S2002，对缺乏末次强震发生时间记录的断层段，进行末次强震离逝时间期望值的估计。
具体的，图4示出了本发明实施例提供的方法中末次强震离逝时间期望值估计的流程图。
其中，在本发明提供的优选实施方式中，所述末次强震离逝时间期望值估计具体包括以下步骤：
步骤S20021，依据区域研究结果确定各断层段地震目录完备的起始时间。
步骤S20022，根据公式
Figure BDA0003782402270000091
估计末次强震离逝时间的期望值，其中ξ(τ|τ≥TH)为末次强震离逝时间的期望值，TH为地震目录完备的起始时间，F为累积概率密度函数，分别采用对数正态和布朗过程时间两种模型进行计算。
步骤S20023，将基于对数正态和布朗时间过程两种模型的强震离逝时间期望值结果进行平均，以该平均值作为最终结果。
进一步的，所述一种活动地块边界带断层震级确定方法还包括以下步骤：
步骤S300，用于根据各断层段的断层位移加载速率和末次强震离逝时间，计算各断层段累积位移量。
在本步骤中，具体的，各断层段累积位移量由公式S＝v*T计算，其中S为各断层段累积位移量，v为各断层段的断层位移加载速率，T为各断层段末次强震离逝时间。
进一步的，所述一种活动地块边界带断层震级确定方法还包括以下步骤：
步骤S400，用于基于位错可预测的强震复发模型，依据震级与位错的经验公式确定主要断层段的发震震级。
在本步骤中，具体的，用于基于位错可预测的强震复发模型确定震级，该模型假设每次地震会将应力彻底释放至稳定值，强震离逝时间越长，断层加载速率越高，则累积位移量越大，预期发生的地震震级就越高，其计算结果是时间依赖的。依据现有技术中震级与位错的经验公式确定主要断层段的发震震级：
M＝a+b log10(AD)，
其中，M为预测震级，AD为断层段累积位移量；a和b为线性回归系数，a取6.93±0.05，b取0.82±0.1。
进一步的，图5示出了本发明实施例提供的系统的应用架构图。
其中，在本发明提供的又一个优选实施方式中，一种活动地块边界带断层震级确定系统，包括：
断层位移加载速率反演模块100，用于通过反演给出我国大陆活动地块边界带主要断层段现今断层位移加载速率。
具体的，图6示出了本发明实施例提供的系统中断层位移加载速率反演模块的结构框图。
其中，在本发明提供的优选实施方式中，所述断层位移加载速率反演模块100具体包括：
断层模型构建单元1001，用于基于各断层段几何产状参数和闭锁深度参数构建用于反演计算的断层模型。
观测资料约束单元1002，用于将大地测量观测的地壳运动速率作为反演约束进行引入，并且将速度误差超过第一预设值以及距离断层距离小于第二预设值的点位数据舍弃。
反演计算单元1003，用于基于弹性回跳和块体运动理论及断层连接元模型，反演计算断层位移加载速率，主要公式为：
x＝[ATC-1A+BTD-1B]-1ATC-1y，
其中，x为断层错动分量构成的向量，y是大地测量观测的东向和北向速度分量构成的速度场向量，A为各个断层错动分量对某一测站速度分量的贡献构成的传播矩阵，B为包含了结点处位移连续性条件和拉张分量限制条件的断层系统信息的传播矩阵，C为各个测站速度分量的方差-协方差构成的矩阵，D为各结点位移连续性误差的方差矩阵。
进一步的，所述一种活动地块边界带断层震级确定系统还包括：
末次强震离逝时间确定模块200，用于根据断层段末次强震发生时间的记录以及区域内地震目录完备的起始时间获得各断层段末次强震离逝时间。
具体的，图7示出了本发明实施例提供的系统中末次强震离逝时间确定模块的结构框图。
其中，在本发明提供的优选实施方式中，所述末次强震离逝时间确定模块200具体包括：
直接记录单元2001，用于对存在末次强震发生时间记录的断层段，直接使用记录计算获得末次强震离逝时间。
期望估计单元2002，用于对缺乏末次强震发生时间记录的断层段，进行末次强震离逝时间期望值的估计。
具体的，图8示出了本发明实施例提供的系统中期望估计单元的结构框图。
其中，在本发明提供的优选实施方式中，所述期望估计单元2002具体包括：
目录完备时间子单元20021，用于依据区域研究结果确定各断层段地震目录完备的起始时间。
期望值单独求解子单元20022，用于根据公式
Figure BDA0003782402270000111
估计末次强震离逝时间的期望值，其中ξ(τ|τ≥TH)为末次强震离逝时间的期望值，TH为地震目录完备的起始时间，F为累积概率密度函数，分别采用对数正态和布朗过程时间两种模型进行计算。
期望值平均子单元20023，用于将基于对数正态和布朗时间过程两种模型的强震离逝时间期望值结果进行平均，以该平均值作为最终结果。
进一步的，所述一种活动地块边界带断层震级确定系统还包括：
累积位移量计算模块300，用于根据各断层段的断层位移加载速率和末次强震离逝时间，计算各断层段累积位移量。
在本系统中，累积位移量计算模块300对于各断层段累积位移量由公式S＝v*T计算，其中S为各断层段累积位移量，v为各断层段的断层位移加载速率，T为各断层段末次强震离逝时间。
进一步的，所述一种活动地块边界带断层震级确定系统还包括：
震级确定模块400，用于基于位错可预测的强震复发模型，依据震级与位错的经验公式确定主要断层段的发震震级。
在本系统中，震级确定模块400基于位错可预测的强震复发模型确定震级，该模型假设每次地震会将应力彻底释放至稳定值，强震离逝时间越长，断层加载速率越高，则累积位移量越大，预期发生的地震震级就越高，其计算结果是时间依赖的。依据现有技术中震级与位错的经验公式确定主要断层段的发震震级：
M＝a+b log10(AD)，
其中，M为预测震级，AD为断层段累积位移量；a和b为线性回归系数，a取6.93±0.05，b取0.82±0.1。
在一个实施例中，提出了一种计算机设备，所述计算机设备包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现以下步骤：
通过反演给出我国大陆活动地块边界带主要断层段现今断层位移加载速率；
根据断层段末次强震发生时间的记录以及区域内地震目录完备的起始时间获得各断层段末次强震离逝时间；
根据各断层段的断层位移加载速率和末次强震离逝时间，计算各断层段累积位移量；
基于位错可预测的强震复发模型，依据震级与位错的经验公式确定主要断层段的发震震级。
应该理解的是，虽然本发明各实施例的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，各实施例中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一非易失性计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (12)

1.一种活动地块边界带断层震级确定方法，其特征在于，所述方法包括：
通过反演给出我国大陆活动地块边界带主要断层段现今断层位移加载速率；
根据断层段末次强震发生时间的记录以及区域内地震目录完备的起始时间获得各断层段末次强震离逝时间；
根据各断层段的断层位移加载速率和末次强震离逝时间，计算各断层段累积位移量；
基于位错可预测的强震复发模型，依据震级与位错的经验公式确定主要断层段的发震震级。
2.根据权利要求1所述的活动地块边界带断层震级确定方法，其特征在于，所述通过反演给出我国大陆活动地块边界带主要断层段现今断层位移加载速率的步骤中：
基于各断层段几何产状参数和闭锁深度参数构建用于反演计算的断层模型；
将大地测量观测的地壳运动速率作为反演约束进行引入，并且将速度误差超过第一预设值以及距离断层距离小于第二预设值的点位数据舍弃；
基于弹性回跳和块体运动理论及断层连接元模型，反演计算断层位移加载速率，主要公式为：
x＝[ATC-1A+BTD-1B]-1ATC-1y，
其中，x为断层错动分量构成的向量，y是大地测量观测的东向和北向速度分量构成的速度场向量，A为各个断层错动分量对某一测站速度分量的贡献构成的传播矩阵，B为包含了结点处位移连续性条件和拉张分量限制条件的断层系统信息的传播矩阵，C为各个测站速度分量的方差-协方差构成的矩阵，D为各结点位移连续性误差的方差矩阵。
3.根据权利要求1所述的活动地块边界带断层震级确定方法，其特征在于，所述根据断层段末次强震发生时间的记录以及区域内地震目录完备的起始时间获得各断层段末次强震离逝时间的步骤中，对于存在末次强震发生时间记录的断层段，直接使用记录计算获得末次强震离逝时间。
4.根据权利要求1所述的活动地块边界带断层震级确定方法，其特征在于，所述根据断层段末次强震发生时间的记录以及区域内地震目录完备的起始时间获得各断层段末次强震离逝时间的步骤中，对于缺乏末次强震发生时间记录的断层段，进行末次强震离逝时间期望值的估计。
5.根据权利要求4所述的活动地块边界带断层震级确定方法，其特征在于，所述进行末次强震离逝时间期望值的估计的过程，具体包括：
依据区域研究结果确定各断层段地震目录完备的起始时间；
根据公式
Figure FDA0003782402260000021
估计末次强震离逝时间的期望值，其中ξ(τ|τ≥TH)为末次强震离逝时间的期望值，TH为地震目录完备的起始时间，F为累积概率密度函数，分别采用对数正态和布朗过程时间两种模型进行计算；
以计算结果平均值作为最终结果。
6.根据权利要求1所述的活动地块边界带断层震级确定方法，其特征在于，各断层段累积位移量由公式S＝v*T计算，其中S为各断层段累积位移量，v为各断层段的断层位移加载速率，T为各断层段末次强震离逝时间。
7.根据权利要求1所述的活动地块边界带断层震级确定方法，其特征在于，通过基于位错可预测的强震复发模型确定震级，在该基于位错可预测的强震复发模型种，假设每次地震会将应力彻底释放至稳定值，强震离逝时间越长，断层加载速率越高，则累积位移量越大，预期发生的地震震级就越高，其计算结果是时间依赖的。
8.根据权利要求1所述的活动地块边界带断层震级确定方法，其特征在于，所述震级与位错的经验公式为：
M＝a+blog10(AD)，
其中，M为预测震级，AD为断层段累积位移量；a和b为线性回归系数，a取6.93±0.05，b取0.82±0.1。
9.一种活动地块边界带断层震级确定系统，其特征在于，所述系统包括：
断层位移加载速率反演模块，用于通过反演给出我国大陆活动地块边界带主要断层段现今断层位移加载速率；
末次强震离逝时间确定模块，用于根据断层段末次强震发生时间的记录以及区域内地震目录完备的起始时间获得各断层段末次强震离逝时间；
累积位移量计算模块，用于根据各断层段的断层位移加载速率和末次强震离逝时间，计算各断层段累积位移量；
震级确定模块，用于基于位错可预测的强震复发模型，依据震级与位错的经验公式确定主要断层段的发震震级。
10.根据权利要求9所述的活动地块边界带断层震级确定系统，其特征在于，所述断层位移加载速率反演模块包括：
断层模型构建单元，用于基于各断层段几何产状参数和闭锁深度参数构建用于反演计算的断层模型。
观测资料约束单元，用于将大地测量观测的地壳运动速率作为反演约束进行引入，并且将速度误差超过第一预设值以及距离断层距离小于第二预设值的点位数据舍弃。
反演计算单元，用于基于弹性回跳和块体运动理论及断层连接元模型，反演计算断层位移加载速率，主要公式为：
x＝[ATC-1A+BTD-1B]-1ATC-1y，
其中，x为断层错动分量构成的向量，y是大地测量观测的东向和北向速度分量构成的速度场向量，A为各个断层错动分量对某一测站速度分量的贡献构成的传播矩阵，B为包含了结点处位移连续性条件和拉张分量限制条件的断层系统信息的传播矩阵，C为各个测站速度分量的方差-协方差构成的矩阵，D为各结点位移连续性误差的方差矩阵。
11.根据权利要求9所述的活动地块边界带断层震级确定系统，其特征在于，所述末次强震离逝时间确定模块包括：
直接记录单元，用于对存在末次强震发生时间记录的断层段，直接使用记录计算获得末次强震离逝时间。
期望估计单元，用于对缺乏末次强震发生时间记录的断层段，进行末次强震离逝时间期望值的估计。
12.根据权利要求11所述的活动地块边界带断层震级确定系统，其特征在于，所述期望估计单元包括：
目录完备时间子单元，用于依据区域研究结果确定各断层段地震目录完备的起始时间。
期望值单独求解子单元，用于根据公式
Figure FDA0003782402260000041
估计末次强震离逝时间的期望值，其中ξ(τ|τ≥TH)为末次强震离逝时间的期望值，TH为地震目录完备的起始时间，F为累积概率密度函数，分别采用对数正态和布朗过程时间两种模型进行计算。
期望值平均子单元，用于将基于对数正态和布朗时间过程两种模型的强震离逝时间期望值结果进行平均，以该平均值作为最终结果。
CN202210932909.5A 2022-08-04 2022-08-04 一种活动地块边界带断层震级确定方法及系统 Pending CN115343751A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210932909.5A CN115343751A (zh) 2022-08-04 2022-08-04 一种活动地块边界带断层震级确定方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210932909.5A CN115343751A (zh) 2022-08-04 2022-08-04 一种活动地块边界带断层震级确定方法及系统

Publications (1)

Publication Number Publication Date
CN115343751A true CN115343751A (zh) 2022-11-15

Family

ID=83950593

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210932909.5A Pending CN115343751A (zh) 2022-08-04 2022-08-04 一种活动地块边界带断层震级确定方法及系统

Country Status (1)

Country Link
CN (1) CN115343751A (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014092454A (ja) * 2012-11-02 2014-05-19 Takeda Engineering Consultant:Kk 大地震と巨大地震の予知方法、予知装置、予知プログラム及び記録媒体
JP2018036254A (ja) * 2016-09-01 2018-03-08 公益財団法人鉄道総合技術研究所 早期地震警報のための地震規模決定方法
CN111382908A (zh) * 2020-03-11 2020-07-07 中国地震局地球物理研究所 一种考虑大震时间相关性的地震随机事件集模拟方法
US20210302603A1 (en) * 2020-03-30 2021-09-30 Qingdao university of technology Performance-level seismic motion hazard analysis method based on three-layer dataset neural network
CN114675325A (zh) * 2022-03-13 2022-06-28 中国地震台网中心 基于Mindlin解的地震地表破裂永久位移的估计方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014092454A (ja) * 2012-11-02 2014-05-19 Takeda Engineering Consultant:Kk 大地震と巨大地震の予知方法、予知装置、予知プログラム及び記録媒体
JP2018036254A (ja) * 2016-09-01 2018-03-08 公益財団法人鉄道総合技術研究所 早期地震警報のための地震規模決定方法
CN111382908A (zh) * 2020-03-11 2020-07-07 中国地震局地球物理研究所 一种考虑大震时间相关性的地震随机事件集模拟方法
US20210302603A1 (en) * 2020-03-30 2021-09-30 Qingdao university of technology Performance-level seismic motion hazard analysis method based on three-layer dataset neural network
CN114675325A (zh) * 2022-03-13 2022-06-28 中国地震台网中心 基于Mindlin解的地震地表破裂永久位移的估计方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
朱亚戈 等: "基于GPS资料约束的2021年玛多地震发震断层的滑动速率", 中国科学, vol. 51, no. 10, pages 1790 - 1791 *
李献瑞: "青藏高原东缘三维构造应力场数值模拟及地震危险性分析", 中国博士学位论文全文数据库 基础科学辑, no. 03, pages 99 *
王阎昭 等: "基于 GPS 资料约束反演川滇地区主要断裂现今活动速率", 中国科学, vol. 38, no. 5, pages 584 - 586 *
贺鹏超等: "基于2001年Mw7. 8可可西里地震震后形变模拟研究藏北地区岩石圈流变学结构", 地球物理学报, vol. 61, no. 2, pages 537 *

Similar Documents

Publication Publication Date Title
Ghosh et al. Consideration of time-evolving capacity distributions and improved degradation models for seismic fragility assessment of aging highway bridges
Ryu et al. Developing fragilities for mainshock-damaged structures through incremental dynamic analysis
Hatzigeorgiou et al. Evaluation of maximum seismic displacements of SDOF systems from their residual deformation
Val et al. Reliability evaluation in nonlinear analysis of reinforced concrete structures
Pnevmatikos et al. Fragility curves for mixed concrete/steel frames subjected to seismic excitation
Bojórquez et al. Structural reliability of reinforced concrete buildings under earthquakes and corrosion effects
Jiang et al. Seismic vulnerability analysis for shallow-buried underground frame structure considering 18 existing subway stations
Walsh et al. Geometric characterisation and out-of-plane seismic stability of low-rise unreinforced brick masonry buildings in Auckland, New Zealand
Kia et al. Rapid seismic vulnerability assessment by new regression-based demand and collapse models for steel moment frames
Han et al. Seismic loss estimation with consideration of aftershock hazard and post-quake decisions
Hu et al. Integrated multi‐type sensor placement and response reconstruction method for high‐rise buildings under unknown seismic loading
CN115343751A (zh) 一种活动地块边界带断层震级确定方法及系统
Kia et al. Closed-form fragility analysis of the steel moment resisting frames
Bioud et al. ESTIMATING THE FUNDAMENTAL PERIOD OF INFILLED RC FRAME STRUCTURES VIA DEEP LEARNING.
Wang et al. Correlations between structural damage and ground motion parameters during the Ms8. 0 Wenchuan earthquake
Liu et al. Implications of bidirectional interaction on nonlinear seismic response of steel piers
Argyroudis et al. Seismic fragility curves of shallow tunnels considering SSI and aging effects
Ródenas et al. Advances in seismic vulnerability assessment of reinforced concrete buildings applied to the experience of Lorca (Spain) 2011 earthquake
JP7323468B2 (ja) 対象物測位システム、及び対象物測位方法
Shekhar et al. Influence of corrosion on failure modes and lifetime seismic vulnerability assessment of low‐ductility RC frames
Billah et al. Mainshock-aftershock damage assessment of concrete bridge reinforced with shape memory alloy rebar
Shafei et al. Assessment of extent of capacity loss in deteriorated highway bridges
Thinley et al. Seismic performance of reinforced concrete buildings in Thimphu, Bhutan
Cabrera et al. Empirical fragility curves for masonry houses using data from two earthquakes in Chile
Al Shamsi Seismic risk assessment of buildings in Dubai, United Arab Emirates

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination