CN106772600B - 一种双对双差地震定位方法及装置 - Google Patents
一种双对双差地震定位方法及装置 Download PDFInfo
- Publication number
- CN106772600B CN106772600B CN201611192581.9A CN201611192581A CN106772600B CN 106772600 B CN106772600 B CN 106772600B CN 201611192581 A CN201611192581 A CN 201611192581A CN 106772600 B CN106772600 B CN 106772600B
- Authority
- CN
- China
- Prior art keywords
- pair
- station
- seismic
- double
- event
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 239000011159 matrix material Substances 0.000 claims description 34
- 238000004364 calculation method Methods 0.000 claims description 28
- 230000009977 dual effect Effects 0.000 claims description 19
- 239000000126 substance Substances 0.000 claims description 10
- 238000010276 construction Methods 0.000 claims description 5
- 238000010219 correlation analysis Methods 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 4
- 230000010355 oscillation Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000009189 diving Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
Images
Classifications
-
- 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. for interpretation or for event detection
- G01V1/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/12—Signal generation
- G01V2210/129—Source location
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供的双对双差地震定位方法及装置,提取观测的地震事件对和台站对的到时差数据,进而提取观测的双对到时差数据;计算各个地震事件对双差、各个台站对双差,以及地震事件对到台站对的双对双差;双对双差可以消除发震时刻的影响,无需引入台站校正项,同时消除了源区外和源区内的速度模型的影响,得到高精度的相对位置,而台站对双差可以得到高精度的绝对位置。因此,根据台站对双差与地震位置的线性关系,以及双对双差与地震位置的线性关系,构建反演***,并采取分级加权的方式来反演各个地震事件的位置,可以得到高精度的地震绝对位置以及地震间的相对位置。
Description
技术领域
本发明涉及地震定位技术领域,更具体地,涉及一种双对双差地震定位方法及装置。
背景技术
地震信号的精确定位对于监测地震、研究活动断层、俯冲板块、火山,以及其他与地震活动相关的结构有重要意义。
目前,最流行的地震定位算法是Waldhauser&Ellsworth(2000)提出的事件对双差定位算法,该算法利用两个相邻事件到同一个台站的地震目录或者互相关走时差信息,消除源区外的速度模型影响,可以明显提高地震间的相对位置。但该方法仍受到源区内的速度模型影响,并且其绝对位置较差。
Zhang et al.(2010)提出了台站对双差定位算法。该方法利用同一个事件到两个台站的走时差数据,可以消除发震时刻以及源区内的速度模型的影响,可以得到比较好的绝对位置。但该方法受源区外的速度模型影响很大,地震间的相对位置精度较差。
因此,当前的地震定位方法得到的地震位置精度有限,需要进一步的提高。
发明内容
有鉴于此,本发明提供了一种双对双差地震定位方法及装置,利用地震事件对到台站对的双对双差,消除发震时刻的影响,也无需引入台站校正项,同时消除了源区外和源区内的速度模型的影响,得到高精度的地震间相对位置;并通过联合台站对双差以及双对双差来同时得到高精度的地震绝对位置和地震间相对位置。
具体技术方案如下:
一种双对双差地震定位方法,所述方法包括:
基于地震的初至波绝对到时数据或者地震波形的互相关数据,按预设规则,提取观测的地震事件对(i,j)到同一台站的到时差数据以及同一地震到台站对(k,l)的到时差数据,进而提取观测的双对到时差数据;所述双对为所述地震事件对(i,j)和所述台站对(k,l)的组合,其中,i和j分别表示地震事件对中的两个地震事件,地震事件i和地震事件j的位置相邻,l和k分别表示台站对中的两个台站;
根据台站对双差与地震位置相对于地震初始位置的变化量的线性关系,以及双对双差与地震位置相对于地震初始位置的变化量的线性关系,构建反演***;
采取分级加权的方式来反演得到各个地震事件的地震位置。
根据所述地震事件对(i,j)中地震事件i和地震事件j到台站k的到时和的观测值与预测值的差,得到所述地震事件i和所述地震事件j到台站k的到时残差和并根据所述地震事件i和所述地震事件j到所述台站k的到时残差和得到所述地震事件对(i,j)到所述台站k的事件对双差所述地震事件对(i,j)为多个地震事件对中的任意一个地震事件对,
根据所述地震事件i和所述地震事件j到台站l的到时和的观测值与预测值的差,得到所述地震事件i和所述地震事件j到台站l的到时残差和并根据所述地震事件i和所述地震事件j到所述台站l的到时残差和得到所述地震事件对(i,j)到所述台站l的事件对双差
根据所述地震事件i到所述台站k和所述台站l的到时和的观测值与预测值的差,得到所述地震事件i到所述台站k和所述台站l的到时残差和并根据所述地震事件i到所述台站k和所述台站l的到时残差和得到地震i到所述台站对(k,l)的台站对双差所述台站对(k,l)为多个台站对中的任意一个台站对;
根据所述地震事件j到所述台站k和所述台站l的到时和的观测值与预测值的差,得到所述地震事件j到所述台站k和所述台站l的到时残差和并根据所述地震事件j到所述台站k和所述台站l的到时残差和得到地震j到所述台站对(k,l)的台站对双差
计算地震事件对(i,j)到台站k的事件对双差与地震事件对(i,j)到台站l的事件对双差的差值,或计算地震事件i到台站对(k,l)的台站对双差与地震事件j到台站对(k,l)的台站对双差的差值,得到所述地震事件对(i,j)到所述台站对(k,l)的双对双差所述地震事件对(i,j)与所述台站对(k,l)的双对组合为多个双对中的任意一个双对;
优选的,所述根据台站对双差与地震位置相对于地震初始位置的变化量的线性关系,以及双对双差与地震位置相对于地震初始位置的变化量的线性关系,构建反演***,包括:
将上述两个线性关系公式以矩阵的形式表示,构成反演***:
其中,G为台站对到时差与双对到时差数据对地震位置的偏导数矩阵,所述偏导数矩阵根据预先设定的各个所述地震事件的初始位置计算得到,所述偏导数矩阵的维度是M*3N,M为所述台站对到时差数据数量与双对到时差数据数量的总和,N为所述地震事件数量,m是N个地震事件的位置在x,y,z三个方向上的变化量向量[dx,dy,dz],d为各个所述台站对双差数据与各个所述双对双差数据,W为权重矩阵,λ为阻力因子,I为单位阵;
为所述台站对双差设置第一权重,为所述双对双差设置第二权重。
优选的,所述采取分级加权的方式来反演得到各个地震事件的地震位置,包括:
根据LSQR算法求解所述反演***,得到各个所述地震事件的地震位置,并将各个所述地震事件的地震位置分别作为各个所述地震事件的新的初始位置,开始迭代求解各个所述地震事件的地震位置,具体为:在初期迭代中,使所述台站对双差的第一权重大于所述双对双差的第二权重,直到迭代次数达到第一预设次数;在后期迭代中,使所述双对双差的第二权重大于所述台站对双差的第一权重,直到迭代次数达到第二预设次数,迭代结束,得到最终各个所述地震事件的地震位置。
一种双对双差地震定位装置,所述装置包括:
提取单元,用于基于地震的初至波绝对到时数据或者地震波形的互相关数据,按预设规则,提取观测的地震事件对(i,j)到同一台站的到时差数据以及同一地震到台站对(k,l)的到时差数据,进而提取观测的双对到时差数据;所述双对为所述地震事件对(i,j)和所述台站对(k,l)的组合,其中,i和j分别表示地震事件对中的两个地震事件,地震事件i和地震事件j的位置相邻,l和k分别表示台站对中的两个台站;
构建单元,用于根据台站对双差与地震位置相对于地震初始位置的变化量的线性关系,以及双对双差与地震位置相对于地震初始位置的变化量的线性关系,构建反演***;
第四计算单元,用于采取分级加权的方式来反演得到各个地震事件的地震位置。
优选的,所述第一计算单元包括:
第一计算子单元,用于根据所述地震事件对(i,j)中地震事件i和地震事件j到台站k的到时和的观测值与预测值的差,得到所述地震事件i和所述地震事件j到台站k的到时残差和并根据所述地震事件i和所述地震事件j到所述台站k的到时残差和得到所述地震事件对(i,j)到所述台站k的事件对双差所述地震事件对(i,j)为多个地震事件对中的任意一个地震事件对,
第二计算子单元,用于根据所述地震事件i和所述地震事件j到台站l的到时和的观测值与预测值的差,得到所述地震事件i和所述地震事件j到台站l的到时残差和并根据所述地震事件i和所述地震事件j到所述台站l的到时残差和得到所述地震事件对(i,j)到所述台站l的事件对双差
优选的,所述第二计算单元包括:
第三计算子单元,用于根据所述地震事件i到所述台站k和所述台站l的到时和的观测值与预测值的差,得到所述地震事件i到所述台站k和所述台站l的到时残差和并根据所述地震事件i到所述台站k和所述台站l的到时残差和得到地震i到所述台站对(k,l)的台站对双差所述台站对(k,l)为多个台站对中的任意一个台站对;
第四计算子单元,用于根据所述地震事件j到所述台站k和所述台站l的到时和的观测值与预测值的差,得到所述地震事件j到所述台站k和所述台站l的到时残差和并根据所述地震事件j到所述台站k和所述台站l的到时残差和得到地震j到所述台站对(k,l)的台站对双差
优选的,所述第三计算单元具体用于:计算地震事件对(i,j)到台站k的事件对双差与地震事件对(i,j)到台站l的事件对双差的差值,或计算地震事件i到台站对(k,l)的台站对双差与地震事件j到台站对(k,l)的台站对双差的差值,得到所述地震事件对(i,j)到所述台站对(k,l)的双对双差所述地震事件对(i,j)与所述台站对(k,l)的双对组合为多个双对中的任意一个双对;
优选的,所述构建单元包括:
其中,G为台站对到时差与双对到时差数据对地震位置的偏导数矩阵,所述偏导数矩阵根据预先设定的各个所述地震事件的初始位置计算得到,所述偏导数矩阵的维度是M*3N,M为所述台站对到时差数据数量与双对到时差数据数量的总和,N为所述地震事件数量,m是N个地震事件的位置在x,y,z三个方向上的变化量向量[dx,dy,dz],d为各个所述台站对双差数据与各个所述双对双差数据,W为权重矩阵,λ为阻力因子,I为单位阵;
设置子单元,用于为所述台站对双差设置第一权重,为所述双对双差设置第二权重。
优选的,所述第四计算单元具体用于:根据LSQR算法求解所述反演***,得到各个所述地震事件的地震位置,并将各个所述地震事件的地震位置分别作为各个所述地震事件的新的初始位置,开始迭代求解各个所述地震事件的地震位置,具体为:在初期迭代中,使所述台站对双差的第一权重大于所述双对双差的第二权重,直到迭代次数达到第一预设次数;在后期迭代中,使所述双对双差的第二权重大于所述台站对双差的第一权重,直到迭代次数达到第二预设次数,迭代结束,得到最终各个所述地震事件的地震位置。
相对于现有技术,本发明的有益效果如下:
本发明提供的双对双差地震定位方法及装置,基于地震的初至波绝对到时数据或者地震波形的互相关数据,按预设规则,提取观测的地震事件对和台站对的到时差数据,进而提取观测的双对到时差数据;计算各个地震事件对双差、各个台站对双差,以及地震事件对到台站对的双对双差;双对双差可以消除发震时刻的影响,无需引入台站校正项,同时消除了源区外和源区内的速度模型的影响,得到高精度的相对位置,而台站对双差可以得到高精度的绝对位置。因此,根据台站对双差与地震位置相对于地震初始位置的变化量的线性关系,以及双对双差与地震位置相对于地震初始位置的变化量的线性关系,构建反演***,并采取分级加权的方式来反演各个地震事件的位置,可以得到高精度的地震绝对位置以及地震间的相对位置。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1为本发明实施例公开的一种双对双差定位方法流程图;
图2为本发明实施例公开的一种双对双差定位装置结构示意图;
图3为本发明实施例公开的一种双对双差定位装置结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参阅图1,本实施例公开的一种双对双差地震定位方法,需要说明的是,本实施例中的待处理地震数据为已发生地震的地震记录数据,具体包括以下步骤:
S101:基于地震的初至波绝对到时数据或者地震波形的互相关数据,按预设规则,提取观测的地震事件对(i,j)到同一台站的到时差数据以及同一地震到台站对(k,l)的到时差数据,进而提取观测的双对到时差数据;
具体的,所述双对为所述地震事件对(i,j)和所述台站对(k,l)的组合,其中,i和j分别表示地震事件对中的两个地震事件,地震事件i和地震事件j的位置相邻,l和k分别表示台站对中的两个台站;
需要说明的是,现有技术中的事件对双差定位算法和台站对双差定位算法提供了构建地震事件对和构建台站对的规则,根据实际需要,用户也可以自定义构建地震事件对和台站对的规则,由于构建地震事件对和台站对的规则并不是本实施例研究的重点,在此不再赘述。双对基于地震事件对和台站对进行组合,用户也可以根据实际需要自定义构建双对的规则。
根据所述地震事件对(i,j)中地震事件i和地震事件j到台站k的到时和的观测值与预测值的差,得到所述地震事件i和所述地震事件j到台站k的到时残差和并根据所述地震事件i和所述地震事件j到所述台站k的到时残差和得到所述地震事件对(i,j)到所述台站k的事件对双差所述地震事件对(i,j)为多个地震事件对中的任意一个地震事件对,
具体的,地震事件i的体波到台站k的到时可以表示为其射线路径的积分:其中,地震事件位置,发震时刻τi,射线路径,慢度场u是未知数。到时与地震事件位置参数之间是非线性的,因此我们用一个截断的泰勒展开式来线性表示公式观测到时与预测到时的残差可以与地震位置以及发震时刻参数建立线性联系,即,
根据所述地震事件i和所述地震事件j到台站l的到时和的观测值与预测值的差,得到所述地震事件i和所述地震事件j到台站l的到时残差和并根据所述地震事件i和所述地震事件j到所述台站l的到时残差和得到所述地震事件对(i,j)到所述台站l的事件对双差
根据所述地震事件i到所述台站k和所述台站l的到时和的观测值与预测值的差,得到所述地震事件i到所述台站k和所述台站l的到时残差和并根据所述地震事件i到所述台站k和所述台站l的到时残差和得到地震i到所述台站对(k,l)的台站对双差所述台站对(k,l)为多个台站对中的任意一个台站对;
根据所述地震事件j到所述台站k和所述台站l的到时和的观测值与预测值的差,得到所述地震事件j到所述台站k和所述台站l的到时残差和并根据所述地震事件j到所述台站k和所述台站l的到时残差和得到地震j到所述台站对(k,l)的台站对双差
可见,地震事件的发震时刻项被消除了。
具体的,计算地震事件对(i,j)到台站k的事件对双差与地震事件对(i,j)到台站l的事件对双差的差值,或计算地震事件i到台站对(k,l)的台站对双差与地震事件j到台站对(k,l)的台站对双差的差值,得到所述地震事件对(i,j)到所述台站对(k,l)的双对双差所述地震事件对(i,j)与所述台站对(k,l)的双对组合为多个双对中的任意一个双对;
可见,利用地震事件对到台站对的双对双差,消除发震时刻的影响,也无需引入台站校正项,同时消除了源区外和源区内的速度模型的影响。
S105:根据台站对双差与地震位置相对于地震初始位置的变化量的线性关系,以及双对双差与地震位置相对于地震初始位置的变化量的线性关系,构建反演***;
具体的,步骤S105的具体执行过程如下:
将上述两个线性关系公式以矩阵的形式表示,构成反演***:
其中,G为台站对到时差与双对到时差数据对地震位置的偏导数矩阵,所述偏导数矩阵根据预先设定的各个所述地震事件的初始位置计算得到,所述偏导数矩阵的维度是M*3N,M为所述台站对到时差数据数量与双对到时差数据数量的总和,N为所述地震事件数量,m是N个地震事件的位置在x,y,z三个方向上的变化量向量[dx,dy,dz],d为各个所述台站对双差数据与各个所述双对双差数据,W为权重矩阵,λ为阻力因子,I为单位阵;
为所述台站对双差设置第一权重,为所述双对双差设置第二权重。
S106:采取分级加权的方式来反演得到各个地震事件的地震位置。
具体的,根据LSQR算法求解所述反演***,得到各个所述地震事件的地震位置,并将各个所述地震事件的地震位置分别作为各个所述地震事件的新的初始位置,开始迭代求解各个所述地震事件的地震位置,具体为:在初期迭代中,使所述台站对双差的第一权重大于所述双对双差的第二权重,直到迭代次数达到第一预设次数;在后期迭代中,使所述双对双差的第二权重大于所述台站对双差的第一权重,直到迭代次数达到第二预设次数,迭代结束,得到最终各个所述地震事件的地震位置。
需要说明的是,初期迭代的第一预设次数以及后期迭代的第二预设次数可以根据实际情况进行调整,同时,第一权重和第二权重的值也可以根据实际需要进行设定。
还需要说明的是,实际应用中,用户可以选择获取事件对双差或台站对双差计算双对双差。
本实施例提供的双对双差地震定位方法,基于地震的初至波绝对到时数据或者地震波形的互相关数据,按预设规则,提取观测的地震事件对和台站对的到时差数据,进而提取观测的双对到时差数据;计算各个地震事件对双差、各个台站对双差,以及地震事件对到台站对的双对双差;双对双差可以消除发震时刻的影响,无需引入台站校正项,同时消除了源区外和源区内的速度模型的影响,得到高精度的相对位置,而台站对双差可以得到高精度的绝对位置。因此,根据台站对双差与地震位置相对于地震初始位置的变化量的线性关系,以及双对双差与地震位置相对于地震初始位置的变化量的线性关系,构建反演***,并采取分级加权的方式来反演各个地震事件的位置,可以得到高精度的地震绝对位置以及地震间的相对位置。
基于上述实施例公开的一种双对双差地震定位方法,本实施例对应公开了一种双对双差地震定位装置,请参阅图2,所述装置具体包括:
提取单元101,用于基于地震的初至波绝对到时数据或者地震波形的互相关数据,按预设规则,提取观测的地震事件对(i,j)到同一台站的到时差数据以及同一地震到台站对(k,l)的到时差数据,进而提取观测的双对到时差数据;所述双对为所述地震事件对(i,j)和所述台站对(k,l)的组合,其中,i和j分别表示地震事件对中的两个地震事件,地震事件i和地震事件j的位置相邻,l和k分别表示台站对中的两个台站;
所述第三计算单元104具体用于:计算地震事件对(i,j)到台站k的事件对双差与地震事件对(i,j)到台站l的事件对双差的差值,或计算地震事件i到台站对(k,l)的台站对双差与地震事件j到台站对(k,l)的台站对双差的差值,得到所述地震事件对(i,j)到所述台站对(k,l)的双对双差所述地震事件对(i,j)与所述台站对(k,l)的双对组合为多个双对中的任意一个双对;
构建单元105,用于根据台站对双差与地震位置相对于地震初始位置的变化量的线性关系,以及双对双差与地震位置相对于地震初始位置的变化量的线性关系,构建反演***;
第四计算单元106,用于采取分级加权的方式来反演得到各个地震事件的地震位置。
所述第四计算单元106具体用于:根据LSQR算法求解所述反演***,得到各个所述地震事件的地震位置,并将各个所述地震事件的地震位置分别作为各个所述地震事件的新的初始位置,开始迭代求解各个所述地震事件的地震位置,具体为:在初期迭代中,使所述台站对双差的第一权重大于所述双对双差的第二权重,直到迭代次数达到第一预设次数;在后期迭代中,使所述双对双差的第二权重大于所述台站对双差的第一权重,直到迭代次数达到第二预设次数,迭代结束,得到最终各个所述地震事件的地震位置。
需要说明的是,实际计算时用户可以选择第一计算单元102或第二计算单元103计算得出计算单元104所需的数据。
请参阅图3,所述双对双差地震定位装置中的第一计算单元102包括:
第一计算子单元107,用于根据所述地震事件对(i,j)中地震事件i和地震事件j到台站k的到时和的观测值与预测值的差,得到所述地震事件i和所述地震事件j到台站k的到时残差和并根据所述地震事件i和所述地震事件j到所述台站k的到时残差和得到所述地震事件对(i,j)到所述台站k的事件对双差所述地震事件对(i,j)为多个地震事件对中的任意一个地震事件对,
第二计算子单元108,用于根据所述地震事件i和所述地震事件j到台站l的到时和的观测值与预测值的差,得到所述地震事件i和所述地震事件j到台站l的到时残差和并根据所述地震事件i和所述地震事件j到所述台站l的到时残差和得到所述地震事件对(i,j)到所述台站l的事件对双差
第二计算单元103包括:
第三计算子单元109,用于根据所述地震事件i到所述台站k和所述台站l的到时和的观测值与预测值的差,得到所述地震事件i到所述台站k和所述台站l的到时残差和并根据所述地震事件i到所述台站k和所述台站l的到时残差和得到地震i到所述台站对(k,l)的台站对双差所述台站对(k,l)为多个台站对中的任意一个台站对;
第四计算子单元110,用于根据所述地震事件j到所述台站k和所述台站l的到时和的观测值与预测值的差,得到所述地震事件j到所述台站k和所述台站l的到时残差和并根据所述地震事件j到所述台站k和所述台站l的到时残差和得到地震j到所述台站对(k,l)的台站对双差
所述构建单元105包括:
其中,G为台站对到时差与双对到时差数据对地震位置的偏导数矩阵,所述偏导数矩阵根据预先设定的各个所述地震事件的初始位置计算得到,所述偏导数矩阵的维度是M*3N,M为所述台站对到时差数据数量与双对到时差数据数量的总和,N为所述地震事件数量,m是N个地震事件的位置在x,y,z三个方向上的变化量向量[dx,dy,dz],d为各个所述台站对双差数据与各个所述双对双差数据,W为权重矩阵,λ为阻力因子,I为单位阵;
设置子单元112,用于为所述台站对双差设置第一权重,为所述双对双差设置第二权重。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (8)
1.一种双对双差地震定位方法,其特征在于,所述方法包括:
基于地震的初至波绝对到时数据或者地震波形的互相关数据,按预设规则,提取观测的地震事件对(i,j)到同一台站的到时差数据以及同一地震到台站对(k,l)的到时差数据,进而提取观测的双对到时差数据;按所述预设规则,所述双对基于地震事件对和台站对进行组合;所述双对为所述地震事件对(i,j)和所述台站对(k,l)的组合,其中,i和j分别表示地震事件对中的两个地震事件,地震事件i和地震事件j的位置相邻,l和k分别表示台站对中的两个台站;
根据台站对双差与地震位置相对于地震初始位置的变化量的线性关系,以及双对双差与地震位置相对于地震初始位置的变化量的线性关系,构建反演***:
其中,G为台站对到时差与双对到时差数据对地震位置的偏导数矩阵,所述偏导数矩阵根据预先设定的各个地震事件的初始位置计算得到,所述偏导数矩阵的维度是M*3N,M为台站对到时差数据数量与双对到时差数据数量的总和,N为所述地震事件数量,m是N个地震事件的位置在x,y,z三个方向上的变化量向量[dx,dy,dz],d为各个台站对双差数据与各个双对双差数据,W为权重矩阵,λ为阻力因子,I为单位阵;
为所述台站对双差设置第一权重,为所述双对双差设置第二权重;
根据LSQR算法求解所述反演***,得到各个所述地震事件的地震位置,并将各个所述地震事件的地震位置分别作为各个所述地震事件的新的初始位置,开始迭代求解各个所述地震事件的地震位置,具体为:在初期迭代中,使所述台站对双差的第一权重大于所述双对双差的第二权重,直到迭代次数达到第一预设次数;在后期迭代中,使所述双对双差的第二权重大于所述台站对双差的第一权重,直到迭代次数达到第二预设次数,迭代结束,得到最终各个所述地震事件的地震位置,其中所述地震位置包括:地震绝对位置以及地震间的相对位置。
根据所述地震事件对(i,j)中地震事件i和地震事件j到台站k的到时和的观测值与预测值的差,得到所述地震事件i和所述地震事件j到台站k的到时残差和并根据所述地震事件i和所述地震事件j到所述台站k的到时残差和得到所述地震事件对(i,j)到所述台站k的事件对双差所述地震事件对(i,j)为多个地震事件对中的任意一个地震事件对,
根据所述地震事件i和所述地震事件j到台站l的到时Tl i和Tl j的观测值与预测值的差,得到所述地震事件i和所述地震事件j到台站l的到时残差rl i和rl j,并根据所述地震事件i和所述地震事件j到所述台站l的到时残差rl i和rl j,得到所述地震事件对(i,j)到所述台站l的事件对双差drl ij,
根据所述地震事件i到所述台站k和所述台站l的到时和Tl i的观测值与预测值的差,得到所述地震事件i到所述台站k和所述台站l的到时残差和rl i,并根据所述地震事件i到所述台站k和所述台站l的到时残差和rl i,得到地震i到所述台站对(k,l)的台站对双差所述台站对(k,l)为多个台站对中的任意一个台站对;
根据所述地震事件j到所述台站k和所述台站l的到时和Tl j的观测值与预测值的差,得到所述地震事件j到所述台站k和所述台站l的到时残差和rl j并根据所述地震事件j到所述台站k和所述台站l的到时残差和rl j,得到地震j到所述台站对(k,l)的台站对双差
5.一种双对双差地震定位装置,其特征在于,所述装置包括:
提取单元,用于基于地震的初至波绝对到时数据或者地震波形的互相关数据,按预设规则,提取观测的地震事件对(i,j)到同一台站的到时差数据以及同一地震到台站对(k,l)的到时差数据,进而提取观测的双对到时差数据;按所述预设规则,所述双对基于地震事件对和台站对进行组合;所述双对为所述地震事件对(i,j)和所述台站对(k,l)的组合,其中,i和j分别表示地震事件对中的两个地震事件,地震事件i和地震事件j的位置相邻,l和k分别表示台站对中的两个台站;
第一计算单元,用于计算所述地震事件对(i,j)到台站k的事件对双差和所述地震事件对(i,j)到台站l的事件对双差drl ij;所述事件对双差为所述地震事件对(i,j)的到时差数据的观测值与预测值的差;
构建单元,用于根据台站对双差与地震位置相对于地震初始位置的变化量的线性关系,以及双对双差与地震位置相对于地震初始位置的变化量的线性关系,构建反演***:
其中,G为台站对到时差与双对到时差数据对地震位置的偏导数矩阵,所述偏导数矩阵根据预先设定的各个地震事件的初始位置计算得到,所述偏导数矩阵的维度是M*3N,M为台站对到时差数据数量与双对到时差数据数量的总和,N为所述地震事件数量,m是N个地震事件的位置在x,y,z三个方向上的变化量向量[dx,dy,dz],d为各个台站对双差数据与各个双对双差数据,W为权重矩阵,λ为阻力因子,I为单位阵;
第四计算单元,用于根据LSQR算法求解所述反演***,得到各个所述地震事件的地震位置,并将各个所述地震事件的地震位置分别作为各个所述地震事件的新的初始位置,开始迭代求解各个所述地震事件的地震位置,具体为:在初期迭代中,使所述台站对双差的第一权重大于所述双对双差的第二权重,直到迭代次数达到第一预设次数;在后期迭代中,使所述双对双差的第二权重大于所述台站对双差的第一权重,直到迭代次数达到第二预设次数,迭代结束,得到最终各个所述地震事件的地震位置,其中所述地震位置包括:地震绝对位置以及地震间的相对位置;
其中,所述构建单元至少包括:
设置子单元,用于为所述台站对双差设置第一权重,为所述双对双差设置第二权重。
6.根据权利要求5所述的装置,其特征在于,所述第一计算单元包括:
第一计算子单元,用于根据所述地震事件对(i,j)中地震事件i和地震事件j到台站k的到时和的观测值与预测值的差,得到所述地震事件i和所述地震事件j到台站k的到时残差和并根据所述地震事件i和所述地震事件j到所述台站k的到时残差和得到所述地震事件对(i,j)到所述台站k的事件对双差所述地震事件对(i,j)为多个地震事件对中的任意一个地震事件对,
第二计算子单元,用于根据所述地震事件i和所述地震事件j到台站l的到时Tl i和Tl j的观测值与预测值的差,得到所述地震事件i和所述地震事件j到台站l的到时残差rl i和rl j,并根据所述地震事件i和所述地震事件j到所述台站l的到时残差rl i和rl j,得到所述地震事件对(i,j)到所述台站l的事件对双差drl ij,
其中,dx为地震位置的变化量,dτ为发震时刻的变化量,(Tl i-Tl j)obs为观测的事件对(i,j)到时差数据,(Tl i-Tl j)cal为预测的事件对(i,j)到时差数据。
7.根据权利要求5所述的装置,其特征在于,所述第二计算单元包括:
第三计算子单元,用于根据所述地震事件i到所述台站k和所述台站l的到时和Tl i的观测值与预测值的差,得到所述地震事件i到所述台站k和所述台站l的到时残差和rl i,并根据所述地震事件i到所述台站k和所述台站l的到时残差和rl i,得到地震i到所述台站对(k,l)的台站对双差所述台站对(k,l)为多个台站对中的任意一个台站对;
第四计算子单元,用于根据所述地震事件j到所述台站k和所述台站l的到时和Tl j的观测值与预测值的差,得到所述地震事件j到所述台站k和所述台站l的到时残差和rl j并根据所述地震事件j到所述台站k和所述台站l的到时残差和rl j,得到地震j到所述台站对(k,l)的台站对双差
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611192581.9A CN106772600B (zh) | 2016-12-21 | 2016-12-21 | 一种双对双差地震定位方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611192581.9A CN106772600B (zh) | 2016-12-21 | 2016-12-21 | 一种双对双差地震定位方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106772600A CN106772600A (zh) | 2017-05-31 |
CN106772600B true CN106772600B (zh) | 2020-08-28 |
Family
ID=58893641
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611192581.9A Active CN106772600B (zh) | 2016-12-21 | 2016-12-21 | 一种双对双差地震定位方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106772600B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112731534A (zh) * | 2020-12-18 | 2021-04-30 | 重庆大学 | 一种考虑p波初至***误差的双声发射事件联合定位方法、***、电子设备及可读介质 |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111190740B (zh) * | 2020-01-02 | 2023-04-18 | 江南大学 | 一种用于波形互相关双差定位算法的并行优化方法 |
CN111257928B (zh) * | 2020-02-12 | 2021-08-13 | 中国科学技术大学 | 一种震源的定位方法及定位*** |
CN111352160B (zh) * | 2020-03-19 | 2020-11-10 | 中国科学院地质与地球物理研究所 | 一种海底地震仪自动重定位装置及方法 |
CN111580165A (zh) * | 2020-05-27 | 2020-08-25 | 中国科学院地质与地球物理研究所 | 一种海底地震仪到时差定位装置及方法 |
CN114200512A (zh) * | 2021-12-14 | 2022-03-18 | 四川旷谷信息工程有限公司 | 一种用于铁路施工期关键工点的地震烈度预警方法及*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105510880A (zh) * | 2014-09-23 | 2016-04-20 | 中国石油化工股份有限公司 | 一种基于双差法的微地震震源定位方法 |
CN105589100A (zh) * | 2014-10-21 | 2016-05-18 | 中国石油化工股份有限公司 | 一种微地震震源位置和速度模型同时反演方法 |
CN105759311A (zh) * | 2016-01-25 | 2016-07-13 | 西南交通大学 | 一种近实时地震震源位置定位方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10345463B2 (en) * | 2014-02-21 | 2019-07-09 | Ion Geophysical Corporation | Methods and systems for using known source events in seismic data processing |
-
2016
- 2016-12-21 CN CN201611192581.9A patent/CN106772600B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105510880A (zh) * | 2014-09-23 | 2016-04-20 | 中国石油化工股份有限公司 | 一种基于双差法的微地震震源定位方法 |
CN105589100A (zh) * | 2014-10-21 | 2016-05-18 | 中国石油化工股份有限公司 | 一种微地震震源位置和速度模型同时反演方法 |
CN105759311A (zh) * | 2016-01-25 | 2016-07-13 | 西南交通大学 | 一种近实时地震震源位置定位方法 |
Non-Patent Citations (2)
Title |
---|
Development of double-pair double difference earthquake location algorithm for improving earthquake locations;Hao Guo 等;《Geophysical Journal International》;20161020;第208卷(第1期);第333-348页 * |
基于长江安徽段***主动源实验庐江台阵的台站对双差地震速度成像;胡景 等;《中国地震》;20160630;第32卷(第2期);第343-355页 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112731534A (zh) * | 2020-12-18 | 2021-04-30 | 重庆大学 | 一种考虑p波初至***误差的双声发射事件联合定位方法、***、电子设备及可读介质 |
Also Published As
Publication number | Publication date |
---|---|
CN106772600A (zh) | 2017-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106772600B (zh) | 一种双对双差地震定位方法及装置 | |
Ebrahimian et al. | Bayesian optimal estimation for output‐only nonlinear system and damage identification of civil structures | |
CN105759311B (zh) | 一种近实时地震震源位置定位方法 | |
CN106054214B (zh) | 一种可削弱多路径效应的gnss数据处理方法 | |
CN103808349B (zh) | 矢量传感器的误差校正方法和装置 | |
CN111551895B (zh) | 基于加权多维标度和拉格朗日乘子的运动源tdoa和fdoa定位方法 | |
CN104066179A (zh) | 一种改进的自适应迭代ukf的wsn节点定位方法 | |
CN104992054B (zh) | 基于时间序列二维化的电离层垂直总电子含量预报方法 | |
Guasch | A direct transmissibility formulation for experimental statistical energy analysis with no input power measurements | |
CN102540252A (zh) | 基于互相关的高精度中值叠加方法 | |
CN108717201A (zh) | 一种隧道围岩微震源定位方法 | |
Chen et al. | Two‐stage damage detection algorithms of structure using modal parameters identified from recursive subspace identification | |
CN108802814A (zh) | 一种隧道围岩微震波速的获取方法 | |
CN103543026A (zh) | 基于振动传递率函数和支持向量机的结构损伤识别方法 | |
CN112098947B (zh) | 一种基于ctls的声发射源定位方法、***及存储介质 | |
Na et al. | Automated assessment of building damage from seismic events using smartphones | |
WO2020255796A1 (ja) | 受信端末、測位方法及びプログラム | |
CN110221245B (zh) | 联合估计目标位置和非视距误差的鲁棒tdoa定位方法 | |
CN111736208A (zh) | 变权重联合P波和S波初至数据的微震事件Bayes定位方法、***及介质 | |
CN109581504B (zh) | 地面微地震记录中直达纵波极性计算方法及*** | |
CN104597461A (zh) | 一种高效的glonass信号模拟器卫星星座仿真方法 | |
CN112130202B (zh) | 一种正交各向异性速度反演的方法及*** | |
CN112255647A (zh) | 弹道测量信息事后融合处理优化方法 | |
Canclini et al. | Distributed 3D source localization from 2D DOA measurements using multiple linear arrays | |
CN106100609A (zh) | 单状态变量和两级Kalman滤波器时间尺度算法 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |