技术领域
[0001] 本
发明属于岩石物理特性分析技术领域，具体涉及一种幅值和相位自校准的低频岩石物理测量系统。
背景技术
[0002] 岩石物理特性研究的重要性在于它能很好地解释地表地球物理观测结果，为地质正反演提供科学的数据资料，标度地球物理参数。岩石的物理性质主要是由组成岩石矿物的性质、岩石所处的热
力学环境（
温度和压力）以及岩石微构造（孔隙、裂纹）三类因素所决定的。通过对岩石样品物理性质的测量，就能获得地下岩层的整体性质，就能获得
地层下油气、矿物和元素分布的信息。
[0003] 迄今为止，基于岩石中波传播性质的
地震方法是目前地球物理勘探中最主要的方法。岩石中可能产生两类弹性波（
频率在2kHz以下），一类是纵波，也就作P波，其质点运动方向与波传播方向平行；另一类是横波，也叫S波，它的质点运动方向与波传播方向垂直。岩石受到外力而在与该力平行和垂直的两个方向上发生的形变就是P波和S波在岩石中传播的宏观体现。
[0004] 实验室测量P波和S波以及
波速衰减的方法很多，但各种方法能够测量的
频率范围不同，每种方法只适用于一定的频率范围。
[0005] 描述岩石弹性参数的量有
拉梅常数λ，
剪切模量μ，体积模量K，
杨氏模量E，泊松比γ，弹性波波速VP、VS等。事实上，这些参数中只有两个量是相互独立的，并且他们均可以由
应力σ、纵向应变ε⊥和横向应变ε||三个量推导出。因此，岩石弹性参数的测量就转化为应力、纵向应变和横向应变的测量。
[0006] 由于岩石并不是理想的弹性体，它还表现出一定的非弹性性质。采用岩石的品质因子Q（与衰减系数α成反比）来表征岩石的非弹性性，Q值越小，非弹性特性就越突出，对于完全弹性体，Q=∞。在低频岩石物理测量中，我们给岩石施加的并非一个恒定的压力，而是一个交变的应力，即σ=σ0sinωt。岩石的非弹性特征就通过纵向应变或横向应变与应力σ之间的
相位差θ表现出来。因此，岩石非弹性参数的测量就化为应力与应变（纵向或横向均可）之间相位差的测量。
[0007] 现有的岩石应力-应变测量方法没有考虑
电子系统在测量一个
信号的幅值和相位时，往往会受诸多因素影响，如温度、时间以及信号频率的改变都会导致测量值与真实值有很大偏差。休斯顿大学岩石物理研究所关于低频岩石物理测量的实验结果表明，固定一个信号的幅值和相位，当信号频率变化时，
电路部分实际测量的幅值和相位也随之变化，尤其是相位甚至发生了反转。
发明内容
[0008] 针对
现有技术的以上
缺陷或改进需求，本发明的目的是在于提供了一种幅值和相位自校准的低频岩石物理测量系统，该系统能消除环境参数变化带来电子测量系统幅值和相位的漂移，提高系统幅值和相位测量
精度。
[0009] 一种幅值和相位自校准的低频岩石物理测量方法，具体为：
[0010] （1）搭建测量环境步骤：
[0011] 将圆柱形岩石置于金属
支撑台上并处于围压容器内，向围压容器内
泵入传压介质，分别在岩石的侧面纵向和横向多点布置
传感器；金属支撑台的上部分为截面直径与岩石相同的圆柱形，在金属支撑台的圆柱形侧面纵向多点布置传感器；
[0012] （2）标定步骤：
[0013] 在2Hz～2kHz频率范围内选取多个频率点，分别将对应这些不同频率但相同相位和幅值的交变
电压信号输入各位移传感器，采集各位移传感器的输出交变电压信号，对输出交变电压信号放大、
模数转换得到标定电压信号；针对每一个位移传感器,计算各频率对应的标定电压信号幅值V1和相位θ1，构建标定步骤频率点与标定电压信号幅值和相位间的映射关系；
[0014] （3）测量步骤：
[0015] 向各位移传感器输入直流电压信号，采用标定步骤所用的多个交变电压信号作为激振信号产生交变力,将交变力施加到岩石上，岩石产生形变，各位移传感器采集岩石形变量并转换为测量电压信号，对测量电压信号进行放大和模数转换；所述交变力的频率、幅值和相位与标定步骤所用的交变电压信号相同；针对每一个位移传感器，计算各频率对应的测量电压信号幅值V2和相位θ2，构建测量步骤频率点与测量电压信号幅值和相位间的映射关系；
[0016] （4）幅值和相位自校准步骤：
[0017] 对于每一个传感器，计算其对应的校准后交变电压信号的幅值Vin＝VrefV2V1及相位θin＝θ2-θ1+θref，其中，Vref和θref为交变电压信号的幅值和相位；
[0018] （5）岩石属性参数计算步骤：
[0019] 依据校准后交变电压信号的幅值和相位计算岩石弹性参数和非弹性参数。
[0020] 进一步地，所述步骤（5）岩石属性参数计算步骤具体为：
[0021] 计算幅值和相位均值：计算金属支撑台的多个传感器校准后幅值均值 和相位均值 岩石侧面纵向设置的多个传感器校准后幅值均值 和相位均值 岩石侧面横向设置的多个传感器校准后幅值均值 和相位均值
[0022] 计算岩石弹性参数：岩石杨氏模量 岩石泊松比岩石体积模量Kr＝Er3(1-2λr)，岩石拉梅常数γr＝Erλr(1+λr)(1-2λr)，岩石纵波速度岩石横波速度 式中，Em为金属支撑
台杨氏模量，ρ为岩石
密度；
[0023] 计算岩石非弹性参数：岩石衰减系数
[0024] 进一步地，所述步骤（1）测量环境步骤中位移传感器设置方式具体为：在岩石
侧壁纵向和横向分别均匀布置四个位移传感器，在金属支撑台的圆柱形侧壁纵向均匀布置四个位移传感器。
[0025] 本发明方法的技术效果体现在：
[0026] 本发明考虑到温度、时间、频率等因素对测量系统电路部分的影响，对多个不同频率点采用实时标定和测量的方式，对测量结果进行幅值和相位自校准，能够准确地测量电压信号的幅值和相位。该方法能消除环境参数变化带来电子测量系统幅值和相位的漂移，提高系统幅值和相位测量精度。
[0027] 进一步地，在岩石属性参数计算步骤中，采用自校准后的交变电压信号的幅值和相位均值作为参数计算依据，提高了测量精度。
[0028] 进一步地，对位移传感器的布置进行了优化，综合考虑岩石在各个方位存在不同程度的
变形以及计算量，优选分别横向和纵向均匀布置多个（优选四个）位移传感器。
[0029] 为实现本发明的另一目的，本发明提供了实现上述方法的低频岩石物理测量系统，包括
[0030] 信号激励模
块，其输入端连接
数据采集处理模块的信号控制输出端，其输出端连接标定与测量切换模块和数据采集处理模块的输入端，用于接收来自数据采集处理模块的频率相位幅值
控制信号，产生相应的交变电压信号，并输出给标定与测量切换模块和数据采集处理模块；
[0031] 机械模块，包括金属支撑台、围压容器、位移传感器和激振器，金属支撑台用于支撑岩石，围压容器用于为岩石和金属支撑台提供密闭带压空间，岩石侧壁纵向和横向多点布置有位移传感器，金属支撑台侧壁纵向多点布置有位移传感器，激振器的激励端连接标定与测量切换模块的输出端；在标定模式下，位移传感器接收来自信号激励模块的交变电压信号，将其作为标定电压信号输出给数据采集处理模块；在测量模式下，激振器以来自信号激励模块的交变电压信号为激励源产生交变应力作用在岩石和金属支撑台上，传感器采集表征岩石形变量信号并转换为测量电压信号传送给数据采集处理模块；
[0032] 标定与测量切换模块，其输入端连接信号激励模块的输出端，其控制端连接数据采集处理模块的模式选择输出端，用于接收来自采集处理模块输出的模式选择信号，以使得机械模块在标定模式和测量模式间切换；
[0033] 数据采集处理模块，其输入端连接机械模块的传感器的输出端，其信号控制输出端连接信号激励模块的输入端，其模式选择输出端连接标定与测量切换模块的控制端，用于向信号激励模块传送频率相位幅值控制信号，向标定与测量切换模块传送模式选择信号，用于接收来自机械模块的标定电压信号或测量电压信号；在标定模式下，对于每个传感器，构建标定模式频率点与标定电压信号幅值和相位间的映射关系；在测量模式下，对于每个传感器，构建测量模式下频率点与测量电压信号对应幅值和相位间的映射关系，结合标定模式和测量模式下构建的映射关系进行幅值和相位的自校准，并依据自校准结果计算岩石属性参数。
[0034] 进一步地，所述标定与测量切换模块包括
电池、继电器驱动模块和两个继电器；电池的输出端连接第一继电器的输入端，第一继电器的输出端连接位移传感器的输入端；第二继电器为双刀双掷式，其第一输出端连接位移传感器的输入端，其第二输出端连接激振器的输入端，其输入端连接取信号激励模块的输出端；两继电器的控制端均连接继电器驱动模块的继电器控制信号输出端，继电器驱动模块的控制端连接数据采集处理模块的模式选择输出端。
[0035] 进一步地，所述数据采集处理模块包括依次连接的前置放大电路、模数转换电路、FPGA和计算机；前置放大电路的第一输入端连接信号激励模块的输出端，前置放大电路的第二输入端连接位移传感器的输出端，FPGA的模式选择输出端连接标定与测量切换模块的控制端，FPGA的信号控制输出端连接信号激励模块的输入端。
[0036] 进一步地，所述计算机内加载有计算系统，计算系统包括
[0037] 标定映射关系构建模块，用于在标定模式下，对于每个传感器,计算各频率对应的标定电压信号幅值V1和相位θ1，构建标定模式频率点与标定电压信号幅值和相位间的映射关系；
[0038] 测量映射关系构建模块，用于在测量模式下，对于每个传感器，计算各频率对应的测量电压信号幅值V2和相位θ2，构建测量模式频率点与测量电压信号幅值和相位间的映射关系；
[0039] 幅值和相位自校准模块，用于对于每一个传感器，计算其对应的校准后交变电压信号的幅值Vin＝VrefV2V1及相位θin＝θ2-θ1+θref，其中，Vref和θref为交变电压信号的幅值和相位；
[0040] 岩石属性参数计算模块，包括
[0041] 幅值和相位均值计算子模块，用于计算金属支撑台的多个传感器校准后幅值均值和相位均值 岩石侧面纵向设置的多个传感器校准后幅值均值 和相位均值岩石侧面横向设置的多个传感器校准后幅值均值 和相位均值
[0042] 岩石弹性参数计算子模块，用于计算岩石杨氏模量 岩石泊松比 岩石体积模量Kr＝Er/3(1-2λr)，岩石拉梅常数γr＝Erλr/
(1+λr)(1-2λr)，岩石纵波速度 岩石横波速度
式中，Em为金属支撑台杨氏模量，ρ为岩石密度；
[0043] 岩石非弹性参数计算子模块，用于计算岩石衰减系数
[0044] 本发明测量系统的技术效果体现在：
[0045] 本发明能够模拟地层条件，即在低频段（2Hz～2kHz）频率范围内对岩石样品的弹性和非弹性参数进行测量：激振器产生一定频率的纵向交变应力，使岩石样品和基准块（即金属支撑台）产生微小形变，位于岩石样品和基准块表面的光纤式位移传感器将微小形变转化为微小
电信号，最后经多通道数据采集系统获取该频率点的测量结果。对多个不同频率点采用实时标定和测量的方式，对测量结果进行幅值和相位自校准，能够准确地测量电压信号的幅值和相位。本发明利用标定与测量切换模块方便地在标定与测量模式进行切换，通过数据采集处理模块对多个不同频率点进行实时标定和测量，
激励信号由数据采集处理模块控制，并经低通滤波和功率放大，可以实现激励信号的频率可调、相位可调、幅度可调、功率可调等功能。本发明测量系统能够准确地测量电压信号的幅值和相位，消除环境参数变化带来电子测量系统幅值和相位的漂移，提高系统幅值和相位测量精度。
[0046] 进一步地，所述标定与测量切换模块采用两个继电器配合方式实现标定与测量模式的切换，结构简单，容易实现，响应快。
[0047] 进一步地，所述数据采集处理模块以FPGA为核心，主要完成对继电器的
开关控制、DDS电路的频率控制、模数转换电路的时序控制以及和计算机的通信。FPGA与计算机可采用USB2.0、100M Ethernet和RS230多
接口相连，数据传输速率可达数十兆每秒
附图说明
[0048] 图1为一种幅值和相位自校准的低频岩石物理测量系统工作
流程图。
[0049] 图2为一种幅值和相位自校准的低频岩石物理测量系统结构
框图。
[0050] 图3为电动式激振器结构图。
[0051] 图4为光纤式位移传感器结构图。
[0052] 图5为经DPSD
算法处理前后的信号对比图，图5（a）为处理前，5（b）为处理后。
[0053] 其中：1、DDS（直接数字频率合成器）；2、低通
滤波器；3、功率放大电路；4、电池；5、取样电路；6、激振器；7、围压容器；8、岩石样品；9、基准
铝块；10、光纤式位移传感器；11、压力泵；12、第一继电器；13、第二继电器；14、继电器驱动模块；15、前置放大电路；16、模数转换电路；17、FPGA；18、计算机；19、
弹簧；20、壳体；21、顶杆；22、磁
钢；23、磁极；24、线圈；25、
铁芯；26、岩石样品；27、基准铝块；28、金属
套管；29、光纤；30、相关电路。
具体实施方式
[0054] 为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及
实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
[0055] 参见图1,本发明测量方法包括如下步骤：
[0056] （1）搭建测量环境
[0057] 将圆柱形岩石置于金属支撑台上并处于围压容器内，向围压容器内泵入传压介质，分别在岩石和金属支撑台上设置位移传感器。具体而言，由于岩石变形存在
各向异性，因此可在其侧面多点布置传感器，推荐在其侧壁纵向均匀布置四个位移传感器，用于采集岩石表面的纵向形变，在其侧壁横向均匀布置四个位移传感器，用于采集岩石表面的横向形变。在本发明中，金属支撑台作为基准，因此将其上部分加工为直径与岩石相同的圆柱形，按照与岩石相同的方式布置位移传感器。
[0058] （2）标定
[0059] 在2Hz～2kHz频率范围内选取n个频率点fk，k=1,…,n，分别将对应这些不同频率但相同相位和幅值的交变电压信号输入位移传感器，采集位移传感器的输出交变电压信号，对输出交变电压信号放大、模数转换得到标定电压信号。计算各频率对应的标定电压信号幅值和相位，构建未发生形变时频率点与标定电压信号对应幅值和相位间的映射关系，以备后续测量过程中的幅值和相位自校准。
[0060] 频率点选取越多，测量结果越准确，但增加了计算量，实际应用中可根据精度和时间要求进行调整，例如按每十倍频程取十个频率点。
[0061] （3）测量
[0062] 向各位移传感器输入直流电压信号，采用标定所用的n个交变电压信号作为激振信号产生并向岩石施加交变力，岩石产生形变，各位移传感器采集岩石形变信号，对岩石形变信号放大、模数转换得到测量电压信号。交变力的频率、幅值和相位与标定所用的交变电压信号相同。针对每一个传感器，计算各频率对应的标定电压信号幅值和相位，构建岩石形变下频率点与测量电压信号对应幅值和相位间的映射关系。
[0063] （4）幅值和相位自校准
[0064] 令作为激振源的交变电压信号的幅值为Vref和相位为θref，标定电压信号的幅值为V1和相位为θ1，测量电压信号的幅值为V2和相位为θ2，
[0065] 计算校准后交变电压信号的幅值Vin＝VrefV2V1，相位θin＝θ2-θ1+θref。
[0066] 由于Vin和θin表达式中的参量都是系统实时测得的，消除了温度、时间和频率等因素的影响，因而待测信号的幅值和相位实现了自校准。
[0067] 按照上述自校准方式得到每个传感器对应的校准后交变电压信号的幅值和相位。
[0068] （5）计算岩石弹性参数和非弹性参数
[0069] 求取金属支撑台的多个传感器校准后幅值均值 和相位均值 ，
[0070] 岩石侧面纵向设置的多个传感器校准后幅值均值 和相位均值 岩石侧面横向设置的多个传感器校准后幅值均值 和相位均值
[0071] A.弹性参数计算
[0072] 在测量过程中，岩石和金属支撑台在纵向上受到的应力的幅值一样，于是岩石和金属支撑台的杨氏模量为：Er＝σεv，Em＝σεr(Em为金属支撑台的杨氏模量，是已知量)。由于位移传感器
输出电压幅值与所测形变量成正比，即 那么同样地，有岩石的泊松比为： 岩石其它的弹性参数均可
由Er和λr导出：体积模量Kr＝Er3(1-2λr)，拉梅常数γr＝Erλr(1+λr)(1-2λr)，纵波速度 横波速度 ρ为岩石密度。
[0073] B.非弹性参数计算
[0074] 施加的交变应力是周期性的振动波，先通过岩石然后再传送到金属支撑台会产生一定的延时，因此 和 能在一定程度上表征该岩石特性，从而可计算出岩石的衰减系数为：
[0075] 基于上述方法，本发明提供了一种低频岩石物理测量系统，包括信号激励模块A、机械模块B、标定与测量切换模块C和数据采集处理模块D。
[0076] 信号激励模块A用于产生激励交变电压信号，其包括依次电连接的交变信号发生器1、
低通滤波器2、功率放大电路3和取样电路5。DDS1的输入端连接数据采集处理模块D的信号控制输出端，取样电路5的输出端连接数据采集处理模块D的输入端和测量切换模块C的输入端。交变信号发生器1接收来自数据采集处理模块D的频率相位幅值控制信号，产生相应的交变电压信号，交变电压信号经过低通滤波器2滤波整形、功率放大电路放大后输出，取样电路5对其取样后输出给数据采集处理模块D。
[0077] 按照本发明的一种实施方式，交变信号发生器1采用直接数字频率合成器（DDS），DDS具体选型是ADI公司的AD9850。AD9850是一款高度集成的器件，采用先进的DDS技术，内置一个高速、高性能
数模转换器和比较器，共同构成完整的数字可编程频率合成器和时钟发生器。DDS输出的
正弦波信号实际上是片内DAC输出的阶梯正弦波，再加上参考时钟频率和器件本身噪声等等许多因素的影响，高频杂散噪声很大。为了得到一
频谱纯度很高的信号，采用5阶截止频率为42MHz椭圆低通滤波器。功率放大电路用来将正弦波信号作进一步的功率放大以驱动激振器和光纤式位移传感器。采用甲类功率放大电路，可以使信号的非线性失真最小。
[0078] 机械模块B包括激振器6、围压容器7、金属支撑台9、位移传感器10和压力泵11。金属支撑台9用于放置岩石8，围压容器7连接压力泵11，用于对金属支撑台9和岩石8形成一个密闭带压的空间，以模拟地层高压封闭环境。在金属支撑台9和岩石8侧壁设置位移传感器10，位移传感器10的输出端连接
信号处理控
制模块，用于将采集的岩石形变信号传送给信号处理
控制模块。激振器6的输入端连接标定与测量切换模块C的输出端，若激振器6接收到标定与测量切换模块C的启动信号，则表明进入测量模式，向岩石产生交变力，否则，表明进入标定模式，暂停工作。金属支撑台9可采用铝、铝
合金等金属制备。激振器
7可采用电动式的永
磁铁激振器，如图2所示，当线圈24通过经功率放大后的交变
电流时，根据
磁场中载
流体受力的原理，线圈将受到与电流成正比的点动力作用，此力通过顶杆传到被测对象，即为激振力。激振力作为纵向的交变应力施加在岩石样品和基准铝块上。激振器只在测量模式下工作。
液压泵和围压容器内用
煤油作为传压液体，最高可以产生1GPa的围压。岩石样品、基准铝块和光纤式位移传感器的光纤都放置在围压容器内。
[0079] 按照本发明的一种实施方式，激振器7采用的是电动式永磁铁激振器，根据图3所示，当线圈通过经功率放大后的交变电流时，根据磁场中载流体受力的原理，线圈将受到与电流成正比的点动力作用，此力通过顶杆传到被测对象，即为激振力。激振力作为纵向的交变应力施加在岩石样品和基准铝块上。激振器只在测量模式下工作。
[0080] 位移传感器10可采用光纤式位移传感器、金属应变片等，采用光纤式位移传感器测量岩石样品和基准铝块的形变，可以承受更高的围压。光纤式位移传感器的机构参见图4，主要由金属套管28和光纤29构成。按照本发明的一种实施方式，光纤式位移传感器一共有12个，如图3所示，分别测量岩石在0°、90°、180°、270°四个方向上的纵向应变和横向应变以及铝块在0°、90°、180°、270°四个方向上的纵向应变。光纤式位移传感器在测量模式下的激励信号是电池提供的直流信号，而在标定模式下的激励信号则是信号激励模块提供的交流信号。
[0081] 标定与测量切换模块C用于标定与测量模式的切换，包括电池4、继电器驱动模块14和两个继电器12，13。继电器12的输出端连接位移传感器的输入端；继电器13为双刀双掷式，其第一输出端连接位移传感器的输入端，第二输出端连接激振器6的输入端，输入端连接取样电路5的输出端；两继电器12，13的控制端均连接继电器驱动模块14的继电器控制信号输出端，继电器驱动模块14的控制端连接数据采集处理模块D的模式选择输出端。电池4的输出端连接继电器12的输入端，用于给继电器12供电。继电器驱动模块本质上是一个功率放大电路，数据采集处理模块D输出的继电器控制信号经继电器驱动模块作用在继电器的线圈上，使继电器实现断开或合上的操作。继电器控制信号表明当前为标定模式，则驱使继电器12断开，继电器13与位移传感器10之间连通，继电器13与激振器之间断开，取样电路5采集的交变电压信号为位移传感器10供电；继电器控制信号表明当前为测量模式，则驱使继电器12闭合，继电器13与位移传感器10之间断开，继电器13与激振器6之间连通，取样电路5采集的交变电压信号作为激振器6的激励源，激振器6相应产生交变应力。
[0082] 数据采集处理模块D作为整个系统的核心枢纽，用于各模块的协调控制和
数据处理。其包括依次连接的前置放大电路15、模数转换电路（ADC）16、FPGA17和计算机18。前置放大电路15的第一输入端连接取样电路5的输出端，前置放大电路15的第二输入端连接位移传感器10的输出端，FPGA17的模式选择输出端连接继电器驱动模块14的输入端，FPGA17的信号控制输出端连接DDS1的输入端，将特定频率相位幅值的控制信号传送给DDS1。
[0083] 按照本发明的一种实施方式，前置放大电路共有13路，分别接收来自取样电路5的参考信号和12路来自光纤式位移传感器的
输出信号，并将它们放大1000倍以便ADC16采集。运放芯片采用ADI公司的
仪表放大器AD8429，AD8429是一款超低噪声（1nV/√Hz）
仪表放大器，增益为1000时仍能保证150kHz的带宽和90dB的共模抑制比。其优异的不失真性能，能够用在振动分析等要求苛刻的应用中。
[0084] 按照本发明的一种实施方式，模数转换电路同样有13路，分别将13路前置放大电路的输出信号转化为
数字信号供FPGA17缓存。模数转换芯片采用TI公司的ADS1252，ADS1252是一款高精度、24位串行Σ-Δ型ADC。每个通道采用24位串行Σ-Δ型ADC
采样，能实时改变
采样频率。根据输入不同的时钟频率可以实现从10Hz到40kHz的数据输出速率，并且其内置
数字滤波器，因此该芯片非常适合本发明低频（2Hz～2kHz）测量系统的数据转换。
[0085] FPGA作为整个测量系统的核心，与模数转换电路、继电器驱动电路、DDS电路和计算机相连，主要完成对继电器的开关控制、DDS电路的频率控制、模数转换电路的时序控制以及和计算机的通信。按照本发明的一种实施方式，FPGA芯片是ACTEL公司的A3PE3000，A3PE3000含有112个4608bit的RAM块，根据本测量系统的需要，可以将其配置成13个深度1536×宽度24位的RAM，每个RAM可以缓存来自模数转换电路的3072个数据，然后通过USB、100M Ethernet或RS230接口上传至计算机，数据上传完后开始下一次的采集。FPGA与计算机可采用USB2.0、100M Ethernet和RS230多接口相连，数据传输速率可达数十兆每秒。
[0086] 计算机主要完成的工作是接收FPGA采集的数据，并对数据进行相关算法（DPSD）处理，采用DPSD对采集数据进行处理，将微弱信号从较强背景噪声中提取出来时，其幅值和相位仍能保持较高精度。在
信噪比为-40dB时，经DPSD处理前后的信号如图5（a）和（b）所示，可以明显看出DPSD有很好地噪声抑制能力。处理完毕后，得到每个通道31个频率点的信号幅值和相位，计算出本测量系统需求的弹性参数和非弹性参数并绘制出各个参数2Hz～2kHz的频率特性曲线。
[0087] 本发明系统的工作过程具体为：
[0088] FPGA发出模式选择信号切换继电器，将功率放大电路3通过取样电路5连接位移传感器10，此时进入标定模式，功率放大电路3输出的交变电压信号输入位移传感器10，位移传感器10输出的交变电压信号通过放大、模数转换后得到标定电压信号。计算机18计算各频率对应的标定电压信号幅值和相位，构建未发生形变时频率点与标定电压信号对应幅值和相位间的映射关系，以备后续测量过程中的幅值和相位自校准。
[0089] 再次切换继电器，将功率放大电路3通过取样电路5连接激振器6，此时进入测量模式，激振器给岩石和基准铝块提供一个纵向的交变应力，使其发生微小形变，光纤式位移传感器将此微小形变转化为电压信号，经数据采集系统采集并上传至计算机，待当前频率点采集完毕，这就完成了当前频率点测量。接着按每十倍频程取十个频率点的方法，对其余频率点（到f=2kHz为止）进行测量，所有频率点的数据都采集完毕后，计算机对各个频率点的数据进行DPSD（数字相敏检波）算法处理，计算出各个频率点的幅值和相位，然后根据同一个频率点的标定值和测量值，求出相应的弹性参数和非弹性参数并绘制它们的频率特性曲线。
[0090] 本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何
修改、等同替换和改进等，均应包含在本发明的保护范围之内。
