CN107240228A - 一种基于倾斜仪和gps的火山监测与预警系统及方法 - Google Patents

一种基于倾斜仪和gps的火山监测与预警系统及方法 Download PDF

Info

Publication number
CN107240228A
CN107240228A CN201710370031.XA CN201710370031A CN107240228A CN 107240228 A CN107240228 A CN 107240228A CN 201710370031 A CN201710370031 A CN 201710370031A CN 107240228 A CN107240228 A CN 107240228A
Authority
CN
China
Prior art keywords
mrow
mtd
msub
mtr
gps
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
CN201710370031.XA
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.)
East China Jiaotong University
Original Assignee
East China Jiaotong University
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 East China Jiaotong University filed Critical East China Jiaotong University
Priority to CN201710370031.XA priority Critical patent/CN107240228A/zh
Publication of CN107240228A publication Critical patent/CN107240228A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G08—SIGNALLING
    • G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B21/00—Alarms responsive to a single specified undesired or abnormal condition and not otherwise provided for
    • G08B21/02—Alarms for ensuring the safety of persons
    • G08B21/10—Alarms for ensuring the safety of persons responsive to calamitous events, e.g. tornados or earthquakes
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B7/00—Measuring arrangements characterised by the use of electric or magnetic techniques
    • G01B7/16—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring the deformation in a solid, e.g. by resistance strain gauge
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C9/00—Measuring inclination, e.g. by clinometers, by levels
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position
    • G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Business, Economics & Management (AREA)
  • Emergency Management (AREA)
  • Environmental & Geological Engineering (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明公开了一种基于倾斜仪和GPS的火山监测与预警系统及方法，系统包括倾斜测量子系统、GPS监测子系统、地控预警子系统；倾斜测量子系统，用于根据地表结构因受到岩浆的压力或质量再分布而发生形变而引起的倾斜测量子系统形变，获得的地表倾斜角度，从而判断火山地表运动的微小形变，进行预警；GPS监测子系统，用于获取原始观测数据，判断火山运动引起的地表位移量，进行相应的火山运动状态监测；地控预警子系统，用于接收倾斜测量子系统和GPS监测子系统传输的倾斜测量数据和GPS数据，并进行数据处理，分析得到预警数据，并根据设置的预警值，发出预警信号，从而实现对火山状态进行检测。本发明对火山的监测更加准确，保证了数据的可靠性。

Description

一种基于倾斜仪和GPS的火山监测与预警系统及方法
技术领域
本发明属于灾害预警技术领域，涉及一种地质灾害预警系统及方法，具体涉及一种基于倾斜仪和GPS的火山监测与预警系统及方法。
背景技术
现代活动的四个火山岩带主要是环太平洋火山带、大西洋火山带、地中海火山带和东非火山带，对这几个火山带的监测安装了大量的传感器，设置诸多的地震站，如圣海伦斯火山周围有13个，夏威夷基拉韦亚火山周围有47个，印尼默拉皮火山周围设6个，苏弗里埃尔火山在1978年4月大爆发前，可感地震每小时达15次，火山活动的过程常造成许多微小地震，大爆发更能引起强烈地震。地震的发生也常导致火山活动。目前，对火山的监测有诸多的方案，向生建等人利用地震监测设备、山体滑坡监测设备、泥石流监测设备等，连接客户端与报警器，远程管理监测数据；柯福阳利用采集每个监测站点和连续运行的参考站以及周边的观测数据分析处理进行预警；张科元根据地震周期和现有监测地震和火山的监测数据在地壳内压力增大时人工引爆地球上现有的休眠火山小规模分批喷发排放地球内部的高压物质以减轻地球内部压力；但是，上述等人的监测系统，需求设备量大，数据的精准度不高，需求人工人力的支持量强度大，未考虑前期火山微小地表的运动状态的变化，基于此，针对现有对火山监测和预警系统的现状，需要提供一种低成本、便利又测量精度较高的火山监测预警系统及方法。
发明内容
为了解决上述技术问题，本发明提供了一种基于倾斜仪和GPS的火山监测与预警系统及方法。
本发明的系统所采用的技术方案是：一种基于倾斜仪和GPS的火山监测与预警系统，其特征在于：包括倾斜测量子系统、GPS监测子系统、地控预警子系统；
所述倾斜测量子系统，用于根据地表结构因受到岩浆的压力或质量再分布而发生形变而引起的倾斜测量子系统形变，获得的地表倾斜角度，从而判断火山地表运动的微小形变，并根据预设的限值，进行预警；
所述GPS监测子系统，用于获取原始观测数据，判断火山运动引起的地表位移量，进行相应的火山运动状态监测；
所述地控预警子系统，用于接收倾斜测量子系统和GPS监测子系统传输的倾斜测量数据和GPS数据，并进行数据处理，分析得到预警数据，并根据设置的预警值，发出预警信号，从而实现对火山状态进行检测。
本发明的方法所采用的技术方案是：一种基于倾斜仪和GPS的火山监测与预警方法，其特征在于，包括以下步骤：
步骤1：测量地表倾斜度；
根据地表结构因受到岩浆的压力或质量再分布而发生形变而引起的倾斜仪形变，获得的地表倾斜角度θ；
步骤2：测量地表位移量；
利用设置在火山口两侧的GPS接收机，获取GPS基线向量、载波相位观测值、星历历元数据；对星历历元数据进行GPSRTK处理，组设的载波相位双差观测方程以及双差观测方程的误差方程模型，精准处理观测数据，实时观测数据实时处理数据，判断基线长以及地表监测点的三维坐标是否有变形，坐标位移量是否在预警数据精度范围；
步骤3：预警；
接收处理后的倾斜测量子系统和GPS监测子系统传输的倾斜测量数据和GPS数据，分析得到预警数据，并根据设置的预警值，发出预警信号；
如果不在预警数据范围内，则对当前倾斜测量子系统和GPS监测子系统观测值进行最小二乘估计，然后重复迭代，全天候、实时的连续进行观测。
相对于现有技术，本发明的有益效果是：本发明对数据处理严格把关，消除粗差，对数据的高精度、高要求使得对火山的监测更加准确，保证了数据的可靠性，实现对火山的全天候监测和数据的实时更新，通过对数据的分析能够及时对危险信息发出预报。
附图说明
图1是本发明实施例的倾斜仪结构示意图；
图2是本发明实施例的数据传输流程图；
图3是本发明实施例的方法流程图；
图4是本发明实施例的数据处理流程图。
具体实施方式
为了便于本领域普通技术人员理解和实施本发明，下面结合附图及实施例对本发明作进一步的详细描述，应当理解，此处所描述的实施示例仅用于说明和解释本发明，并不用于限定本发明。
本发明提供的一种基于倾斜仪和GPS的火山监测与预警系统，包括倾斜测量子系统、GPS监测子系统、地控预警子系统；
倾斜测量子系统主要是由布设在火山口附近的倾斜仪组成，利用对垂向位移变化比较敏感，灵敏度高的优势进行观测。倾斜测量系统的基本原理如下：当火山内部的物质(如岩浆等)在不规则的运动，或内部压力促使岩浆向火山口流动，监测站附近地表结构因受到岩浆的压力或质量再分布而发生形变，从而使得布设在火山口附近的倾斜仪发生形变，通过高灵敏度的倾斜仪获得的地倾斜角度θ，来判断火山地表运动的微小形变，并根据预设的限值，进行预警。
请见图1，本实施例的倾斜仪能够及时的观测到火山监测站点的结构物的变形，ELT-10型倾斜仪是由倾斜传感器、安装支架、信号传输电缆等组成，结构物产生的倾斜变形，通过安装支架传递给倾斜传感器，传感器内装有电解液和导电触点，当传感器发生倾斜变化时，电解液的液面始终处于水平，但液面相对触点的部位发生了改变，也同时引起了输出电量的改变，倾斜仪随结构物的倾斜变形量与输出的电量呈对应关系，以此可测出被测结构物的倾斜角度θ，同时它的测量值可显示出以零点为基准值的倾斜角变化的正负方向，倾斜仪可以回收重复使用，并且可方便实现倾斜测量的自动化。
GPS监测子系统由设置在火山口两侧的GPS接收机组成，观测所得的数据包括GPS基线向量，载波相位观测值、星历历元数据等，在地控服务中心的数据处理模块对GPS导航系统获得的星历数据进行GPSRTK处理，组设的载波相位双差观测方程以及观测误差模型，精准处理观测数据，实时观测数据实时处理数据，判断基线长以及地表监测点的三维坐标是否有变形，坐标位移量是否在预警数据精度范围，数据之间通过GPRS及时传输到数据处理模块对数据进行接收、储存。GPS的全天候的测量优势大大减轻了人力的消耗，以免因观测的缺失而导致的危险。
地控预警子系统包括地控服务中心和动态监控系统中心。地控服务中心通过DTU数据传输模块，接收传输的倾斜测量数据和GPS数据，利用实时数据处理模块进行高精度数据处理；动态监控系统中心利用接收到的数据，分析得到预警数据，并根据设置的预警值，发出预警信号，从而实现对火山状态进行检测。
DTU数据传输模块，通过GPRS远程传输数据，GPRS无线通信模块主要是起数据之间的无线传输作用，系统中应用无线网络进行各种数据的传输以及工作上的交流，GPRS功能强大，GPRS终端设备永久在线，无线传输数据的延时为秒级范围，GPRS-DTU与数据中心的数据通信速率一般在10kbps-60kbps之间，按照系统要求，设置数据平均通信量在30kbps以内(不同时段的要求根据情况而设定。)传输速率快，极大的方便和满足了工作网络需求。
请见图2，本实施例的DTU数据传输模块传输的观测数据主要有卫星数据、GPS接收输出的数据、倾斜仪观测数据，利用DTU数据传输模块将数据实时传输，DTU数据传输模块内部集成TCP/IP协议栈，封装了PPP拨号协议以及TCP/IP协议栈并且具有嵌入式操作系统，提供了串行通信接口，在设计上将串口数据设计成“透明转换”的方式，即将串口上的原始数据转换成TCP/IP数据包进行传送，不需要改变原有的数据通信内容，也就是将未经处理的原始数据经过GPRS传输到数据处理中心，做处理，经过一系列的数据处理之后再次利用无线传输网络将数据传输到预警系统接收，数据传输各自部件有各自的传输接口，根据数据量的大小，对不同的数据在传输的同时并行存储，且GPRS终端设备永久在线，确保数据的实时传输不断。
请见图3和图4，本发明提供的一种基于倾斜仪和GPS的火山监测与预警方法，包括以下步骤：
步骤1：测量地表倾斜度；
根据地表结构因受到岩浆的压力或质量再分布而发生形变而引起的倾斜仪形变，获得的地表倾斜角度θ；
ELT-10型的倾斜仪测量范围为±15°，灵敏度小于9(〃/F)测量精度为±0.1％(F.S)，耐水压大于等于1(MPa)，绝缘电阻大雨等于50(MΩ)，当被测结构物体发生倾斜变形时，其倾斜角度θ与输出的电量读数F可用如下计算公式：
θ＝a+b*F+c*F+d*F；
式中：θ为被测结构物的倾斜角度，单位为度(°)；F为倾斜仪的实时电量测量值，单位为F；a﹑b﹑c﹑d为倾斜仪的标定系数。
步骤2：测量地表位移量；
利用设置在火山口两侧的GPS接收机，获取GPS基线向量、载波相位观测值、星历历元数据；对星历历元数据进行GPSRTK处理，组设的载波相位双差观测方程以及双差观测方程的误差方程模型，精准处理观测数据，实时观测数据实时处理数据，判断基线长以及地表监测点的三维坐标是否有变形，坐标位移量是否在预警数据精度范围；
GPS基线向量的估值为：
mi-1是函数独立的基线向量，bi,k为第K条基线向量；
基线向量估值的验后方差-协方差阵为：
是基线向量k、l间的协方差子阵；对GPS接收机获得的星历历元数据进行GPSRTK处理，组设的载波相位双差观测方程为：
式中是两测站两卫星之间的双差载波相位，是两测站两卫星间双差整周模糊度参数，是两测站两卫星间双差星站距离，ε是测量噪声，c为光速，f是信号频率；
双差观测方程的误差方程模型为：
式中vk为观测值误差，δx2、δy2、δz2是监测站k的坐标近似值向量的改正数向量，为系数矩阵，ΔΔlk为监测站K的常数项误差。
步骤3：预警；
接收处理后的倾斜测量子系统和GPS监测子系统传输的倾斜测量数据和GPS数据，分析得到预警数据，并根据设置的预警值，发出预警信号；
如果不在预警数据范围内，则对当前倾斜测量子系统和GPS监测子系统观测值进行最小二乘估计，然后重复迭代，全天候、实时的连续进行观测。
本实施例的系统包括倾斜测量子系统、GPS监测子系统、地控预警子系统。倾斜测量主要是利用对垂向位移变化比较敏感，灵敏度高的优势进行观测，获取地倾斜角度θ，判断火山微小地表的运动状态，并根据预设的限值，进行预警；GPS监测子系统由设置在火山口两侧的GPS接收机组成。地控预警子系统包括地控服务中心和动态监控系统中心；系统通过对布设在火山扣附近的倾斜仪装置获取相应的观测资料即获取倾斜仪的倾斜角度θ来判断监测点是否变形，另一方面，通过对安置在火山口两侧的GPS接收机获得的观测数据进行实时解算，获得相应的基线向量，通过基线向量的变化值来判断火山运动引起的地表位移量，即通过基线向量的三维坐标变化量是否在预警数据范围内，进行相应的火山运动状态监测；三大子系统利用GPRS无线通信模块，将所有的测量数据经数据处理模块处理后的平差数据无线传输给预警系统的地控服务中心，经过系统分析评估，达到数据预警精度范围，则由此发出预警信号。如果不在预警数据范围内，则对当前观测值进行最小二乘估计，然后重复由获得的原始数据的流程继续迭代进行，全天候、实时的连续进行观测。
应当理解的是，本说明书未详细阐述的部分均属于现有技术。
应当理解的是，上述针对较佳实施例的描述较为详细，并不能因此而认为是对本发明专利保护范围的限制，本领域的普通技术人员在本发明的启示下，在不脱离本发明权利要求所保护的范围情况下，还可以做出替换或变形，均落入本发明的保护范围之内，本发明的请求保护范围应以所附权利要求为准。

Claims (7)

1.一种基于倾斜仪和GPS的火山监测与预警系统，其特征在于：包括倾斜测量子系统、GPS监测子系统、地控预警子系统；
所述倾斜测量子系统，用于根据地表结构因受到岩浆的压力或质量再分布而发生形变而引起的倾斜测量子系统形变，获得的地表倾斜角度，从而判断火山地表运动的微小形变，并根据预设的限值，进行预警；
所述GPS监测子系统，用于获取原始观测数据，判断火山运动引起的地表位移量，进行相应的火山运动状态监测；
所述地控预警子系统，用于接收倾斜测量子系统和GPS监测子系统传输的倾斜测量数据和GPS数据，并进行数据处理，分析得到预警数据，并根据设置的预警值，发出预警信号，从而实现对火山状态进行检测。
2.根据权利要求1所述的基于倾斜仪和GPS的火山监测与预警系统，其特征在于：所述倾斜测量子系统是由若干布设在火山口附近的倾斜仪组成。
3.根据权利要求1所述的基于倾斜仪和GPS的火山监测与预警系统，其特征在于：所述GPS监测子系统由设置在火山口两侧的GPS接收机组成。
4.根据权利要求1所述的基于倾斜仪和GPS的火山监测与预警系统，其特征在于：所述地控预警子系统包括地控服务中心和动态监控系统中心；
所述地控服务中心通过DTU数据传输模块，用于接收倾斜测量子系统和GPS监测子系统传输的倾斜测量数据和GPS数据，并进行数据处理；所述动态监控系统中心利用接收到的数据，分析得到预警数据，并根据设置的预警值，发出预警信号，从而实现对火山状态进行检测。
5.一种基于倾斜仪和GPS的火山监测与预警方法，其特征在于，包括以下步骤：
步骤1：测量地表倾斜度；
根据地表结构因受到岩浆的压力或质量再分布而发生形变而引起的倾斜仪形变，获得的地表倾斜角度θ；
步骤2：测量地表位移量；
利用设置在火山口两侧的GPS接收机，获取GPS基线向量、载波相位观测值、星历历元数据；对星历历元数据进行GPSRTK处理，组设的载波相位双差观测方程以及双差观测方程的误差方程模型，精准处理观测数据，实时观测数据实时处理数据，判断基线长以及地表监测点的三维坐标是否有变形，坐标位移量是否在预警数据精度范围；
步骤3：预警；
接收处理后的倾斜测量子系统和GPS监测子系统传输的倾斜测量数据和GPS数据，分析得到预警数据，并根据设置的预警值，发出预警信号；
如果不在预警数据范围内，则对当前倾斜测量子系统和GPS监测子系统观测值进行最小二乘估计，然后重复迭代，全天候、实时的连续进行观测。
6.根据权利要求5所述的基于倾斜仪和GPS的火山监测与预警方法，其特征在于，步骤1中倾斜角度θ与输出的电量读数F关系为：
θ＝a+b*F+c*F+d*F；
式中，θ为被测结构物的倾斜角度，单位为度(°)；F为倾斜仪的实时电量测量值，单位为F；a﹑b﹑c﹑d为倾斜仪的标定系数。
7.根据权利要求5所述的基于倾斜仪和GPS的火山监测与预警方法，其特征在于，步骤2中，GPS基线向量的估值为：
<mrow> <msub> <mi>B</mi> <mi>i</mi> </msub> <mo>=</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>b</mi> <mrow> <mi>i</mi> <mo>,</mo> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>b</mi> <mrow> <mi>i</mi> <mo>,</mo> <mn>2</mn> </mrow> </msub> </mtd> <mtd> <mn>......</mn> </mtd> <mtd> <msub> <mi>b</mi> <mrow> <mi>i</mi> <mo>,</mo> <msub> <mi>m</mi> <mi>i</mi> </msub> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>;</mo> </mrow>
mi-1是函数独立的基线向量，bi,k为第K条基线向量；
基线向量估值的验后方差-协方差阵为：
<mrow> <msub> <mi>D</mi> <mrow> <mi>B</mi> <mi>i</mi> </mrow> </msub> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <msub> <mi>d</mi> <mrow> <msub> <mi>b</mi> <mrow> <mi>i</mi> <mo>,</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>b</mi> <mrow> <mi>i</mi> <mo>,</mo> <mn>1</mn> </mrow> </msub> </mrow> </msub> </mtd> <mtd> <msub> <mi>d</mi> <mrow> <msub> <mi>b</mi> <mrow> <mi>i</mi> <mo>,</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>b</mi> <mrow> <mi>i</mi> <mo>,</mo> <mn>2</mn> </mrow> </msub> </mrow> </msub> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msub> <mi>d</mi> <mrow> <msub> <mi>b</mi> <mrow> <mi>i</mi> <mo>,</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>b</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>d</mi> <mrow> <msub> <mi>b</mi> <mrow> <mi>i</mi> <mo>,</mo> <mn>2</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>b</mi> <mrow> <mi>i</mi> <mo>,</mo> <mn>1</mn> </mrow> </msub> </mrow> </msub> </mtd> <mtd> <msub> <mi>d</mi> <mrow> <mi>b</mi> <mn>2</mn> <mo>,</mo> <msub> <mi>b</mi> <mrow> <mi>i</mi> <mo>,</mo> <mn>2</mn> </mrow> </msub> </mrow> </msub> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msub> <mi>d</mi> <mrow> <mi>b</mi> <mn>2</mn> <mo>,</mo> <msub> <mi>b</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <mtable> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> </mtable> </mtd> <mtd> <mtable> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> </mtable> </mtd> <mtd> <mtable> <mtr> <mtd> <mrow></mrow> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> </mtr> </mtable> </mtd> <mtd> <mtable> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> </mtable> </mtd> </mtr> <mtr> <mtd> <msub> <mi>d</mi> <mrow> <msub> <mi>b</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>b</mi> <mrow> <mi>i</mi> <mo>,</mo> <mn>1</mn> </mrow> </msub> </mrow> </msub> </mtd> <mtd> <msub> <mi>d</mi> <mrow> <msub> <mi>b</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>b</mi> <mrow> <mi>i</mi> <mo>,</mo> <mn>1</mn> </mrow> </msub> </mrow> </msub> </mtd> <mtd> <mo>...</mo> </mtd> <mtd> <msub> <mi>d</mi> <mrow> <msub> <mi>b</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>b</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>m</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow>
是基线向量k、l间的协方差子阵；对GPS接收机获得的星历历元数据进行GPSRTK处理，组设的载波相位双差观测方程为：
<mrow> <mo>&amp;dtri;</mo> <msubsup> <mi>&amp;Delta;&amp;Phi;</mi> <mn>12</mn> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mi>f</mi> <mi>c</mi> </mfrac> <mo>&amp;dtri;</mo> <msubsup> <mi>&amp;Delta;&amp;rho;</mi> <mn>12</mn> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <mo>&amp;dtri;</mo> <msubsup> <mi>&amp;Delta;N</mi> <mn>12</mn> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msubsup> <mo>+</mo> <mi>&amp;epsiv;</mi> <mo>;</mo> </mrow>
式中是两测站两卫星之间的双差载波相位，是两测站两卫星间双差整周模糊度参数，是两测站两卫星间双差星站距离，ε是测量噪声，c为光速，f是信号频率；
双差观测方程的误差方程模型为：
<mrow> <msup> <mi>v</mi> <mi>k</mi> </msup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <mfrac> <mi>f</mi> <mi>c</mi> </mfrac> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>&amp;Delta;l</mi> <mn>2</mn> <mi>k</mi> </msubsup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <msubsup> <mi>&amp;Delta;m</mi> <mn>2</mn> <mi>k</mi> </msubsup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <msubsup> <mi>&amp;Delta;n</mi> <mn>2</mn> <mi>k</mi> </msubsup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msub> <mi>&amp;delta;x</mi> <mn>2</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&amp;delta;y</mi> <mn>2</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&amp;delta;z</mi> <mn>2</mn> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <msup> <mi>&amp;Delta;&amp;Delta;N</mi> <mi>k</mi> </msup> <mo>+</mo> <msup> <mi>&amp;Delta;&amp;Delta;l</mi> <mi>k</mi> </msup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
式中vk为观测值误差，δx2、δy2、δz2是监测站k的坐标近似值向量的改正数向量，为系数矩阵，△△lk为监测站K的常数项误差。
CN201710370031.XA 2017-05-23 2017-05-23 一种基于倾斜仪和gps的火山监测与预警系统及方法 Pending CN107240228A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710370031.XA CN107240228A (zh) 2017-05-23 2017-05-23 一种基于倾斜仪和gps的火山监测与预警系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710370031.XA CN107240228A (zh) 2017-05-23 2017-05-23 一种基于倾斜仪和gps的火山监测与预警系统及方法

Publications (1)

Publication Number Publication Date
CN107240228A true CN107240228A (zh) 2017-10-10

Family

ID=59985854

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710370031.XA Pending CN107240228A (zh) 2017-05-23 2017-05-23 一种基于倾斜仪和gps的火山监测与预警系统及方法

Country Status (1)

Country Link
CN (1) CN107240228A (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108955770A (zh) * 2018-07-25 2018-12-07 成都天地量子科技有限公司 一种基于多光谱卫星影像的火山活动监测方法
CN108981659A (zh) * 2018-07-27 2018-12-11 中国十七冶集团有限公司 一种全自动gps位移观测设备的测量方法
CN110230990A (zh) * 2019-06-27 2019-09-13 北京讯腾智慧科技股份有限公司 燃气管网地质形变监测系统和方法
CN111797491A (zh) * 2020-04-30 2020-10-20 中国空间技术研究院 华北平原地壳垂向位移季节性和时空变化分析方法和系统
CN112444187A (zh) * 2019-08-28 2021-03-05 千寻位置网络有限公司 形变监测方法及装置
CN116304541A (zh) * 2023-05-25 2023-06-23 山东大学 一种基于多源地球物理参数的火山爆发分析方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006317257A (ja) * 2005-05-12 2006-11-24 National Research Institute For Earth Science & Disaster Provention 火山活動の自動監視システム
CN106441174A (zh) * 2016-09-09 2017-02-22 桂林电子科技大学 一种高边坡形变监测方法及系统

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006317257A (ja) * 2005-05-12 2006-11-24 National Research Institute For Earth Science & Disaster Provention 火山活動の自動監視システム
CN106441174A (zh) * 2016-09-09 2017-02-22 桂林电子科技大学 一种高边坡形变监测方法及系统

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
R.SCARPA等: "《火山监测与减灾》", 30 September 2001 *
余学祥等: "《GPS测量与数据处理》", 31 October 2013 *
周云等: "《土木工程防灾减灾概论》", 30 November 2005 *
杜玉柱等: "《GNSS测量技术》", 28 February 2013 *
杨清福等: "长白山天池火山监测预警系统建设架构(II)——火山预警平台建设", 《自然灾害学报》 *
杨清福等: "长白山天池火山预警系统建设架构(I)——火山监测网络建设", 《自然灾害学报》 *
潘懋等: "《环境地质学》", 31 October 2003 *
绝望的童子鸡: "《百度文库》", 15 September 2015 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108955770A (zh) * 2018-07-25 2018-12-07 成都天地量子科技有限公司 一种基于多光谱卫星影像的火山活动监测方法
CN108981659A (zh) * 2018-07-27 2018-12-11 中国十七冶集团有限公司 一种全自动gps位移观测设备的测量方法
CN110230990A (zh) * 2019-06-27 2019-09-13 北京讯腾智慧科技股份有限公司 燃气管网地质形变监测系统和方法
CN112444187A (zh) * 2019-08-28 2021-03-05 千寻位置网络有限公司 形变监测方法及装置
CN111797491A (zh) * 2020-04-30 2020-10-20 中国空间技术研究院 华北平原地壳垂向位移季节性和时空变化分析方法和系统
CN111797491B (zh) * 2020-04-30 2024-04-12 中国空间技术研究院 华北平原地壳垂向位移季节性和时空变化分析方法和系统
CN116304541A (zh) * 2023-05-25 2023-06-23 山东大学 一种基于多源地球物理参数的火山爆发分析方法及系统
CN116304541B (zh) * 2023-05-25 2023-09-01 山东大学 一种基于多源地球物理参数的火山爆发分析方法及系统

Similar Documents

Publication Publication Date Title
CN107240228A (zh) 一种基于倾斜仪和gps的火山监测与预警系统及方法
CN106441174B (zh) 一种高边坡形变监测方法及系统
CN106255065B (zh) 智能手机室内外无缝定位系统及其方法
CN103822629B (zh) 基于多方向偏振光导航传感器的定位系统及其定位方法
CN108318904B (zh) 基于北斗定位技术的高压输电线路舞动监测系统及方法
CN109061696A (zh) 一种确定导航卫星轨道和钟差的方法
CN108508469A (zh) 一种基于北斗前置解算的电力铁塔形变监测系统及其监测方法
CN105467415B (zh) 一种基于差分气压高度约束的小型无人机rtk相对定位方法
CN103217177B (zh) 一种电波折射修正方法、装置及系统
CN203908551U (zh) 一种基站天线工参的自动获取装置
CN207601192U (zh) 姿态可感知雷暴云三维电场探空仪
CN105527642B (zh) 一种单星定位装置与方法
US11608170B2 (en) Buoy position monitoring method and buoy position monitoring system
CN103399335A (zh) 一种移动平台测试系统及误差补偿算法
US20180017392A1 (en) Apparatus, system, and method for traffic
CN107064979A (zh) 一种分组式rtk定位方法及系统
CN109507695A (zh) 一种基于gnss卫星信号的铁塔运行状态监控系统
CN104913743A (zh) 基于惯性测量的电力铁塔变形监测方法
CN114879222A (zh) 一种基于自适应随机模型的全球电离层建模方法
CN117194928B (zh) 基于gnss的地理形变监测系统
CN112632473B (zh) 一种融合地基和空基gnss大气可降水量的计算方法
CN107678046A (zh) 一种便携式高精度全球卫星导航系统接收机
CN116679331A (zh) 基于轮询北斗rtk定位与lstm辅助解算的杆塔偏移监测方法及系统
CN205539539U (zh) 一种低空数字探空仪
CN204464473U (zh) 定位天线装置

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20171010