CN101826873A - 数据获取设备和方法 - Google Patents
数据获取设备和方法 Download PDFInfo
- Publication number
- CN101826873A CN101826873A CN201010112930A CN201010112930A CN101826873A CN 101826873 A CN101826873 A CN 101826873A CN 201010112930 A CN201010112930 A CN 201010112930A CN 201010112930 A CN201010112930 A CN 201010112930A CN 101826873 A CN101826873 A CN 101826873A
- Authority
- CN
- China
- Prior art keywords
- data
- interpolation
- sampling
- value
- frequency
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 16
- 238000005259 measurement Methods 0.000 claims abstract description 16
- 238000009434 installation Methods 0.000 claims abstract description 3
- 238000005070 sampling Methods 0.000 claims description 78
- 238000004364 calculation method Methods 0.000 claims description 7
- 241001269238 Data Species 0.000 claims description 5
- 238000013481 data capture Methods 0.000 claims description 3
- 230000006870 function Effects 0.000 description 22
- 238000012952 Resampling Methods 0.000 description 10
- 101150071403 INP1 gene Proteins 0.000 description 7
- 101150016601 INP2 gene Proteins 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 230000014509 gene expression Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000005265 energy consumption Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000008676 import Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005039 memory span Effects 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012546 transfer 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/32—Transforming one recording into another or one representation into another
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2200/00—Details of seismic or acoustic prospecting or detecting in general
- G01V2200/10—Miscellaneous details
- G01V2200/12—Clock synchronization-related issues
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)
- Compression, Expansion, Code Conversion, And Decoders (AREA)
- Complex Calculations (AREA)
Abstract
本发明涉及数据获取设备和方法。根据本发明,提供了:第一装置,在由本地时钟(2)提供的不完整采样频率(FE)下提供一系列第一数字采样数据(X),所述第一装置包括sigma-delta调制装置(3),用于提供所述一系列第一数字sigma-delta调制数据;第二校准装置(50),用于根据参考时钟(6)测量本地时钟(2)的时钟误差;第三装置,用于校正至少基于测量到的频率误差的第一数据,第三装置至少包括:内插装置(4),至少基于所述第一系列的内插(X)产生第二数字数据(Y)并补偿测量到的频率误差(FD);以及抽选滤波器(7),用于基于所述第二数字数据(Y)产生第三数字数据(A)。
Description
技术领域
本发明涉及数据获取设备和方法。本发明特别涉及用于地震数据获取的设备和方法,该设备用于形成地震传感器网络的节点。更详细地,本发明涉及地震数据获取设备,该设备由数字信号处理装置构成,该数字信号处理装置连接至均产生一系列数字地震采样或地位移测量的数字地震传感器或与模/数转换器相关联的模拟地震传感器。
背景技术
所有的节点必须同步从而所有节点可以具有相同的采样时间。
本领域普通技术人员将意识到在数字获取网络内实现的多个数据获取设备和方法。
在数据获取网络是由地震传感器组件形成的情况下,这些装置通常由通过线或电缆以组的形式链接在一起的传感器构成,作为网络的“节点”。这些节点通常在“集中器”周围以组布置,并且一组节点与它的集中器形成子网,集中器和节点之间的链路也可以是有线链路。集中器局部地管理节点的子网,提供电源并且集中数据。集中器连接至本地计算机网络,类似于有线网络,驱动组件并记录源于子网的数据的中央单元连接至该网络。这些现有的解决方案对于本领域普通技术人员是公知的,因此不会在下面详细描述。
本领域普通技术人员已经提出了无线数据获取网络,其具有避免使用多个电缆的明显优点。在地震应用的情况下,这些网络也依照上面所述根据组织由节点、集中器和中央单元形成,网络的这些不同元素之间的通信通过RF以射频模式执行。
传统技术是向所有节点仅发送一个采样参考时钟。由于控制本地振动器的锁相环(PLL)系统,每个节点提取和保持该采样时钟。由PLL控制的VCXO(压控晶体振荡器)的使用具有高成本和高能耗,并且涉及高的硬件需求用于发送时钟。
在文档GB2428799中已知另一个数据获取设备。该数据获取设备包括模/数转换器,该模/数转换器在由传感器获取的本地时钟数据提供的不完整(imperfect)的频率处采样;时间戳模块,该时间戳模块通过根据通用时间确定采样数据来测量本地时钟的频率误差,并且包括具有参考时钟的基于卫星的定位系统,该时间戳模块能够将本地时钟校准到基于卫星的定位系统的参考时钟上;当已经执行了本地时钟的校准时,用于关闭基于卫星的定位系统的装置;用于校准采样和定时的数据的重新采样模块。该重新采样模块基于一系列采样的和定时的数据来生成一系列重新采样和重新定时的数据,这就是说,该系列数据的相位误差的采样频率误差均可以被补偿。该重新采样模块包括第一存储器和第二存储器,第一存储器存储在对应于本地时钟的频率的频率处采样和定时的系列数据,第二存储器存储根据通用时间确定的采样数据。重新采样模块包括重新采样器,在该重新采样器的一侧连接包含系列采样和定时数据的第一存储器,在重新采样器的另一侧连接内插滤波器,该内插滤波器本身连接至包含根据通用时间确定的采样数据的第二存储器,并且内插滤波器还连接至参考滤波器,内插滤波器通过参考滤波器动态地计算要被计算的该系列采样的每个采样,由此该内插滤波器产生校正内插因子,从而可以在重新采样器中产生系列重新采样的数据。
所述公知设备在实际制造中是困难并昂贵的。并且还存在需要高计算负载的确定。该内插滤波器具有高阶,即,具有高数目的因子,并且必须对每个采样重新计算。这涉及存储器容量、连续采样之间流逝的时间允许的短时间内的大量计算,以及由此高能耗的严格限制。
发明内容
由此,本发明的目的在于提供一种设备和方法以减轻本领域的上述不便。
本发明的另一个目的在于提供一种设备和方法,能够获得正确的获取数据同时最小化能量消耗、计算负载和存储器负载。
本发明框架中通过数据获取设备可以实现这些目的,所述数据获取设备包括:
第一装置,在第一本地时钟(2)提供的不完整采样频率(FE)下提供一系列第一数字采样定时数据(X),
第二校准装置(50),用于通过第二参考时钟(6)测量所述第一本地时钟(2)的频率误差,
第三装置,用于至少基于测量的频率误差来校准第一数据,
其特征在于,
第一装置(3)包括sigma-delta调制装置(3),设置为产生在由所述第一本地时钟(2)提供的所述第一不完整采样频率(FE)下采样的所述一系列第一数字sigma-delta调制数据(X),
第三装置至少包括内插装置(4),至少基于所述第一系列数据采样(X)的内插来产生第二数字数据(Y)并且补偿所述测量到的频率误差(FD),并且至少包括抽选滤波器(7),用于基于所述第二数字数据(Y)产生第三数字数据(A)。
根据本发明的实施例,第一装置至少包括产生模拟测量数据的模拟传感器(102),并且sigma-delta调制装置(3)是sigma-delta模数转换器(3)的形式,用于将传感器(102)的所述模拟测量数据转换为所述系列第一数字sigma-delta调制数据(X)。
根据本发明的实施例,所述第一装置至少包括数字传感器,具有所述sigma-delta调制装置以产生所述一系列第一数字sigma-delta调制数据(X)。
根据本发明的实施例,该设备包括控制装置(42),设置为至少基于所述测量到的频率误差(FD)设置所述内插。
根据本发明的实施例,内插装置(4)具有内插函数(F),所述内插函数(F)具有预设的固定阶以及至少具有参数(C0、C1、C2),所述第三控制装置(42)包括至少基于所述测量到的频率误差(FD)设置内插函数(F)的所述参数(C0、C1、C2)的装置。
根据本发明的实施例,内插函数(F)的阶小于6。
根据本发明的实施例,内插函数(F)是线性的,并且控制装置(42)根据当前第一采样Xn的以下内插公式设置内插函数(F):
Yn=C0Xn+C1Xn-1
其中C0=Pn
C1=1-Pn
其中Xn-1和Xn分别是连续的第一采样(X),
Yn是对于Xn-1和Xn内插的第三数据(Y),
其中控制装置(42)根据以下公式计算内插因子(P)的值Pn:
Pn=Pn-1+FD
其中,Pn-1和Pn分别是内插因子P的连续值,并且FD是频率漂移。
根据本发明的实施例,内插函数(F)具有等于2的阶,并且控制装置(42)根据当前第一采样Xn的以下内插公式来设置内插函数(F):
Yn=C0Xn+C1Xn-1+C2Xn-2
其中
其中Xn-2、Xn-1和Xn分别是连续的第一采样(X),
Yn是对于Xn-2、Xn-1和Xn的内插后的第三数据(Y),
其中第三控制装置(42)根据以下公式计算内插因子(P)的值Pn:
Pn=Pn-1+FD
其中Pn-1和Pn分别是内插因子P的连续值,并且FD是频率漂移。
根据本发明的实施例,控制装置(42)包括根据方程Pn=Pn-1+FD检查内插因子(P)的值Pn是否大于或等于0且小于或等于1的装置(43、44),并且只要内插因子大于等于0且小于等于1,控制装置(42)根据具有内插因子(P)的值Pn的Xn的内插方程设置内插函数(F)。
根据本发明的实施例,控制装置(42)包括根据方程Pn=Pn-1+FD检查内插因子(P)的值Pn对于当前第一采样Xn是否大于1的装置(44),并且如果根据方程Pn=Pn-1+FD对于当前第一采样Xn计算的内插因子(P)的值Pn大于1,那么控制装置(42)不通过Pn=Pn-1+FD将对于当前第一采样Xn根据内插公式计算的对应的内插后的值Yn作为第三数据(Y),而是通过根据Pn=Pn-1+FD-1计算的内插因子的值Pn,产生对于当前第一采样Xn根据内插公式计算的值Yn作为与当前第一采样Xn相关的第三数据。
根据本发明的实施例,控制装置(42)包括预测装置(43),根据方程Pn+1=Pn+FD检查对于当前第一采样Xn内插因子(P)的下一个值Pn+1是否小于0,并且如果根据方程Pn=Pn-1+FD对于当前第一采样Xn计算的内插因子(P)的值Pn小于0,那么控制装置(42)通过Pn=Pn-1+FD对于当前第一采样Xn使用根据内插公式计算的对应的内插后的值Yn作为第三数据(Y),并且通过根据Pn=Pn-1+FD+1计算的内插因子的值Pn,产生对于当前第一采样Xn根据内插公式计算的另一个值作为第三数据(Y)。
根据本发明的实施例,所述系列第一数字sigma-delta调制数据(X)具有至少为1并且小于等于4的第一预设比特精度(BR1),第二数据(Y)具有所述第一采样频率(FE)以及高于第一比特精度(BR1)的第二预设比特精度(BR2),并且第三数据(A)具有低于第一采样频率(FE)的第三预设采样频率(F3)以及高于第一比特精度(BR1)的第三预设比特精度(BR3)。
根据本发明的实施例,预设第一本地时钟频率FE的值,使得根据过采样率OR有至少一个传感器提供的信号的过采样高于1,其中过采样率OR被定义为等于:
OR=FE/2FU
其中FU是由所述传感器提供的所述信号的最高有用频率。
在本发明框架中通过由计算装置执行的数据获取方法可以实现这些目的,该数据获取方法包括以下步骤:
在第一本地时钟(2)提供的不完整采样频率(FE)下提供一系列第一数字采样定时数据(X),
通过根据第二参考时钟(6)校准第一本地时钟(100)来根据第二参考时钟(6)测量第一本地时钟(100)的频率误差,
至少基于测量到的频率误差来校正第一数据,
其特征在于,
由sigma-delta调制(3)产生在由所述第一本地时钟(2)提供的所述第一不完整采样频率(FE)下采样的所述系列第一数字sigma-delta调制数据(X),
为了补偿测试到的频率误差(FD),对于所述一系列第一数字sigma-delta调制数据(X)进行内插以产生第二数字数据(Y),
对所述第二数字数据(Y)进行抽选以产生第三数字数据(A)。
附图说明
结合附图读取下面的以不受限的方式给出的具体实施方式时,本发明的其他特征、目的和优点将变得显而易见,其中
图1表示根据本发明的数据获取设备的实施例的框图;
图2表示根据本发明的数据获取设备中可以使用的sigma-delta转换器的实例的框图;
图3表示根据本发明的内插的数据的计时图的实例;
图4表示在第一种情况下根据本发明的内插过程中数据流水线(pipeline)的计时图;
图5表示在第二种情况下根据本发明的内插过程中数据流水线的计时图;
图6表示根据本发明具有多个数据获取设备的系统的框图;
图7表示根据本发明的数据获取设备的实施例。
具体实施方式
下面在一个或多个模拟传感器102产生模拟测量数据的情况下描述本发明。模拟传感器的例子是地震检波器、转换传感器、速度传感器、加速计、压力传感器。当然,本发明还可以应用至产生数字测量数据的一个或多个传感器102。
在图1中,设备1专门用于获取由这是一个传感器102测量的地震模拟数据,并且将其发送至第一模拟数据输入10。设备1包括第一本地时钟2,例如由TCXO振荡器(温度补偿晶体振荡器)形成;以及模数转换器3,用于在第一模拟数据输入10上表示的模拟数据E转换为在其转换器输出11上的第一数字数据采样X。第一数字数据采样X具有第一比特精度BR1。转换器3具有连接至第一本地时钟2的第二频率输入12。第一频率输入12设置由转换器3生成的第一数字数据采样X的频率。第一数据采样X是比特串的形式。
转换器3预设输入12上的本地时钟频率FE的值,引起了模拟数据输入10的过采样。将过采样速率OR定位为等于:
OR=FE/2.FU
其中FU是输入10的信号的最高有用频率。采样频率FE高于2FU,这意味着
OR>1
例如
FE>128.FU
OR>64
可以在本地时钟2和第二频率输入12之间插入频率适配器20,从而将本地时钟2提供的信号的频率适配到第二频率输入12上的期望频率FE。例如,频率适配器20具有分频器,从而在频率FE低于时钟2的频率Ftcxo的情况下,将时钟2的频率Ftcxo划分为第二频率输入12上的频率FE。在一个例子中,本地时钟2产生具有若干MHz的频率Ftcxo的信号,例如8MHz,并且FE=256kHz。
提供过采样的模数转换器3包括sigma-delta调制器。在这种情况下是模拟sigma-delta转换器3。Sigma-delta转换器的操作和结构对于本领域普通技术人员来说是公知的,因此在此不再详述。过采样转换器在它的输出OUT11上产生具有1比特精度BR1的信号。
作为图2上的不受限的例子,sigma-delta转换器3主要具有积分器INT,积分器INT的输出连接至比较器COMP第一输入INP1,并且积分器INT的输入连接至具有加输入和减输入的减法器或加法器SUM的输出OUT3,加法器或减法器SUM是所谓的sigma元素,其将两个加和减输入相加。比较器COMP用于将具有预设的固定值的第一输入INP1和第二输入INP2之间的差值转换为一比特输出OUT,其就是所谓的delta输出,即,当第一输入INP1高于第二输出INP2时,在比较器输出OUT上提供比特1,并且当第一输入INP1低于第二输出INP2时,在比较器输出OUT上提供比特0。通过频率输入12施加的频率来设置比较器输出OUT的每个比特专用的周期时隙。比较器输出OUT通过数模转换器CNA被反馈回减法器SUM的减输入。转换器CNA在它的连接至减输入的输出OUT2上产生用于比特1的+Vref信号和用于比特0的-Vref信号,其中Vref是不同于零的固定值。减法器SUM的加输入连接至接收模拟信号V(t)的模拟输入10。结果,加法器SUM在它的输出OUT3上提供模拟输入V(t)和转换器CNA之间的差值,其为OUT3=V(t)-OUT2。当在比较器COMP处,输入INP1的信号高于输入INP2的信号,比较器输出OUT变成比特1,并且将使得数模转换器CNA在输出OUT2上输出+Vref信号。由于连接至数模转换器输出OUT2的加法器输入上的减号,加法器SUM将产生信号OUT3,使得通过积分器INT的输入INP1降低直到输入INP1小于输入INP2。当V(t)>OUT2时,积分器INT将确定积分以迫使OUT2更接近于V(t)。当V(t)<OUT2时,将通过向积分器INT输入应用负信号来降低OUT2。由此,模拟输入信号V(t)越高,越多的比特1将出现在转换器输出OUT,11。输出OUT,11的平均值表示输入信号V(t)的平均值。
在数字传感器102产生数字测量数据的情况下,传感器在输出11上产生所述数字数据X,该数字数据X在传感器102中是sigma-delta调制的。结果,数字传感器102的外壳中集成了产生所述数字数据X的sigma-delta调制器。数字传感器102的例子是数字加速计。
当然,输出OUT,11可以具有多于一个比特,并且例如是一到四比特。
下面描述模数sigma-delta转换器3产生第一个1比特数据采样X的情况下的本发明。
Sigma-delta转换器输出11连接至内插装置4的数据输入40,从而将第一个1比特数据采样X发送至数据输入40。内插装置4例如是内插滤波器的形式。内插装置4包含在内插输出41上产生第二数据Y的装置,第二数据Y具有第一采样频率FE以及高于第一比特精度BR1的第二预设比特精度BR2。通过数据输入40的第一个1比特数据采样X的内插功能计算在内插输出41上的第二数据Y。Xn-1,Xn代表在各个时间Tn和Tn-1上的X的连续1比特数据采样,其中n-1和n代表以速率FE进行的数据采样的连续数目。Yn-1,Yn代表Y的连续数据采样,其中n-1和n代表数据采样的连续数目。结果,
Tn=Tn-1+1/FE
在内插装置4中,计算装置42通过使用本地设置2的频率漂移FD来对每个第二数据Y计算内插函数F。内插函数F具有预设的固定阶(degree)。内插装置4可以实现为FIR滤波器。
控制模块5的计算装置50计算频率漂移或频率误差FD,该控制模块5具有第一频率输入51和第二频率输入51,第一频率输入51连接至本地时钟2用于接收由具有频率FE的本地时钟2提供的第一不完整周期时钟信号TLC,并且第二频率输入52连接至参考时钟输入60,用于接收周期参考时钟信号REF。例如,提供参考时钟的参考装置6提供参考时钟信号REF,该参考装置6例如是基于卫星的定位系统。周期时钟信号TLC由通常被称为在频率FE重复并且能够被检测的周期时间模式的脉冲或前沿或其他信号模式表示。周期参考时钟信号REF由通常被称为在频率FREF重复的周期时间模式的脉冲或前沿或另一信号模式或方波信号表示,频率FREF可以等于或不等于频率FE,并且可以被检测。通过具有例如电缆的传输装置61从参考时钟6向第二频率输入52发送参考时钟信号REF。
参考装置6用作同步源,以提供作为脉冲的时间戳的信号TLC。
下面描述计算装置50的实施例。
参考装置6提供参考频率FREF和参考相位ΦREF,参考相位ΦREF是在世界标准时间(UTC)上的时间戳的脉冲的相位。FREF低于FE,FE低于本地时钟2的频率F2:
FREF<<FE<F2
参考相位ΦREF对应于参考时钟信号REF通过传输装置从参考时钟6到输入52的传播时间。根据节点距离同步源的位置(其间的电缆长度和中继器的数目),ΦREF对于每个节点可以不同。该传播时间是在获取开始之前测量的。
计算装置50基于参考频率FREF测量频率FE。例如,为了测量FE,计数器对于由FREF的周期数目NFREF定义的时间间隔过程中FE的周期数目NFE进行计数。然后如下获得FE:
FE=FREF.NFE/NFREF
然后,基于测量的FE获得频率漂移FD
FD=(FE-FE-SET)/FE-SET
其中FE-SET是对于FE的频率的固定和存储的设置值。
FD表示本地时钟2的测量的频率误差。
计算装置50测量本地设置2的频率漂移FD。FD被表示为时间的TCXO频率比上通常值的百分比数值。频率漂移FD由百万分之几(ppm)的表示。
如下所述提供补偿频率漂移FD的装置。
对于依赖于内插因子P的每个输入采样X设置内插装置4的内插函数F。所述内插因子P具有连续的值Pn-1,Pn。内插因子Pn对于每个采样Xn不同,并且直接依赖于频率漂移FD。
0≤P≤1
在第一实施例中,内插函数F是线性的:
Yn=Xn-1+Pn·(Xn-Xn-1)
其具有下面的转移函数:
其中Z-1表示对于采样Xn的一个时隙TSn。
这样的内插函数F具有两个参数:
对于Xn,C0=Pn
对于Xn-1,C1=Pn-1
所以第一阶(线性)内插函数可以被看做具有两个因子C0,C1的数字滤波器。
对于每个采样Xn,因子Pn增加FD。
所以对于每个采样Xn,内插因子Pn增加对于连续的采样Xn已经计算的频率漂移FD的相同值。
根据以下的公式计算内插因子P的值Pn:
Pn=Pn-1+FD
其中Pn-1是对于在前的采样Xn-1的内插因子P的值。
控制模块5还包括对内插因子P进行初始化的装置54。
装置54测量FE的相位误差ΦFE。在图示的实施例中,通过根据下面的公式,在FE脉冲的上升沿和FREF脉冲的上升沿之间定义的时间中计数本地时钟2的信号Ftcxo的周期Ntcxo的数目来测量FE的相位误差ΦFE:
ΦFE=Ntcxo.1/Ftcxo+ΦFEF
对于初始时间T0,装置54将内插因子P初始化为初始值P0。在获取的开始T0=0时,P被初始化为:
P0=modulus(ΦFE/TE)
其中取模意味着在除以TE之后ΦFE的余数,并且
TE=1/FE
例如,对于ΦFE=3.4μs和FE=1,P0=0.4μs
与时钟信号处理中的通常推荐相反,通常的推荐实际上禁止线性内插,因为由此产生了太多错误,这些错误由于具有高的过采样的转换器3而变得可以忽略了。
图3和7表示对于具有下面的连续采样作为例子的第一个1比特数据采样X,由第一实施例的线性内插装置进行的内插:
X1=0,
X2=0,
X3=1,
X4=1,
X5=0,
X6=0。
图4表示在内插因子变得高于1的情况下,这对应于数据获取设备的本地时钟太快的情况,对应的信号Xn、Xn+1、Yn、TLC和进一步的说明信号FILTLC的对应的计时图。在本地时钟信号TLC中的频率FE上的周期时间模式重复的是前沿。
线性内插产生了下面的连续的第二数据Y。
对于n=2,Y2是基于X1,X2和P2而内插的,由此Y2=0。
对于n=3,Y3是基于X2,X3和P3而内插的,P3=P2+FD。
例如,P3=0.96并且FD=0.03,那么P4=0.99并且Y4=0.01X3+0.99X4。
对于n=4,Y4是基于X3,X4和P4而内插的,P4=P3+FD。
P3=0.96并且FD=0.03,那么P4=0.99并且Y4=0.01X3+0.99X4
假设频率漂移FD计算的值是整数并且是非零的,那么P通过每个进一步采样X增长。
P3和P4已经接近于1(并且小于1)。
提供装置44以检查内插因子P是否大于1。
如果在图3上,由装置44提供的计算得到的内插因子P5大于1,即:
P5=P4+FD>1
那么对于X4和X5在第二数据Y中不使用对应的内插值Y5=X4+P5·(X5-X4)。
对于n=5,P4=0.99并且FD=0.03,那么P5=P4+FD=1.02并且不使用Y5=(1-P5)X4+P5X5=-0.02X4+1.02X5。
如下所述,除了内插因子低于0或高于1的情况,提供装置45以对时钟信息TLC进行滤波,从而产生等于TLC的滤波后的信号FILTLC。
如图4所示,除了在为X4和X5分配的时隙FILTP45之外,滤波后的信号FILTLC等于本地时钟2的周期时钟信号TLC,在时隙FILTP45中取消了TLC的周期时间模式TP45(前沿)。然后,关于滤波后的信号FILTLC中的时隙FILTP45,Y4的时隙持续时钟信号TLC的两个时隙TSn,而Y2、Y3、Y4中每一个的时隙TSn仅持续时钟信号TLC的一个时隙。被内插了错误的内插因子P5(即P5>1)的Y5由在先的内插第三数据(0≤P5≤1)取代。
然后,将Y5延迟一个时钟时隙TSn,并且基于下一个Xn=X6和下一个Xn-1=X5对于它们的下一个时隙重新计算Y5。在这种情况下,从大于1的错误的内插因子P5=P4+FD中减去1,然后变成P5=P4+FD-1。在内插公式Yn=Xn-1+Pn(Xn-Xn-1)中,这意味着Xn由下一个Xn+1取代,并且Xn-1由下一个Xn取代。
然后Y5=P5X6+P5(1-X5),其中P5=P4+FD-1。
P4=0.99并且FD=0.03,那么P5=0.02并且Y5=0.02X6+0.98X5。
然后,只要0≤Pn≤1,使用内插公式Yn=Xn-1+Pn(Xn-Xn-1)计算下一个Yn,其中该内插公式应用至Xn=Xn+1以及Xn-1=Xn,其中Pn=Pn-1+FD。
图5表示在由于负的计算后的频率漂移FD使得内插因子变得小于0的情况下,对应的信号Xn、Xn+1、Yn、TLC和FILTLC的对应计时图,这对应于数据获取设备的本地时钟2太慢的情况。
假设对于Y13,内插因子P13满足0≤P13≤1的要求,使用内插公式Y13=(1-P13)X12+P13X13进行Y13的计算。
例如,P13=0.04并且FD=-0.03。
然后P14=0.01并且Y14=(1-P14)X13+P14X14=0.99X13+0.01X14.
假设在图4中,P15变得小于0,即P15=P14+FD<0。
提供预测装置43来在TLC的每个时钟时隙TSn检查根据公式Pn+1=Pn+FD对于下一个时隙TSn+1计算的内插因子Pn+1是否小于0。所以内插装置在时隙TSn过程中计算对下一个时隙TSn+1可用的内插因子Pn+1。
如果预测装置43为是,即,对于TLC的时隙TS14,P15=-0.02<0,那么时钟滤波装置将TLC的所述时隙TS14划分为在滤波后的时钟信号FILTLC中的第一部分时隙TS14-PART1和在滤波后的时钟信号FILTLC中的第二部分时隙TS14-PART2,不再保持数据Y14,并且基于应用至新的内插因子Pn+1bis=P15bis=P15+1的内插公式Yn=Xn-1+Pn(Xn-Xn-1)计算新的数据Y15,在该内插公式中Xn由Xn+1取代,Xn-1由Xn取代,这意味着计算基于对应于TLC的下一个时隙TSn+1的采样Xn+1和Xn。这意味着将1加入小于0的错误的内插因子Pn+1=Pn+FD,从而得到新的内插因子Pn+1bis=Pn+1+1,该新的内插因子将大于0。
所以P15bis=0.98,并且Y15=(1-P15bis)X15+P15bisX16=0.02X15+0.98X16
然后,只要0≤Pn≤1,使用内插公式Yn=Xn-1+Pn(Xn-Xn-1)计算下一个Yn,其中该内插公式应用至Xn=Xn+1以及Xn-1=Xn,其中Pn=Pn-1+FD(以P15bis开始)。
在第二实施例中,内插函数F的阶是2:
Yn=C0Xn+C1Xn-1+C2Xn-2
可以通过三个因子数字滤波器C0、C1和C2来运行二次内插函数。
对于每个时隙TSn,内插因子Pn增加FD。根据下面的公式计算内插因子P的值Pn:
Pn=Pn-1+FD
对于Yn,通过以下计算内插函数:
第二实施例与第一实施例相比,可以更好地衰减错误。在第二实施例中表示了参考图4和5描述的装置和操作。
在其他实施例中,内插函数的阶可以是3、4或5,并且可以在这些进一步的实施例中表示参考图4和5描述的装置和操作。
内插之后,将从装置4或的第二数据Y发送给抽选滤波器7。产生第二数据Y的内插装置4的输出41连接至抽选滤波器7的输入71。基于第二数据Y,抽选滤波器7在其输出72产生第三数字数据A,该第三数字数据A具有低于第一采样频率FE的第三预设采样频率F3和具有高于第一比特精度BR1的第三预设比特精度BR3,即,在sigma-delta模数转换器3具有所述第一1比特精度BR1的情况下,第三预设比特精度BR3高于1比特。在上述sigma-delta模数转换器3具有所述第一1比特精度BR1的情况下,所述第三比特精度BR3是16比特,并且所述第二比特精度BR2是24比特,所述第三采样频率F3是2kHz,对于采样X和数据Y来说FE=256kHz。抽选滤波器7也具有至少一个输入73,用于接收滤波后的时钟信号FILTLC,其是用于输入数据Y的时钟。
根据本发明的方法能够使用内插进行连续的获取而不具有持续时间限制。在数字抽选滤波器中同步数据流水线。
图6表示的系统具有若干单元100,也被称为节点,每个节点集成了与一个传感器或一组传感器相关联的数据获取设备1,其中每个传感器产生对其相关的传感器的数据输入。在图6中由盒子102示出了传感器,但是传感器可以在每个单元100的外壳内部或者在每个单元100的外壳之外。单元100被称为现场数字化单元(Field Digitilizer Unit,FDU),因为对于其中放置了相关的传感器102的空间现场提供单元100,并且与不同单元100相关联的传感器102以非零的距离分开,即,由与单元100相关联的传感器覆盖的空间场不同于由与另一个单元100相关联的另一个传感器102覆盖的空间场。至少提供电缆200来将参考时钟REF带给单元100。电缆200例如是在单元100之间延伸的电缆的形式,并且具有连接至参考时钟6的一端201。
然后,避免了涉及对每个单元100的大的硬件需求的PLL(锁相环)的使用。并非控制本地时钟,在每个单元100中测量本地时钟2的默认值,并且通过内插采样来数字地补偿该默认值。内插函数的复杂性和消耗变得可以忽略。
同步具有来自输出11的高过采样率的采样X。为了确定精确的FE值,基于具有任意已知频率FREF的参考时钟6来校准采样时钟2、20。为了进行异步采样率转换,对于来自输出11的数据X进行简单内插。
Claims (14)
1.一种数据获取设备,包括:
第一装置,在第一本地时钟(2)提供的不完整采样频率(FE)下提供一系列第一数字采样定时数据(X),
第二校准装置(50),用于通过第二参考时钟(6)测量所述第一本地时钟(2)的频率误差,
第三装置,用于至少基于测量的频率误差来校准所述第一数据,其特征在于,
第一装置(3)包括sigma-delta调制装置(3),用于产生在第一本地时钟(2)提供的所述第一不完整采样频率(FE)下采样的所述一系列第一数字sigma-delta调制数据(X),
第三装置包括至少一内插装置(4)和至少一抽选滤波器(7),内插装置(4)至少基于所述第一系列数据采样(X)的内插来产生第二数字数据(Y)并且补偿测量到的频率误差(FD),并且,抽选滤波器(7)用于基于所述第二数字数据(Y)产生第三数字数据(A)。
2.根据权利要求1所述的数据获取设备,其中所述第一装置包括至少一产生模拟测量数据的模拟传感器(102),并且sigma-delta调制装置(3)是sigma-delta模数转换器(3)的形式,用于将传感器(102)的所述模拟测量数据转换为所述一系列第一数字sigma-delta调制数据(X)。
3.根据权利要求1所述的数据获取设备,其中所述第一装置包括至少一数字传感器,该数字传感器具有所述sigma-delta调制装置以产生所述一系列第一数字sigma-delta调制数据(X)。
4.根据前述权利要求任一项所述的数据获取设备,其中该设备包括控制装置(42),控制装置(42)用于至少基于所述测量到的频率误差(FD)设置所述内插。
5.根据权利要求4所述的数据获取设备,其中,内插装置(4)具有内插函数(F),所述内插函数(F)具有预设的固定阶以及至少一参数(C0、C1、C2),所述第三控制装置(42)包括至少基于所述测量到的频率误差(FD)设置内插函数(F)的所述参数(C0、C1、C2)的装置。
6.根据权利要求5所述的数据获取设备,其中,内插函数(F)的阶小于6。
7.根据权利要求5所述的数据获取设备,其中,内插函数(F)是线性的,并且控制装置(42)根据当前第一采样Xn的以下内插公式设置内插函数(F):
Yn=C0·Xn+C1·Xn-1
其中C0=Pn
C1=1-Pn
其中Xn-1和Xn分别是连续的第一采样(X),
Yn是对于Xn-1和Xn的内插的第三数据(Y),
其中控制装置(42)根据以下公式计算内插因子(P)的值Pn:
Pn=Pn-1+FD
其中,Pn-1和Pn分别是内插因子P的连续值,并且FD是频率漂移。
8.根据权利要求5所述的数据获取设备,其中,内插函数(F)具有等于2的阶,并且控制装置(42)根据当前第一采样Xn的以下内插公式来设置内插函数(F):
Yn=C0·Xn+C1·Xn-1+C2·Xn-2
其中
其中Xn-2、Xn-1和Xn分别是连续的第一采样(X),
Yn是对于Xn-2、Xn-1和Xn的内插后的第三数据(Y),
其中第三控制装置(42)根据以下公式计算内插因子(P)的值Pn:
Pn=Pn-1+FD
其中Pn-1和Pn分别是内插因子P的连续值,并且FD是频率漂移。
9.根据权利要求7或8所述的数据获取设备,其中,控制装置(42)包括根据公式Pn=Pn-1+FD检查内插因子(P)的值Pn是否大于或等于0且小于或等于1的装置(43、44),并且只要内插因子大于或等于0且小于或等于1,控制装置(42)根据内插因子(P)为值Pn的Xn的内插方程设置内插函数(F)。
10.根据权利要求9所述的数据获取设备,其中,控制装置(42)包括根据公式Pn=Pn-1+FD检查内插因子(P)的值Pn对于当前第一采样Xn是否大于1的装置(44),并且如果根据公式Pn=Pn-1+FD对于当前第一采样Xn计算的内插因子(P)的值Pn大于1,那么控制装置(42)通过Pn=Pn-1+FD不将根据内插公式计算的对应的内插后的值Yn用作当前第一采样Xn的第三数据(Y),而是通过根据Pn=Pn-1+FD-1计算的内插因子的值Pn,产生对于当前第一采样Xn根据内插公式计算的值Yn作为与当前第一采样Xn相关的第三数据。
11.根据权利要求9或10所述的数据获取设备,其中,控制装置(42)包括预测装置(43),根据公式Pn+1=Pn+FD检查对于当前第一采样Xn的内插因子(P)的下一个值Pn+1是否小于0,并且如果根据公式Pn=Pn-1+FD对于当前第一采样Xn计算的内插因子(P)的值Pn小于0,那么控制装置(42)通过Pn=Pn-1+FD对于当前第一采样Xn使用根据内插公式计算的对应的内插后的值Yn作为第三数据(Y),并且通过根据Pn=Pn-1+FD+1计算的内插因子的值Pn,产生对于当前第一采样Xn根据内插公式计算的另一个值作为第三数据(Y)。
12.根据前述权利要求任一项所述的数据获取设备,其中,所述一系列第一数字sigma-delta调制数据(X)具有至少为1并且小于等于4的第一预设比特精度(BR1),第二数据(Y)具有所述第一采样频率(FE)以及高于第一比特精度(BR1)的第二预设比特精度(BR2),并且第三数据(A)具有低于第一采样频率(FE)的第三预设采样频率(F3)以及高于第一比特精度(BR1)的第三预设比特精度(BR3)。
13.根据前述权利要求任一项所述的数据获取设备,其中,预设第一本地时钟频率FE的值,使得根据过采样率OR有至少一个传感器提供的信号的过采样高于1,其中过采样率OR被定义为等于:
OR=FE/2·FU
其中FU是由所述传感器提供的所述信号的最高有用频率。
14.一种由计算装置执行的数据获取方法,包括以下步骤:
在第一本地时钟(2)提供的不完整采样频率(FE)下提供一系列第一数字采样定时数据(X),
通过根据第二参考时钟(6)校准第一本地时钟(100)来根据第二参考时钟(6)测量第一本地时钟(100)的频率误差,
至少基于测量到的频率误差来校正第一数据,
其特征在于,
由sigma-delta调制(3)产生在所述第一本地时钟(2)提供的所述第一不完整采样频率(FE)下采样的所述一系列第一数字sigma-delta调制数据(X),
为了补偿测量到的频率误差(FD),对于所述一系列第一数字sigma-delta调制数据(X)进行内插以产生第二数字数据(Y),
对所述第二数字数据(Y)进行抽选以产生第三数字数据(A)。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP09305102A EP2214034B1 (en) | 2009-02-03 | 2009-02-03 | Data acquisition apparatus and method |
EP09305102.7 | 2009-02-03 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101826873A true CN101826873A (zh) | 2010-09-08 |
CN101826873B CN101826873B (zh) | 2014-04-30 |
Family
ID=40790870
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201010112930.8A Active CN101826873B (zh) | 2009-02-03 | 2010-02-03 | 数据获取设备和方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US8260580B2 (zh) |
EP (1) | EP2214034B1 (zh) |
CN (1) | CN101826873B (zh) |
CA (1) | CA2691160C (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105068121A (zh) * | 2015-07-30 | 2015-11-18 | 吉林大学 | 一种井下多节点信号采集同步误差校正装置及方法 |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2546613B (en) * | 2014-06-24 | 2021-01-06 | Cirrus Logic Inc | Systems and methods for clock synchronization in a data acquisition system |
EP2990832A1 (en) | 2014-08-29 | 2016-03-02 | Sercel | Data acquisition apparatus using one single local clock |
EP3345358B1 (en) * | 2015-09-05 | 2023-08-30 | Schlumberger Technology B.V. | Electromagnetic wave pulse synchronization |
US9823368B2 (en) | 2015-12-15 | 2017-11-21 | Sercel | Average clock adjustment for data acquisition system and method |
EP3253001B1 (en) * | 2016-06-02 | 2018-11-14 | Siemens Industry Software NV | Synchronized data transmission from a data acquisition system to a real time digital bus |
EA201992064A1 (ru) * | 2017-03-03 | 2020-02-04 | Шлюмбергер Текнолоджи Б.В. | Система сейсмических датчиков с тактовым генератором на основе микроэлектромеханических систем ("мэмс") |
EP3718089A1 (en) * | 2018-02-13 | 2020-10-07 | Wärtsilä Finland Oy | Apparatus, device and computer implemented method for providing marine vessel data of marine vessel with plurality of sensor devices |
GB201914919D0 (en) * | 2019-10-15 | 2019-11-27 | Bp Exploration Operating Co Ltd | Clock drift |
CN114895353B (zh) | 2022-05-27 | 2023-03-10 | 中国矿业大学 | 井地一体微震监测系统监测单元采集数据授时对齐方法 |
Family Cites Families (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2237552A5 (zh) | 1973-07-10 | 1975-02-07 | Sercel Rech Const Elect | |
FR2318430A1 (fr) | 1975-07-18 | 1977-02-11 | Sercel Rech Const Elect | Dispositif d'acquisition sur le terrain de donnees sismiques avec numeration a pesee aleatoire |
US5355654A (en) * | 1993-04-23 | 1994-10-18 | Stanley Kenneth M | Simulated solid wood slab |
WO1995008221A1 (en) | 1993-09-13 | 1995-03-23 | Analog Devices, Inc. | Digital to analog conversion using nonuniform sample rates |
WO1995031860A1 (en) * | 1993-09-13 | 1995-11-23 | Analog Devices, Inc. | Digital-to-digital conversion using nonuniform sample rates |
US5712635A (en) | 1993-09-13 | 1998-01-27 | Analog Devices Inc | Digital to analog conversion using nonuniform sample rates |
US5512897A (en) * | 1995-03-15 | 1996-04-30 | Analog Devices, Inc. | Variable sample rate DAC |
FR2747795B1 (fr) | 1996-04-19 | 1998-08-14 | Soc Et Rech Et Const Electroni | Procede de compression de donnees d'acquisition sismique |
FR2769369B1 (fr) | 1997-10-08 | 1999-12-24 | Sercel Rech Const Elect | Accelerometre a plaque mobile, avec moteur electrostatique de contre-reaction |
US6957147B2 (en) | 2002-03-12 | 2005-10-18 | Sercel, Inc. | Data management for seismic acquisition using variable compression ratio as a function of background noise |
KR100547831B1 (ko) * | 2003-06-18 | 2006-01-31 | 삼성전자주식회사 | 가변 데이터 전송률에 대응이 가능한 클럭 및 데이터 복원장치 |
FR2889389B1 (fr) | 2005-07-28 | 2008-02-01 | Sercel Sa | Reseau d'acquisition de donnees mixte sans fil et cable |
FR2889385B1 (fr) | 2005-07-28 | 2008-09-26 | Sercel Sa | Reseau d'acquisition de donnees sans fils |
FR2889331B1 (fr) * | 2005-07-28 | 2008-02-01 | Sercel Sa | Appareil et procede de compensation de derive d'une horloge locale utilisee comme frequence d'echantillonnage |
FR2889386B1 (fr) | 2005-07-28 | 2007-10-19 | Sercel Sa | Dispositif et procede de connexion a un reseau sans fil |
FR2938928B1 (fr) | 2008-11-25 | 2011-01-21 | Sercel Rech Const Elect | Systeme d'acquisition de donnees sismiques comprenant des modules associes a des boitiers relies a des capteurs, les modules etant autonomes en alimentation, synchronisation et memorisation. |
-
2009
- 2009-02-03 EP EP09305102A patent/EP2214034B1/en active Active
-
2010
- 2010-01-26 CA CA2691160A patent/CA2691160C/en active Active
- 2010-01-29 US US12/696,152 patent/US8260580B2/en active Active
- 2010-02-03 CN CN201010112930.8A patent/CN101826873B/zh active Active
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105068121A (zh) * | 2015-07-30 | 2015-11-18 | 吉林大学 | 一种井下多节点信号采集同步误差校正装置及方法 |
Also Published As
Publication number | Publication date |
---|---|
EP2214034A1 (en) | 2010-08-04 |
US8260580B2 (en) | 2012-09-04 |
CA2691160A1 (en) | 2010-08-03 |
US20100198561A1 (en) | 2010-08-05 |
CA2691160C (en) | 2016-12-20 |
EP2214034B1 (en) | 2013-04-03 |
CN101826873B (zh) | 2014-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101826873B (zh) | 数据获取设备和方法 | |
CN1904642B (zh) | 用于补偿用作采样频率的本地时钟的漂移的装置和方法 | |
US6377037B1 (en) | Watt-hour meter with digital per-phase power factor compensation | |
US7398411B2 (en) | Self-calibrating time code generator | |
CN105549379A (zh) | 一种基于高精度时间基准触发的同步测量装置及方法 | |
CN103605023A (zh) | 一种合并单元时间特性测量方法及测量装置 | |
CN103178779A (zh) | 一种具有幅度补偿功能的信号发生器及其方法 | |
CN108535772A (zh) | 一种用于井下多节点采集系统时间同步的补偿方法及装置 | |
CN103487648A (zh) | sigma-delta PLL频率测量电路及方法 | |
US20160061972A1 (en) | Data acquisition apparatus using one single local clock | |
EP1756678B1 (en) | Measuring clock jitter | |
CN108474863B (zh) | 用于数据采集系统和方法的平均时钟调整 | |
JP2010164346A (ja) | データロガーの測定時刻補正方法およびデータロガーの測定時刻補正システム | |
CN102636189A (zh) | 用于校正测量信号的测量信号校正装置和方法 | |
US9151687B2 (en) | Measuring apparatus | |
US6373415B1 (en) | Digital phase compensation methods and systems for a dual-channel analog-to-digital converter | |
Tamura et al. | Scale factor calibration of a superconducting gravimeter at Esashi Station, Japan, using absolute gravity measurements | |
US20040094338A1 (en) | System and method for generating a high speed estimated position output for a position encoder | |
Henry | On-line compensation in a digital Coriolis mass flow meter | |
US6304202B1 (en) | Delay correction system and method for a voltage channel in a sampled data measurement system | |
CA2133837C (en) | Full and partial cycle counting apparatus and method | |
US11038510B2 (en) | Oscillator with time error correction | |
CN109274366A (zh) | 配电终端的时钟测试仪及时钟测试方法 | |
US20240171369A1 (en) | An incoherent clocking method | |
Asch | Seismic recording systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |