CN113805223A - 一种地震勘探系统及其采集数据的处理方法 - Google Patents
一种地震勘探系统及其采集数据的处理方法 Download PDFInfo
- Publication number
- CN113805223A CN113805223A CN202110936146.7A CN202110936146A CN113805223A CN 113805223 A CN113805223 A CN 113805223A CN 202110936146 A CN202110936146 A CN 202110936146A CN 113805223 A CN113805223 A CN 113805223A
- Authority
- CN
- China
- Prior art keywords
- data
- time
- gps
- gps satellite
- acquisition
- 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
Links
- 238000003672 processing method Methods 0.000 title abstract description 6
- 239000013078 crystal Substances 0.000 claims abstract description 17
- 238000000034 method Methods 0.000 claims abstract description 8
- 238000001914 filtration Methods 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 24
- 239000013598 vector Substances 0.000 claims description 15
- 238000012545 processing Methods 0.000 claims description 11
- 238000005070 sampling Methods 0.000 claims description 7
- 238000007689 inspection Methods 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 4
- 230000010355 oscillation Effects 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 4
- 241001652065 Trigonopeltastes delta Species 0.000 claims description 3
- 239000000969 carrier Substances 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 3
- 238000013016 damping Methods 0.000 claims description 3
- 238000013500 data storage Methods 0.000 claims description 3
- 238000003745 diagnosis Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 230000008439 repair process Effects 0.000 claims description 3
- 230000035945 sensitivity Effects 0.000 claims description 3
- 230000001360 synchronised effect Effects 0.000 claims description 3
- 230000009897 systematic effect Effects 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 abstract description 2
- 239000005433 ionosphere Substances 0.000 abstract description 2
- 239000005436 troposphere Substances 0.000 abstract description 2
- 238000003491 array Methods 0.000 abstract 1
- 238000013480 data collection Methods 0.000 abstract 1
- 230000003313 weakening effect Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003208 petroleum Substances 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
-
- 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
- G01S19/44—Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method
-
- 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
Abstract
本发明具体涉及一种地震勘探系统及其采集数据的处理方法,在时间同步方面,本发明利用低电源电压和高度可靠的TCXO晶体振荡器,对整个地震数据进行时时标记,其准确性仅取决于GPS精度和数据采集系统的时钟周期,且其稳定性不会随采集节点的工作时间和空间的变化而减弱;在定位方面,利用载波相位和伪距观测双差组合的方法削弱卫星钟差、接收机钟差、电离层延迟误差和对流层延迟误差,并利用卡尔曼滤波算法消除随机干扰,提高了定位精度,可用于高密度地震勘探阵列。
Description
技术领域
本发明属于地震勘探领域,具体涉及一种地震勘探系统及其采集数据的处理方法。
背景技术
地震勘探是目前最主要的石油勘探方法,工作时,通常将传感器节点密集地布设在数十公里的山区和其他复杂的地形勘探区域中,震源激发多处地震波,所有传感器节点实时记录地震数据,形成了密集的地震节点观测系统,实现了对地下地质构造的精细探测。采集任务结束时,所有节点进行集中地震数据下载,然后对大量地震数据进行统一处理和解释。而同步能力差的地震数据将在后期地震数据的处理中会引起相当大的错误,一方面,随着采集时间增加,采集站上晶体振荡器的频率漂移也会带来大量的时间误差累积;另一方面,同步误差随着节点数量的增加而增加。此外,定位信息的结果也极大的影响着地震勘探数据的处理,因此,需要研究用于高密度地震阵列测量的多节点时间同步和高精度定位系统。
吉林大学自主研制的GElWSK地震仪通过内置高精度GPS接收模块为采集站提供时间服务,以解决分布式传感器节点之间同步精度差的问题,但GPS设备的长期运行会消耗大量能量。另外,在某些复杂区域中并没有GPS信号,也会给地震勘探工作带来困难。Fairfield Nodal公司推出的Z-land节点地震仪是一种小型、无电缆的便携式存储式地震采集系统,采用GPS接收器和本地时钟系统实现地震数据的时间标记,以实现低能耗情况下的多节点时间同步。但由于工作温度、工作状态等的变化,长时间工作时,仪器的本地时钟会漂移并导致较大的误差。此外,上述地震仪在定位处理上定位信息精度不高。
发明内容
本发明针对现有技术中的不足,提供一种地震勘探系统及其采集数据的处理方法,本发明的技术方案如下:
一种地震勘探系统,包括控制中心、GPS卫星、移动站、基准站和若干采集节点,所述采集节点包括中央控制单元MCU、GPS模块、TCXO晶体振荡器、存储模块、电源、DAC模块和ADC模块;
所述电源负责为整个采集节点供电;
所述中央控制单元MCU负责管理节点的数据采集任务;
所述TCXO晶体振荡器为DAC模块与ADC模块提供时间服务和标准时间基准;
所述DAC模块和ADC模块为数据采集单元,用于采集地震数据;
所述GPS模块协同GPS卫星用于获取定位信息和时间信息;
所述存储模块包括FIFO存储器和SD卡,FIFO存储器用于零时缓存数据,SD卡用于存储最终的节点地震数据;
地震勘探时,将若干采集节点布设于野外各区域,并开启电源进行节点自检;自检合格后GPS模块接收GPS卫星同步信号,开始采集地震数据,中央控制单元MCU将采集到的地震数据存入FIFO存储器,当FIFO存储器中的数据存储达到一定容量后,中央控制单元MCU从FIFO存储器中读取数据并以数据块文件的形式存储于SD卡中;勘探结束后回收采集节点,并将各采集节点SD卡中的数据下载聚合到控制中心的数据库中进行后期数据处理。
进一步地,所述节点自检的具体流程为:采集节点执行节点自检,自检指标包括检波器电阻、灵敏度、阻尼、自然频率、排列噪声和温度;发送包含测试参数的自检消息到控制中心;控制中心接收每个采集节点的自检消息并判断采集节点是否正常;若采集节点正常则自检合格,控制中心将为每个采集节点配置采集参数并启动采集任务,否则给出警示以提醒相关人员进行故障诊断并手动修复。
进一步地,所述GPS卫星同步信号包括GPS定位信号和GPS时间信号,采集节点根据GPS定位信号获取伪距单点定位观测数据,并将包含伪距单点定位观测数据的地震数据存储于SD卡中,在不需要高精度定位的情况下,可直接使用该伪距单点定位观测数据获取定位信息,采集节点根据GPS时间信号获取GPS时间,并将获取到的GPS时间作为SD卡中时间同步文件夹的名称,SD卡中每个地震数据块都有时间戳。
进一步地,所述ADC模块采用型号为ADS1282的32位高分辨率ADC,ADS1282的数据输出引脚DRDY与中央控制单元MCU的中断引脚相连,中央控制单元MCU在DRDY输出信号的下降沿读取数据;ADS1282的同步输入引脚SYNC施加有连续时钟信号且其频率为数据输出速率的整数倍;地震数据就绪后中央控制单元MCU产生中断,ADC将32位数据发送到中央控制单元MCU;当DRDY发生中断时,一个数据点就是输出,中央控制单元MCU统计ADC数据点的数量,并在采样点达到一定数量时拍摄GPS时间和ADC数据点数量的快照用于对齐数据;当发生时间同步事件时,ADC生成的数据点数量以及准确的GPS时间将作为时间同步文件存储于SD卡中。
进一步地,所述TCXO晶体振荡器的震荡频率fo为16.384MHz±0.2PPM,时间T内任意两个采集节点的累积时间误差Δtmax为2TΔfmax/fo,Δfmax为晶体振荡器最大偏差;中央控制单元MCU统计FIFO存储器的写入次数,当写入次数达到触发时间时会通知GPS模块串行数据输出中断和PPS中断开始时间同步,触发时间为1/FΔtmax,F为采集节点的数据采样率。
进一步地,当采集到的地震数据大小达到预设的数据块大小时,关闭当前文件并创建下一个数据文件,同时将ADC数据点的计数和FIFO存储器的写入次数计数置零。
基于所述地震勘探系统的采集数据的处理方法,包括高精度定位信息获取和数据单元解析,具体步骤如下:
步骤1:根据采集到的数据,建立如下单差观测方程:
式中,c为光速,分别为移动站、基准站接收载波关于GPS卫星j的伪距观测值, 分别为移动站、基准站接收载波关于GPS卫星j的载波相位观测值,分别为移动站、基准站的位置与GPS卫星j的距离,δtr、δtb分别为移动站、基准站的接收机钟差,λ1、f1 j分别为接收载波的波长和频率,分别为移动站、基准站关于GPS卫星j的整周模糊度,Δε1P,rb为伪距单差测量噪声,Δε1Φ,rb为载波相位单差测量噪声;将单差观测方程进行伪距、载波相位结合,得进而获得单差模糊度浮点解;
步骤2:根据单差观测方程,构建如下双差观测方程:
式中,表示t时刻GPS卫星i与GPS卫星j的载波相位观测值的星间差分,分别表示t时刻移动站、基准站接收载波关于GPS卫星i的载波相位观测值,分别表示t时刻移动站、基准站接收载波关于GPS卫星j的载波相位观测值,表示t时刻GPS卫星i与GPS卫星j的距离的星间差分,分别表示t时刻移动站、基准站与GPS卫星i的距离,分别表示t时刻移动站、基准站与GPS卫星j的距离,表示t时刻GPS卫星i与GPS卫星j的整周模糊度的星间差分,分别表示t时刻移动站、基准站关于GPS卫星i的整周模糊度,分别表示t时刻移动站、基准站关于GPS卫星j的整周模糊度;
步骤3:采用扩展卡尔曼滤波算法将双差观测方程线性化,扩展卡尔曼滤波算法的模型方程为:
式中,Xk、Xk分别表示系统的当前状态和上一时刻的状态,Φ为状态转移矩阵,B为单位矩阵,wk为系统误差向量且为双差观测方程状态向量的协方差矩阵,Zk表示系统的当前观测值,H为双差观测方程的系数矩阵,ek为观测误差向量且ek~N(0,Rk),Rk为载波相位-伪距测量误差的协方差矩阵;双差观测方程线性化得:
式中,Q为预测误差的协方差矩阵,K为卡尔曼增益,v1=Zk-HTXk,Xp为双差浮点解,Σp为状态向量的协方差矩阵;
步骤4:在获得双差浮点解和状态向量的协方差矩阵后,基于自反馈因子的蚁群算法求解整周模糊度,获取最终的定位信息;
步骤5:解析地震数据单元,并分析获得地震勘探结果;控制中心的服务器根据GPS时间信息和时间同步文件对地震数据进行更新和标记,并自动分析每个采集节点的时钟偏移,以保证采集数据的时间准确性,达到多采集节点时间同步的目的。
本发明的有益效果是:在时间方面,本发明利用在温度范围内频率的稳定性,使用低电源电压和高度可靠的TCXO晶体振荡器,对整个地震数据进行时时标记;在定位方面,利用载波相位和伪距观测双差组合的方法,以削弱卫星钟差、接收机钟差、电离层延迟误差和对流层延迟误差,获得更高精度的定位。
附图说明
图1为本发明中地震勘探系统的框架示意图。
图2为本发明中采集节点的组成模块示意图。
图3为本发明地震勘探系统的数据采集及处理流程图。
具体实施方式
现在结合附图对本发明作进一步详细的说明。
如图1所示,本发明地震勘探系统主要包括控制中心、GPS卫星、移动站、基准站和若干采集节点,进行地震勘探时,移动站、基准站和若干采集节点被布设于野外。其中,采集节点的组成模块如图2所示,主要包括中央控制单元MCU、GPS模块、TCXO晶体振荡器、存储模块、电源、DAC模块和ADC模块。其中,电源负责为整个采集节点供电;中央控制单元MCU负责管理节点的数据采集任务;TCXO晶体振荡器为DAC模块与ADC模块提供时间服务和标准时间基准;DAC模块和ADC模块为数据采集单元,用于采集地震数据;GPS模块协同GPS卫星用于获取定位信息和时间信息;存储模块包括FIFO存储器和SD卡,FIFO存储器用于零时缓存数据,SD卡用于存储最终的节点地震数据。
图3为本发明地震勘探系统的数据采集及处理流程图,地震勘探时,首先将若干采集节点布设于野外各区域,并开启电源进行节点自检,自检指标包括检波器电阻、灵敏度、阻尼、自然频率、排列噪声和温度。然后将包含测试参数的自检消息发送到控制中心,控制中心接收每个采集节点的自检消息并判断采集节点是否正常。若采集节点正常则自检合格,控制中心将为每个采集节点配置采集参数并启动采集任务,采集节点的GPS模块接收GPS卫星同步信号,开始采集地震数据,中央控制单元MCU将采集到的地震数据存入FIFO存储器,当FIFO存储器中的数据存储达到一定容量后,中央控制单元MCU从FIFO存储器中读取数据并以数据块文件的形式存储于SD卡中。若采集节点异常,则给出警示以提醒相关人员进行故障诊断并手动修复。勘探结束后回收采集节点,并将各采集节点SD卡中的数据下载聚合到控制中心的数据库中进行后期数据处理。
GPS卫星同步信号包括GPS定位信号和GPS时间信号,采集节点根据GPS定位信号获取伪距单点定位观测数据,并将包含伪距单点定位观测数据的地震数据存储于SD卡中,在不需要高精度定位的情况下,可直接使用该伪距单点定位观测数据获取定位信息;采集节点根据GPS时间信号获取GPS时间,并将获取到的GPS时间作为SD卡中时间同步文件夹的名称。由于地震数据块的大小是根据特定数据设置的,以一定的采样率采集时,根据文件夹目录和文件名建立数据文件索引,因此SD卡中每个地震数据块都有时间戳。
本实施例采集节点的ADC模块采用型号为ADS1282的32位高分辨率ADC,ADS1282的数据输出引脚DRDY与中央控制单元MCU的中断引脚相连,中央控制单元MCU在DRDY输出信号的下降沿读取数据;ADS1282的同步输入引脚SYNC施加有连续时钟信号且其频率为数据输出速率的整数倍;地震数据就绪后中央控制单元MCU产生中断,ADC将32位数据发送到中央控制单元MCU;当DRDY发生中断时,一个数据点就是输出,中央控制单元MCU统计ADC数据点的数量,并在采样点达到一定数量时拍摄GPS时间和ADC数据点数量的快照用于对齐数据;当发生时间同步事件时,ADC生成的数据点数量以及准确的GPS时间将作为时间同步文件存储于SD卡中。
模数转换同步的精度取决于系统时钟的精度,ADC和DAC的系统时钟主要依靠晶振。在不同的工作状态下,晶体振荡器经过很长一段时间的振荡后会出现频率漂移,这将导致多采集节点之间的时间异步,进而导致采集的数据出现问题。本实施例所用TCXO晶体振荡器的震荡频率fo为16.384MHz±0.2PPM,时间T内任意两个采集节点的累积时间误差Δtmax为2TΔfmax/fo,Δfmax为晶体振荡器最大偏差;中央控制单元MCU统计FIFO存储器的写入次数,当写入次数达到触发时间时会通知GPS模块串行数据输出中断和PPS中断开始时间同步,触发时间为1/FΔtmax,F为采集节点的数据采样率。当采集到的地震数据大小达到预设的数据块大小时,关闭当前文件并创建下一个数据文件,同时将ADC数据点的计数和FIFO存储器的写入次数计数置零。
本发明地震勘探系统的采集数据的处理包括高精度定位信息获取和数据单元解析,具体的处理步骤如下:
步骤1:根据采集到的数据,建立如下单差观测方程:
式中,c为光速,分别为移动站、基准站接收载波关于GPS卫星j的伪距观测值, 分别为移动站、基准站接收载波关于GPS卫星j的载波相位观测值,分别为移动站、基准站的位置与GPS卫星j的距离,δtr、δtb分别为移动站、基准站的接收机钟差,λ1、f1 j分别为接收载波的波长和频率,分别为移动站、基准站关于GPS卫星j的整周模糊度,Δε1P,rb为伪距单差测量噪声,Δε1Φ,rb为载波相位单差测量噪声;将单差观测方程进行伪距、载波相位结合,得进而获得单差模糊度浮点解。
步骤2:根据单差观测方程,构建如下双差观测方程:
式中,表示t时刻GPS卫星i与GPS卫星j的载波相位观测值的星间差分,分别表示t时刻移动站、基准站接收载波关于GPS卫星i的载波相位观测值,分别表示t时刻移动站、基准站接收载波关于GPS卫星j的载波相位观测值,表示t时刻GPS卫星i与GPS卫星j的距离的星间差分,分别表示t时刻移动站、基准站与GPS卫星i的距离,分别表示t时刻移动站、基准站与GPS卫星j的距离,表示t时刻GPS卫星i与GPS卫星j的整周模糊度的星间差分,分别表示t时刻移动站、基准站关于GPS卫星i的整周模糊度,分别表示t时刻移动站、基准站关于GPS卫星j的整周模糊度。
步骤3:采用扩展卡尔曼滤波算法将双差观测方程线性化,扩展卡尔曼滤波算法的模型方程为:
式中,Xk、Xk分别表示系统的当前状态和上一时刻的状态,Φ为状态转移矩阵,B为单位矩阵,wk为系统误差向量且为双差观测方程状态向量的协方差矩阵,Zk表示系统的当前观测值,H为双差观测方程的系数矩阵,ek为观测误差向量且ek~N(0,Rk),Rk为载波相位-伪距测量误差的协方差矩阵;双差观测方程线性化得:
式中,Q为预测误差的协方差矩阵,K为卡尔曼增益,v1=Zk-HTXk,Xp为双差浮点解,Σp为状态向量的协方差矩阵。
步骤4:在获得双差浮点解和状态向量的协方差矩阵后,基于自反馈因子的蚁群算法求解整周模糊度,获取最终定位信息。
步骤5:解析地震数据单元,并分析获得地震勘探结果;控制中心的服务器根据GPS时间信息和时间同步文件对地震数据进行更新和标记,并自动分析每个采集节点的时钟偏移,以保证采集数据的时间准确性,达到多采集节点时间同步的目的。
以上仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施例,凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理前提下的若干改进和润饰,应视为本发明的保护范围。
Claims (7)
1.一种地震勘探系统,包括控制中心、GPS卫星、移动站、基准站和若干采集节点,其特征在于,所述采集节点包括中央控制单元MCU、GPS模块、TCXO晶体振荡器、存储模块、电源、DAC模块和ADC模块;
所述电源负责为整个采集节点供电;
所述中央控制单元MCU负责管理节点的数据采集任务;
所述TCXO晶体振荡器为DAC模块与ADC模块提供时间服务和标准时间基准;
所述DAC模块和ADC模块为数据采集单元,用于采集地震数据;
所述GPS模块协同GPS卫星用于获取定位信息和时间信息;
所述存储模块包括FIFO存储器和SD卡,FIFO存储器用于零时缓存数据,SD卡用于存储最终的节点地震数据;
地震勘探时,将若干采集节点布设于野外各区域,并开启电源进行节点自检;自检合格后GPS模块接收GPS卫星同步信号,开始采集地震数据,中央控制单元MCU将采集到的地震数据存入FIFO存储器,当FIFO存储器中的数据存储达到一定容量后,中央控制单元MCU从FIFO存储器中读取数据并以数据块文件的形式存储于SD卡中;勘探结束后回收采集节点,并将各采集节点SD卡中的数据下载聚合到控制中心的数据库中进行后期数据处理。
2.如权利要求1所述的一种地震勘探系统,其特征在于,所述节点自检的具体流程为:采集节点执行节点自检,自检指标包括检波器电阻、灵敏度、阻尼、自然频率、排列噪声和温度;发送包含测试参数的自检消息到控制中心;控制中心接收每个采集节点的自检消息并判断采集节点是否正常;若采集节点正常则自检合格,控制中心将为每个采集节点配置采集参数并启动采集任务,否则给出警示以提醒相关人员进行故障诊断并手动修复。
3.如权利要求1所述的一种地震勘探系统,其特征在于,所述GPS卫星同步信号包括GPS定位信号和GPS时间信号,采集节点根据GPS定位信号获取伪距单点定位观测数据,并将包含伪距单点定位观测数据的地震数据存储于SD卡中,在不需要高精度定位的情况下,可直接使用该伪距单点定位观测数据获取定位信息,采集节点根据GPS时间信号获取GPS时间,并将获取到的GPS时间作为SD卡中时间同步文件夹的名称,SD卡中每个地震数据块都有时间戳。
4.如权利要求3所述的一种地震勘探系统,其特征在于,所述ADC模块采用型号为ADS1282的32位高分辨率ADC,ADS1282的数据输出引脚DRDY与中央控制单元MCU的中断引脚相连,中央控制单元MCU在DRDY输出信号的下降沿读取数据;ADS1282的同步输入引脚SYNC施加有连续时钟信号且其频率为数据输出速率的整数倍;地震数据就绪后中央控制单元MCU产生中断,ADC将32位数据发送到中央控制单元MCU;当DRDY发生中断时,一个数据点就是输出,中央控制单元MCU统计ADC数据点的数量,并在采样点达到一定数量时拍摄GPS时间和ADC数据点数量的快照用于对齐数据;当发生时间同步事件时,ADC生成的数据点数量以及准确的GPS时间将作为时间同步文件存储于SD卡中。
5.如权利要求1所述的一种地震勘探系统,其特征在于,所述TCXO晶体振荡器的震荡频率fo为16.384MHz±0.2PPM,时间T内任意两个采集节点的累积时间误差Δtmax为2TΔfmax/fo,Δfmax为晶体振荡器最大偏差;中央控制单元MCU统计FIFO存储器的写入次数,当写入次数达到触发时间时会通知GPS模块串行数据输出中断和PPS中断开始时间同步,触发时间为1/FΔtmax,F为采集节点的数据采样率。
6.如权利要求4或5所述的一种地震勘探系统,其特征在于,当采集到的地震数据大小达到预设的数据块大小时,关闭当前文件并创建下一个数据文件,同时将ADC数据点的计数和FIFO存储器的写入次数计数置零。
7.基于权利要求1所述的一种地震勘探系统的采集数据的处理方法,其特征在于,包括高精度定位信息获取和数据单元解析,具体步骤如下:
步骤1:根据采集到的数据,建立如下单差观测方程:
式中,c为光速,分别为移动站、基准站接收载波关于GPS卫星j的伪距观测值, 分别为移动站、基准站接收载波关于GPS卫星j的载波相位观测值,分别为移动站、基准站的位置与GPS卫星j的距离,δtr、δtb分别为移动站、基准站的接收机钟差,λ1、f1 j分别为接收载波的波长和频率,分别为移动站、基准站关于GPS卫星j的整周模糊度,Δε1P,rb为伪距单差测量噪声,Δε1Φ,rb为载波相位单差测量噪声;将单差观测方程进行伪距、载波相位结合,得进而获得单差模糊度浮点解;
步骤2:根据单差观测方程,构建如下双差观测方程:
式中,表示t时刻GPS卫星i与GPS卫星j的载波相位观测值的星间差分,分别表示t时刻移动站、基准站接收载波关于GPS卫星i的载波相位观测值,分别表示t时刻移动站、基准站接收载波关于GPS卫星j的载波相位观测值,表示t时刻GPS卫星i与GPS卫星j的距离的星间差分,分别表示t时刻移动站、基准站与GPS卫星i的距离,分别表示t时刻移动站、基准站与GPS卫星j的距离,表示t时刻GPS卫星i与GPS卫星j的整周模糊度的星间差分,分别表示t时刻移动站、基准站关于GPS卫星i的整周模糊度,分别表示t时刻移动站、基准站关于GPS卫星j的整周模糊度;
步骤3:采用扩展卡尔曼滤波算法将双差观测方程线性化,扩展卡尔曼滤波算法的模型方程为:
式中,Xk、Xk分别表示系统的当前状态和上一时刻的状态,Φ为状态转移矩阵,B为单位矩阵,wk为系统误差向量且 为双差观测方程状态向量的协方差矩阵,Zk表示系统的当前观测值,H为双差观测方程的系数矩阵,ek为观测误差向量且ek~N(0,Rk),Rk为载波相位-伪距测量误差的协方差矩阵;双差观测方程线性化得:
式中,Q为预测误差的协方差矩阵,K为卡尔曼增益,v1=Zk-HTXk,Xp为双差浮点解,Σp为状态向量的协方差矩阵;
步骤4:在获得双差浮点解和状态向量的协方差矩阵后,基于自反馈因子的蚁群算法求解整周模糊度,获取最终定位信息;
步骤5:解析地震数据单元,并分析获得地震勘探结果;控制中心的服务器根据GPS时间信息和时间同步文件对地震数据进行更新和标记,并自动分析每个采集节点的时钟偏移,以保证采集数据的时间准确性,达到多采集节点时间同步的目的。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110936146.7A CN113805223A (zh) | 2021-08-16 | 2021-08-16 | 一种地震勘探系统及其采集数据的处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110936146.7A CN113805223A (zh) | 2021-08-16 | 2021-08-16 | 一种地震勘探系统及其采集数据的处理方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113805223A true CN113805223A (zh) | 2021-12-17 |
Family
ID=78943045
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110936146.7A Pending CN113805223A (zh) | 2021-08-16 | 2021-08-16 | 一种地震勘探系统及其采集数据的处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113805223A (zh) |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1584626A (zh) * | 2004-06-08 | 2005-02-23 | 河海大学 | Gps和伪卫星组合定位方法 |
CN103529482A (zh) * | 2013-10-25 | 2014-01-22 | 中国人民解放军国防科学技术大学 | 一种高精度确定载体动态加速度的方法 |
CN104749636A (zh) * | 2015-04-09 | 2015-07-01 | 北京中矿大地地球探测工程技术有限公司 | 分布式槽波地震勘探系统及其勘探方法 |
CN105158793A (zh) * | 2015-07-15 | 2015-12-16 | 安徽吉思勘仪器科技有限公司 | 基于高速现场数据收集的无缆地震数据采集系统 |
CN105446848A (zh) * | 2014-06-13 | 2016-03-30 | 北京安兔兔科技有限公司 | 电子设备的数据处理性能的测试方法及装置 |
CN108317949A (zh) * | 2018-02-07 | 2018-07-24 | 桂林电子科技大学 | 一种rtk高精度差分定位形变监测系统及方法 |
CN108802805A (zh) * | 2018-06-07 | 2018-11-13 | 武汉大学 | 六自由度宽频带一体化地震监测系统 |
CN109358487A (zh) * | 2018-10-10 | 2019-02-19 | 武汉大学 | 一种基于gnss精密授时的伪卫星系统及方法 |
CN109557577A (zh) * | 2019-01-16 | 2019-04-02 | 合肥国为电子有限公司 | 适用于节点地震仪的时钟产生装置及方法 |
CN109917356A (zh) * | 2019-03-13 | 2019-06-21 | 武汉际上导航科技有限公司 | 一种机载激光扫描系统误差标定方法 |
CN209296941U (zh) * | 2019-03-01 | 2019-08-23 | 吉林大学 | 一种基于gps触发的航空物探数据采集系统 |
CN110554430A (zh) * | 2019-09-10 | 2019-12-10 | 山东知微智成电子科技有限公司 | 地震勘探数据无线采集节点的ad转换与数据存储系统及方法 |
CN111239787A (zh) * | 2020-02-28 | 2020-06-05 | 同济大学 | 一种集群自主协同中的gnss动态卡尔曼滤波方法 |
CN111751853A (zh) * | 2020-06-20 | 2020-10-09 | 北京华龙通科技有限公司 | 一种gnss双频载波相位整周模糊度解算方法 |
CN112987603A (zh) * | 2019-12-02 | 2021-06-18 | 中国石油化工集团有限公司 | 一种基于gprs的节点地震仪器远程监控系统 |
-
2021
- 2021-08-16 CN CN202110936146.7A patent/CN113805223A/zh active Pending
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1584626A (zh) * | 2004-06-08 | 2005-02-23 | 河海大学 | Gps和伪卫星组合定位方法 |
CN103529482A (zh) * | 2013-10-25 | 2014-01-22 | 中国人民解放军国防科学技术大学 | 一种高精度确定载体动态加速度的方法 |
CN105446848A (zh) * | 2014-06-13 | 2016-03-30 | 北京安兔兔科技有限公司 | 电子设备的数据处理性能的测试方法及装置 |
CN104749636A (zh) * | 2015-04-09 | 2015-07-01 | 北京中矿大地地球探测工程技术有限公司 | 分布式槽波地震勘探系统及其勘探方法 |
CN105158793A (zh) * | 2015-07-15 | 2015-12-16 | 安徽吉思勘仪器科技有限公司 | 基于高速现场数据收集的无缆地震数据采集系统 |
CN108317949A (zh) * | 2018-02-07 | 2018-07-24 | 桂林电子科技大学 | 一种rtk高精度差分定位形变监测系统及方法 |
CN108802805A (zh) * | 2018-06-07 | 2018-11-13 | 武汉大学 | 六自由度宽频带一体化地震监测系统 |
CN109358487A (zh) * | 2018-10-10 | 2019-02-19 | 武汉大学 | 一种基于gnss精密授时的伪卫星系统及方法 |
CN109557577A (zh) * | 2019-01-16 | 2019-04-02 | 合肥国为电子有限公司 | 适用于节点地震仪的时钟产生装置及方法 |
CN209296941U (zh) * | 2019-03-01 | 2019-08-23 | 吉林大学 | 一种基于gps触发的航空物探数据采集系统 |
CN109917356A (zh) * | 2019-03-13 | 2019-06-21 | 武汉际上导航科技有限公司 | 一种机载激光扫描系统误差标定方法 |
CN110554430A (zh) * | 2019-09-10 | 2019-12-10 | 山东知微智成电子科技有限公司 | 地震勘探数据无线采集节点的ad转换与数据存储系统及方法 |
CN112987603A (zh) * | 2019-12-02 | 2021-06-18 | 中国石油化工集团有限公司 | 一种基于gprs的节点地震仪器远程监控系统 |
CN111239787A (zh) * | 2020-02-28 | 2020-06-05 | 同济大学 | 一种集群自主协同中的gnss动态卡尔曼滤波方法 |
CN111751853A (zh) * | 2020-06-20 | 2020-10-09 | 北京华龙通科技有限公司 | 一种gnss双频载波相位整周模糊度解算方法 |
Non-Patent Citations (5)
Title |
---|
刘新华: "基于卫星导航RTK的精密测量技术研究", 杭州电子科技大学 * |
刘新华;尚俊娜;孙建强;张东丽;施浒立: "一种改进蚁群算法的DGPS整周模糊度求解", 航天控制, vol. 39, no. 001 * |
宋玉龙;廉保旺;唐成凯;: "基于多频模糊度解算的BDS/GPS中基线AEKF RTK算法", 系统工程与电子技术, no. 08 * |
尚俊娜;王民顿;刘新华;王奕腾: "改进的PSO与ACO混合搜索GNSS整周模糊度算法", 中国惯性技术学报, no. 003 * |
李慧;赵琳;李亮;丁继成;: "一种双频非组合实时精密定位方法", 哈尔滨工业大学学报, no. 11 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9405012B2 (en) | Advanced global navigation satellite systems (GNSS) positioning using precise satellite information | |
Lebach et al. | Measurement of the solar gravitational deflection of radio waves using very-long-baseline interferometry | |
RU2450255C2 (ru) | Устройство и способ сбора сейсмических данных, система сейсмических наблюдений | |
CN102608625B (zh) | 基于惯性辅助定位接收机的实时形变监测预警系统及方法 | |
CN105116422A (zh) | 一种星载高动态gnss接收机及其导航方法 | |
CN102494699B (zh) | 捷联式航空重力仪测量参数置信度评估方法 | |
CN102116867A (zh) | 一种在动态环境下探测并修复gps载波相位周跳的方法 | |
Häberling et al. | Assessment of high-rate GPS using a single-axis shake table | |
CN106772493A (zh) | 基于北斗差分定位的无人机航向测算系统及其测算方法 | |
CN108958018A (zh) | 一种卫星授时方法和装置、计算机可读存储介质 | |
CN1967282B (zh) | 全站仪与gps单频实时动态组合测量方法及其系统 | |
CN201266237Y (zh) | 高精度gps差分系统 | |
Li et al. | Precise orbit determination for GRACE with zero-difference kinematic method | |
CN109443337B (zh) | 一种基于金刚石内nv色心的定位导航系统与方法 | |
CN102736118A (zh) | 一种用于全球重力场测量的综合型卫星系统 | |
CN102538826A (zh) | 基于地图匹配数据与gps的自适应陀螺误差校正方法 | |
CN113805223A (zh) | 一种地震勘探系统及其采集数据的处理方法 | |
Shuxin et al. | A study of differential GPS positioning accuracy | |
Lipp et al. | Cycle-slip detection and repair in integrated navigation systems | |
Lachapelle et al. | High precision C/A code technology for rapid static DGPS surveys | |
CN106680845B (zh) | 一种卫星轨道综合定权方法 | |
CN203519837U (zh) | 一种导航计算机及导航系统 | |
Glazer | GPS receiver operation | |
Kosarev et al. | The method of cycle-slip detection and repair GNSS meaturements by using receiver with high stability frequency oscillator | |
Li et al. | Monitoring and research on the displacement and deformation of transmission line towers based on beidou monitoring data |
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 |