CN116202570B - 一种多频率面波标定方法及装置 - Google Patents
一种多频率面波标定方法及装置 Download PDFInfo
- Publication number
- CN116202570B CN116202570B CN202310183603.9A CN202310183603A CN116202570B CN 116202570 B CN116202570 B CN 116202570B CN 202310183603 A CN202310183603 A CN 202310183603A CN 116202570 B CN116202570 B CN 116202570B
- Authority
- CN
- China
- Prior art keywords
- groups
- radial
- tangential
- data
- component
- 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 41
- 230000008878 coupling Effects 0.000 claims abstract description 62
- 238000010168 coupling process Methods 0.000 claims abstract description 62
- 238000005859 coupling reaction Methods 0.000 claims abstract description 62
- 238000005553 drilling Methods 0.000 claims abstract description 59
- 238000005259 measurement Methods 0.000 claims abstract description 56
- 230000033001 locomotion Effects 0.000 claims abstract description 34
- 238000001914 filtration Methods 0.000 claims description 9
- 238000004891 communication Methods 0.000 description 11
- 239000011435 rock Substances 0.000 description 8
- 238000010586 diagram Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 239000004568 cement Substances 0.000 description 4
- 239000007787 solid Substances 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 238000009434 installation Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 101100234408 Danio rerio kif7 gene Proteins 0.000 description 1
- 101100221620 Drosophila melanogaster cos gene Proteins 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 101100398237 Xenopus tropicalis kif11 gene Proteins 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 239000003989 dielectric material Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D21/00—Measuring or testing not otherwise provided for
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/32—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring the deformation in a solid
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C1/00—Measuring angles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P3/00—Measuring linear or angular speed; Measuring differences of linear or angular speeds
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/01—Measuring or predicting earthquakes
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Radar, Positioning & Navigation (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请实施例提供一种多频率面波标定方法及装置，包括：获取多组四分量钻孔应变数据和多组水平向震动数据；根据多组四分量钻孔应变数据和钻孔应变仪的方位角，计算多组间接测量值；根据多组水平向震动数据和地震波射线的方位角，计算多组径向地动速度和切向地动速度；根据多组间接测量值、多组径向地动速度和切向地动速度，计算耦合系数。依本申请的标定方法，能够准确测量耦合系数。
Description
技术领域
本申请实施例涉及测量技术领域，尤其涉及一种多频率面波标定方法及装置。
背景技术
钻孔应变仪通过耦合材料(例如，膨胀水泥)与井下围岩进行刚性耦合，测量岩石应力并转换为仪器所观测到的应变。由于钻孔应变测量的耦合效应，应变仪测得的线应变同时存在两个耦合系数，只有在耦合系数确定的前提下，应变仪的测量值才能开展地壳应力场的相关计算，否则其测量值毫无意义。耦合系数反映了应变仪的测量传感器与围岩的耦合程度，其大小与仪器安装和耦合材料固化过程有关，实验室条件下无法测量得到，一般利用理论固体潮信号标定，然而，由于理论固体潮与实际固体潮存在差异，且误差无法确定，因此，理论固体潮标定法所确定的耦合系数并不准确。
发明内容
有鉴于此，本申请实施例的目的在于提出一种多频率面波标定方法及装置，能够准确测量耦合系数。
基于上述目的，本申请实施例提供了一种多频率面波标定方法，包括：
获取多组四分量钻孔应变数据和多组水平向震动数据；
根据多组四分量钻孔应变数据和钻孔应变仪的方位角，计算多组间接测量值；
根据多组水平向震动数据和地震波射线的方位角，计算多组径向地动速度和切向地动速度；
根据多组间接测量值、多组径向地动速度和切向地动速度，计算耦合系数。
可选的，根据多组四分量钻孔应变数据和钻孔应变仪的方位角，计算多组间接测量值，方法为：
sa＝(S1+S2+S3+S4)/2 (3)
其中，S1、S2、S3、S4分别为四分量钻孔应变数据，θ1为钻孔应变仪的其中一个测量单元的方位角，Sa、S0和S45为间接测量值。
可选的，根据多组水平向震动数据和地震波射线的方位角，计算多组径向地动速度和切向地动速度，方法为：
其中，和/>分别为地理坐标系下地震仪测量的北南向分量的震动数据和东西向分量的震动数据，θr为地震波射线的方位角，/>和/>分别为射线坐标系下的径向地动速度和切向地动速度。
可选的，根据多组间接测量值、径向地动速度和切向地动速度，计算耦合系数，方法为：
sr0＝s0cos2θr+s45sin2θr (6)
sr45＝-s0sin2θr+s45cos2θr (7)
其中，cR和cL分别为观测点局部瑞利波和勒夫波的相速度，A、B为耦合系数。
可选的，根据多组间接测量值、多组径向地动速度和切向地动速度，计算耦合系数，包括：
分别对多组间接测量值、多组径向地动速度和切向地动速度进行窄带滤波处理，得到滤波后的多组间接测量值、径向地动速度和切向地动速度；
根据滤波后的多组间接测量值、径向地动速度和切向地动速度，计算相应的多组耦合系数；
计算多组耦合系数的均值作为最终的耦合系数。
可选的，所述获取多组四分量钻孔应变数据和多组水平向震动数据为：
获取钻孔应变仪的多组四分量钻孔应变数据和地震仪的多组水平向震动数据；其中，所述钻孔应变仪和地震仪以观测同一观测点的并列或共点方式安装。
可选的，所述获取多组四分量钻孔应变数据和多组水平向震动数据为：
获取钻孔应变仪的多组四分量钻孔应变数据和地震仪的多组水平向震动数据；其中，所述钻孔应变仪和地震仪之间的台站间距L满足L≤5％λ，λ为面波波长。
本申请实施例还提供一种多频率面波标定装置，包括：
获取模块，用于获取多组四分量钻孔应变数据和多组水平向震动数据；
第一计算模块，用于根据多组四分量钻孔应变数据和钻孔应变仪的方位角，计算多组间接测量值；
第二计算模块，用于根据多组水平向震动数据和地震波射线的方位角，计算多组径向地动速度和切向地动速度；
参数计算模块，用于根据多组间接测量值、多组径向地动速度和切向地动速度，计算耦合系数。
可选的，所述参数计算模块，用于分别对多组间接测量值、多组径向地动速度和切向地动速度进行窄带滤波处理，得到滤波后的多组间接测量值、径向地动速度和切向地动速度；根据滤波后的多组间接测量值、径向地动速度和切向地动速度，计算相应的多组耦合系数；计算多组耦合系数的均值作为最终的耦合系数。
可选的，所述获取模块，用于获取钻孔应变仪的多组四分量钻孔应变数据和地震仪的多组水平向震动数据；其中，所述钻孔应变仪和地震仪以观测同一观测点的并列或共点方式安装，或者，所述钻孔应变仪和地震仪之间的台站间距L满足L≤5％λ，λ为面波波长。
从上面所述可以看出，本申请实施例提供的多频率面波标定方法及装置，通过获取多组四分量钻孔应变数据和多组水平向震动数据，根据多组四分量钻孔应变数据和钻孔应变仪的方位角，计算多组间接测量值，根据多组水平向震动数据和地震波射线的方位角，计算多组径向地动速度和切向地动速度，根据多组间接测量值、多组径向地动速度和切向地动速度，计算耦合系数。依本申请的标定方法，能够准确测量耦合系数。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本申请实施例的方法流程示意图；
图2A为本申请实施例的钻孔应变仪和地震仪的共点安装结构示意图；
图2B为本申请实施例的钻孔应变仪和地震仪的非共点安装结构示意图；
图3为本申请实施例的钻孔应变观测平面应变张量的原理模型示意图；
图4为本申请实施例的钻孔应变仪的四个测量元件的分布示意图；
图5为本申请实施例的装置结构示意图；
图6为本申请实施例的电子设备结构示意图。
具体实施方式
为使本公开的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本公开进一步详细说明。
需要说明的是，除非另外定义，本申请实施例使用的技术术语或者科学术语应当为本公开所属领域内具有一般技能的人士所理解的通常意义。本申请实施例中使用的“第一”、“第二”以及类似的词语并不表示任何顺序、数量或者重要性，而只是用来区分不同的组成部分。“包括”或者“包含”等类似的词语意指出现该词前面的元件或者物件涵盖出现在该词后面列举的元件或者物件及其等同，而不排除其他元件或者物件。“连接”或者“相连”等类似的词语并非限定于物理的或者机械的连接，而是可以包括电性的连接，不管是直接的还是间接的。“上”、“下”、“左”、“右”等仅用于表示相对位置关系，当被描述对象的绝对位置改变后，则该相对位置关系也可能相应地改变。
如图1所示，本申请实施例提供的多频率面波标定方法，包括：
S101：获取多组四分量钻孔应变数据和多组水平向震动数据；
本实施例中，获取钻孔应变仪的四分量钻孔应变数据，其中，钻孔应变仪是埋置于钻孔中，用于测量岩体应变的仪器，如图4所示，钻孔应变仪配置有四个分量的测量单元S1,S2,S3,S4，四个测量单元分别用于测量各自所处方位方向的应变量。获取地震仪的水平向震动数据，地震仪是埋设于特定安装地点，用于测量地面震动的仪器，地震仪配置有东西向、南北向和垂直向三个分向的测量单元，其中，东西向和南北向测量单元为水平向测量单元，水平向测量单元的观测数据为水平向震动数据。
如图2A所示，一些方式中，钻孔应变仪与地震仪并列或共点布设，此种安装方式下，钻孔应变仪和地震仪紧密连接，用于观测同一观测点的地震波。此种安装方式下，获取钻孔应变仪的四个测量单元感测的四分量钻孔应变数据，获取地震仪的水平向测量单元的水平向震动数据。
如图2B所示，另一些方式中，钻孔应变仪与地震仪非共点安装，当满足地震仪与钻孔应变仪之间的台站间距远小于面波波长时，可以根据地震仪的水平向震动数据和钻孔应变仪的四分量钻孔应变数据标定耦合***。一些情况下，当地震仪与钻孔应变仪之间的台站间距L≤5％λ时，可以认为台站间距L远小于面波波长λ，其中，面波为瑞丽波或勒夫波，例如，周期为τ＝20s的瑞丽波，波速为v＝4.0km/s，该瑞丽波的波长为λ＝vτ＝80km，当台站间距L不大于4km时，可以根据地震仪和钻孔应变仪的观测数据标定耦合系数。
S102：根据多组四分量钻孔应变数据和钻孔应变仪的方位角，计算多组间接测量值；
本实施例中，从钻孔应变仪的四个分量的测量单元获取四分量钻孔应变数据后，根据四分量钻孔应变数据和钻孔应变仪的方位角，计算间接测量值，方法为：
sa＝(S1+S2+S3+S4)/2(3)
其中，S1、S2、S3、S4分别为四个测量单元感测的四分量钻孔应变数据，θ1为其中一个测量单元的方位角，为已知量，sa、s0和s45为间接测量值。
S103：根据多组水平向震动数据和地震波射线的方位角，计算多组径向地动速度和切向地动速度；
本实施例中，从地震仪的水平向测量单元获取水平向震动数据后，根据水平向震动数据和地震波射线的方位角，计算径向地动速度和切向地动速度，方法为：
其中，和/>分别为地理坐标系下地震仪测量的北南向分量的震动数据和东西向分量的震动数据，θr为地震波射线的方位角，可由震中经纬度和台站的经纬度确定，/>和/>分别为射线坐标系下的径向地动速度和切向地动速度。
S104：根据多组间接测量值、径向地动速度和切向地动速度，计算耦合系数。
本实施例中，根据钻孔应变仪观测的四分量钻孔应变数据计算出间接测量值、根据地震仪观测的水平向震动数据计算出径向地动速度和切向地动速度之后，根据间接测量值、径向地动速度和切向地动速度，计算两个耦合系数，方法为：
sr0＝s0cos2θr+s45sin2θr (6)
sr45＝-s0sin2θr+s45cos2θr (7)
其中，cR和cL分别为观测点局部瑞利波(Rayleigh)和勒夫波(Love)的相速度。
以下结合附图和具体实施例对本申请的方法原理进行说明。
对于区域范围横向均匀介质传播的近似平面波，即介质材料特性的任何横向变化的尺度超过了主导能量的波长，例如，对于2-40s的周期对应范围约为5-150km。此时，非均质性是“缓慢变化的”，在一点上的变形可以视为介质局部均匀，地震在某一观测点引起的地震动与动态应变量存在如下关系：
其中，ε′rr和ε′rt分别为根据摆式地震仪的观测数据求得的瑞利波和勒夫波的应变分量，是一个中间量，在已知相速度的情况下求得。对于近地表传播的远震面波，公式(8)、(9)所示关系式仍然成立。
当地震波射线的方位角θr为已知的情况下，射线坐标系下径向地动速度和切向地动速度/>可由地震仪的水平分量震动数据(包括北南分量的震动数据/>和东西分量的震动数据/>)通过公式(4)得到。
如图3所示，对于钻孔应变仪的四分量钻孔应变观测，四分量的测量单元的探头通过膨胀水泥与井下围岩耦合，岩石、膨胀水泥和钻孔应变仪的套筒组成一个三层复合双衬套结构，当足够远处有均匀水平主应变ε1和ε2时，套筒内壁在θ方向的孔径相对变化可表示为：
Sθ＝A(ε1+ε2)+B(ε1-ε2)cos2(θ-φ) (10)
式中，耦合系数A、B反映了钻孔应变仪的测量单元与围岩的耦合程度，φ表示水平主应变ε1的方位角，ε2的方位角为φ+90°。根据公式(10)，耦合系数A和B的物理意义是由钢简、水泥与围岩耦合造成的对应变测量值的修正系数，其中，耦合系数A是对面应变(ε1+ε2)测量值的修正系数，耦合系数B是对剪应变(ε1-ε2)测量值的修正系数。
如图4所示，对于钻孔应变仪的四个测量单元，相邻两个测量单元之间的夹角均为45°，若测量单元S1的方位角为θ1，则测量单元S2、S3、S4的方位角分别为θ1、四个测量单元的观测值表示为：
钻孔应变仪安装之后，其中一个测量单元的方位角θ1即为已知量。
根据公式(10)，按照三角函数公式变换，得到：
通过坐标变换，得到钻孔应变仪的测量值在极坐标系下的径向应变量εrr和切向应变量εrt，表示为：
则远震面波的瑞利波可根据径向应变量εrr所确定，勒夫波可根据切向应变量εrt所确定。
当摆式地震仪和钻孔应变仪观测的面波应变分量两两相等，即ε′rr＝εrr，ε′rt＝εrt，根据公式(4)、(8)、(9)、(13)，可求解耦合系数A和B，表示为：
即在已知观测点局部介质相速度的情况下，可由地震仪记录的勒夫波速度(切向地动速度)标定耦合参数B，由瑞利波速度(径向地动速度/>)标定耦合参数A。进一步对公式(14)进行简化，得到公式(5)所示耦合参数的计算公式。
本实施例中，为提高耦合参数的准确性，采用多频率点数据进行面波标定。持续采集多组四分量钻孔应变数据和多组水平向震动数据；根据多组四分量钻孔应变数据计算出对应的多组间接测量值，根据多组水平向震动数据计算出对应的多组径向地动速度和切向地动速度；由于采集的多组四分量钻孔应变数据和多组水平向震动数据频率可能相同也可能不同，因而计算出的多组间接测量值频率可能相同也可能不同，计算出的多组径向地动速度和切向地动速度频率可能相同也可能不同；基于多组频率相同或不同的间接测量值、径向地动速度和切向地动速度，分别进行窄带滤波处理，基于滤波处理后的多组间接测量值、径向地动速度和切向地动速度，利用公式(5)计算出对应的多组耦合系数；最后基于多组耦合系数计算耦合系数均值，将耦合系数均值作为最终确定的耦合系数。这样，基于多样本、多频率点数据的统计标定，能够获得较为准确的耦合系数结果。
一些方式中，对多组、多频率点的间接测量值、径向地动速度和切向地动速度，进行窄带滤波的窄带滤波器的参数为，频率中心点f0i，带通滤波截至频率fHi和fLi，且满足以下关系：
对应不同的中心频率点(f0i，i＝1-n,n个滤波器)设计一系列的窄带滤波器，并对观测分量sr0、sr45和进行窄带滤波。可根据目标研究范围来设定滤波器的频率中心点的数量，例如，目标频带范围为1/20Hz-1/200Hz，可将该频带分为多个子频段，并设计多个频点为滤波器的中心频率。
需要说明的是，本申请实施例的方法可以由单个设备执行，例如一台计算机或服务器等。本实施例的方法也可以应用于分布式场景下，由多台设备相互配合来完成。在这种分布式场景的情况下，这多台设备中的一台设备可以只执行本申请实施例的方法中的某一个或多个步骤，这多台设备相互之间会进行交互以完成所述的方法。
需要说明的是，上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
如图5所示，本申请实施例还提供一种多频率面波标定装置，包括：
获取模块，用于获取多组四分量钻孔应变数据和多组水平向震动数据；
第一计算模块，用于根据多组四分量钻孔应变数据和钻孔应变仪的方位角，计算多组间接测量值；
第二计算模块，用于根据多组水平向震动数据和地震波射线的方位角，计算多组径向地动速度和切向地动速度；
参数计算模块，用于根据多组间接测量值、多组径向地动速度和切向地动速度，计算耦合系数。
为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本申请实施例时可以把各模块的功能在同一个或多个软件和/或硬件中实现。
上述实施例的装置用于实现前述实施例中相应的方法，并且具有相应的方法实施例的有益效果，在此不再赘述。
图6示出了本实施例所提供的一种更为具体的电子设备硬件结构示意图，该设备可以包括：处理器1010、存储器1020、输入/输出接口1030、通信接口1040和总线1050。其中处理器1010、存储器1020、输入/输出接口1030和通信接口1040通过总线1050实现彼此之间在设备内部的通信连接。
处理器1010可以采用通用的CPU(Central Processing Unit，中央处理器)、微处理器、应用专用集成电路(Application Specific Integrated Circuit，ASIC)、或者一个或多个集成电路等方式实现，用于执行相关程序，以实现本说明书实施例所提供的技术方案。
存储器1020可以采用ROM(Read Only Memory，只读存储器)、RAM(Random AccessMemory，随机存取存储器)、静态存储设备，动态存储设备等形式实现。存储器1020可以存储操作***和其他应用程序，在通过软件或者固件来实现本说明书实施例所提供的技术方案时，相关的程序代码保存在存储器1020中，并由处理器1010来调用执行。
输入/输出接口1030用于连接输入/输出模块，以实现信息输入及输出。输入输出/模块可以作为组件配置在设备中(图中未示出)，也可以外接于设备以提供相应功能。其中输入设备可以包括键盘、鼠标、触摸屏、麦克风、各类传感器等，输出设备可以包括显示器、扬声器、振动器、指示灯等。
通信接口1040用于连接通信模块(图中未示出)，以实现本设备与其他设备的通信交互。其中通信模块可以通过有线方式(例如USB、网线等)实现通信，也可以通过无线方式(例如移动网络、WIFI、蓝牙等)实现通信。
总线1050包括一通路，在设备的各个组件(例如处理器1010、存储器1020、输入/输出接口1030和通信接口1040)之间传输信息。
需要说明的是，尽管上述设备仅示出了处理器1010、存储器1020、输入/输出接口1030、通信接口1040以及总线1050，但是在具体实施过程中，该设备还可以包括实现正常运行所必需的其他组件。此外，本领域的技术人员可以理解的是，上述设备中也可以仅包含实现本说明书实施例方案所必需的组件，而不必包含图中所示的全部组件。
上述实施例的电子设备用于实现前述实施例中相应的方法，并且具有相应的方法实施例的有益效果，在此不再赘述。
本实施例的计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。
所属领域的普通技术人员应当理解：以上任何实施例的讨论仅为示例性的，并非旨在暗示本公开的范围(包括权利要求)被限于这些例子；在本公开的思路下，以上实施例或者不同实施例中的技术特征之间也可以进行组合，步骤可以以任意顺序实现，并存在如上所述的本申请实施例的不同方面的许多其它变化，为了简明它们没有在细节中提供。
另外，为简化说明和讨论，并且为了不会使本申请实施例难以理解，在所提供的附图中可以示出或可以不示出与集成电路(IC)芯片和其它部件的公知的电源/接地连接。此外，可以以框图的形式示出装置，以便避免使本申请实施例难以理解，并且这也考虑了以下事实，即关于这些框图装置的实施方式的细节是高度取决于将要实施本申请实施例的平台的(即，这些细节应当完全处于本领域技术人员的理解范围内)。在阐述了具体细节(例如，电路)以描述本公开的示例性实施例的情况下，对本领域技术人员来说显而易见的是，可以在没有这些具体细节的情况下或者这些具体细节有变化的情况下实施本申请实施例。因此，这些描述应被认为是说明性的而不是限制性的。
尽管已经结合了本公开的具体实施例对本公开进行了描述，但是根据前面的描述，这些实施例的很多替换、修改和变型对本领域普通技术人员来说将是显而易见的。例如，其它存储器架构(例如，动态RAM(DRAM))可以使用所讨论的实施例。
本申请实施例旨在涵盖落入所附权利要求的宽泛范围之内的所有这样的替换、修改和变型。因此，凡在本申请实施例的精神和原则之内，所做的任何省略、修改、等同替换、改进等，均应包含在本公开的保护范围之内。
Claims (7)
1.一种多频率面波标定方法，其特征在于，包括：
获取多组四分量钻孔应变数据和多组水平向震动数据；
根据多组四分量钻孔应变数据和钻孔应变仪的方位角，计算多组间接测量值，方法为：
sa＝(S1+S2+S3+S4)/2 (3)
其中，S1、S2、S3、S4分别为四分量钻孔应变数据，θ1为钻孔应变仪的其中一个测量单元的方位角，Sa、S0和S45为间接测量值；
根据多组水平向震动数据和地震波射线的方位角，计算多组径向地动速度和切向地动速度，方法为：
其中，和/>分别为地理坐标系下地震仪测量的北南向分量的震动数据和东西向分量的震动数据，θr为地震波射线的方位角，/>和/>分别为射线坐标系下的径向地动速度和切向地动速度；
根据多组间接测量值、多组径向地动速度和切向地动速度，计算耦合系数，方法为：
sr0＝s0cos2θr+s45sin2θr (6)
sr45＝-s0sin2θr+s45cos2θr (7)
其中，cR和cL分别为观测点局部瑞利波和勒夫波的相速度，A、B为耦合系数。
2.根据权利要求1所述的方法，其特征在于，根据多组间接测量值、多组径向地动速度和切向地动速度，计算耦合系数，包括：
分别对多组间接测量值、多组径向地动速度和切向地动速度进行窄带滤波处理，得到滤波后的多组间接测量值、径向地动速度和切向地动速度；
根据滤波后的多组间接测量值、径向地动速度和切向地动速度，计算相应的多组耦合系数；
计算多组耦合系数的均值作为最终的耦合系数。
3.根据权利要求1所述的方法，其特征在于，所述获取多组四分量钻孔应变数据和多组水平向震动数据为：
获取钻孔应变仪的多组四分量钻孔应变数据和地震仪的多组水平向震动数据；其中，所述钻孔应变仪和地震仪以观测同一观测点的并列或共点方式安装。
4.根据权利要求1所述的方法，其特征在于，所述获取多组四分量钻孔应变数据和多组水平向震动数据为：
获取钻孔应变仪的多组四分量钻孔应变数据和地震仪的多组水平向震动数据；其中，所述钻孔应变仪和地震仪之间的台站间距L满足L≤5％λ，λ为面波波长。
5.一种多频率面波标定装置，其特征在于，包括：
获取模块，用于获取多组四分量钻孔应变数据和多组水平向震动数据；
第一计算模块，用于根据多组四分量钻孔应变数据和钻孔应变仪的方位角，计算多组间接测量值，方法为：
sa＝(S1+S2+S3+S4)/2 (3)
其中，S1、S2、S3、S4分别为四分量钻孔应变数据，θ1为钻孔应变仪的其中一个测量单元的方位角，Sa、S0和S45为间接测量值；
第二计算模块，用于根据多组水平向震动数据和地震波射线的方位角，计算多组径向地动速度和切向地动速度，方法为：
其中，和/>分别为地理坐标系下地震仪测量的北南向分量的震动数据和东西向分量的震动数据，θr为地震波射线的方位角，/>和/>分别为射线坐标系下的径向地动速度和切向地动速度；
参数计算模块，用于根据多组间接测量值、多组径向地动速度和切向地动速度，计算耦合系数，方法为：
sr0＝s0cos2θr+s45sin2θr (6)
sr45＝-s0sin2θr+s45cos2θr (7)
其中，cR和cL分别为观测点局部瑞利波和勒夫波的相速度，A、B为耦合系数。
6.根据权利要求5所述的装置，其特征在于，
所述参数计算模块，用于分别对多组间接测量值、多组径向地动速度和切向地动速度进行窄带滤波处理，得到滤波后的多组间接测量值、径向地动速度和切向地动速度；根据滤波后的多组间接测量值、径向地动速度和切向地动速度，计算相应的多组耦合系数；计算多组耦合系数的均值作为最终的耦合系数。
7.根据权利要求5所述的装置，其特征在于，
所述获取模块，用于获取钻孔应变仪的多组四分量钻孔应变数据和地震仪的多组水平向震动数据；其中，所述钻孔应变仪和地震仪以观测同一观测点的并列或共点方式安装，或者，所述钻孔应变仪和地震仪之间的台站间距L满足L≤5％λ，λ为面波波长。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310183603.9A CN116202570B (zh) | 2023-02-17 | 2023-02-17 | 一种多频率面波标定方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310183603.9A CN116202570B (zh) | 2023-02-17 | 2023-02-17 | 一种多频率面波标定方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116202570A CN116202570A (zh) | 2023-06-02 |
| CN116202570B true CN116202570B (zh) | 2023-09-05 |
Family
ID=86518876
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310183603.9A Active CN116202570B (zh) | 2023-02-17 | 2023-02-17 | 一种多频率面波标定方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116202570B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN202083345U (zh) * | 2011-03-01 | 2011-12-21 | 欧阳祖熙 | 一种用于钻孔应变仪的位移传感装置 |
| CN106918836A (zh) * | 2017-03-31 | 2017-07-04 | 吉林大学 | 基于主成分分析的钻孔应变数据异常提取方法 |
| CN111141434A (zh) * | 2019-12-23 | 2020-05-12 | 中国科学院大学 | 一种基于应力花瓣图确定四分量钻孔应力变化的方法 |
| CN111812746A (zh) * | 2020-07-23 | 2020-10-23 | 中国地震局地壳应力研究所 | 基于地震表面波的钻孔应变仪原位标定方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1774364A2 (en) * | 2004-05-21 | 2007-04-18 | Entre Holdings Company | Full wave seismic recording system |
-
2023
- 2023-02-17 CN CN202310183603.9A patent/CN116202570B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN202083345U (zh) * | 2011-03-01 | 2011-12-21 | 欧阳祖熙 | 一种用于钻孔应变仪的位移传感装置 |
| CN106918836A (zh) * | 2017-03-31 | 2017-07-04 | 吉林大学 | 基于主成分分析的钻孔应变数据异常提取方法 |
| CN111141434A (zh) * | 2019-12-23 | 2020-05-12 | 中国科学院大学 | 一种基于应力花瓣图确定四分量钻孔应力变化的方法 |
| CN111812746A (zh) * | 2020-07-23 | 2020-10-23 | 中国地震局地壳应力研究所 | 基于地震表面波的钻孔应变仪原位标定方法 |
Non-Patent Citations (1)
| Title |
|---|
| 分量式钻孔应变观测耦合系数的计算;张凌空 等;《地球物理学报》;第56卷(第9期);第3029-3037页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116202570A (zh) | 2023-06-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bernauer et al. | Rotational sensors—A comparison of different sensor types | |
| CN101290326B (zh) | 石英挠性加速度计测量组件的参数辨识标定方法 | |
| CN105717553B (zh) | 一种旋转加速度计重力梯度仪标定方法 | |
| Bonvalot et al. | Continuous gravity recording with Scintrex CG-3M meters: a promising tool for monitoring active zones | |
| CN102590553B (zh) | 一种基于小波消噪的加速度计温度补偿方法 | |
| Trifunac et al. | A note on the accuracy of computed ground displacements from strong-motion accelerograms | |
| Scudero et al. | MEMS technology in seismology: A short review | |
| CN113885098B (zh) | 一种重力敏感器低频频率响应误差在线建模及补偿方法 | |
| CN104459826A (zh) | 旋转加速度计重力梯度仪重力梯度信号仿真方法 | |
| CN106768179A (zh) | 基于连续运行gnss站信噪比数据的潮位的测量方法 | |
| Zang et al. | Real-time coseismic deformations from adaptively tight integration of high-rate GNSS and strong motion records | |
| Pucciarelli | Wavelet analysis in volcanology: The case of phlegrean fields | |
| CN116202570B (zh) | 一种多频率面波标定方法及装置 | |
| CN111812746B (zh) | 基于地震表面波的钻孔应变仪原位标定方法 | |
| CN103868527A (zh) | 一种标定捷联惯性组合加速度计组合的方法 | |
| CN201417314Y (zh) | 一种三分量光电混合集成加速度地震检波器 | |
| RU2352961C2 (ru) | Способ определения пространственного положения и параметров движения внутреннего ядра земли | |
| CN105509871A (zh) | 一种应用于振动传感器的自检定装置及其自检定方法 | |
| Kong et al. | A detection method of earthquake precursory anomalies using the four-component borehole strainmeter | |
| Abramov et al. | Key elements of construction and exploitation of a land-based seismic-gravimetric complex | |
| CN208887684U (zh) | 一种具有mems传感器的应力及应变测试仪 | |
| RU2337382C1 (ru) | Способ краткосрочного предсказания землетрясений | |
| Chen et al. | Six‐Component Earthquake Synchronous Observations Across Taiwan Strait: Phase Velocity and Source Location | |
| Dergach et al. | Absolute calibration of seismic sensors using a displacement jump: Theory and practice | |
| RU2581076C2 (ru) | Способ и устройство дистанционного гравиметрического зондирования |
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 |