传统的地震仪在地震监测和研究中发挥着重要作用[1-3]。但在大地震发生时，震源附近的地震仪可能会受到地震波冲击导致振幅满格而难以完整地记录地震波信号。同时，大地震造成的近场地表永久性变形及地面倾斜会导致地震仪记录的地震波信号产生扭曲、变形。而高频GPS观测资料具有能够获取宽频带信号、可以直接获取高精度水平位移信息、不会出现限幅现象以及对地面倾斜不敏感等优势。因此，充分利用高频GPS观测资料中记录的地震信息已被越来越多的学者所重视，GPS地震学[4]随之应运而生。

峰值荷载后在CD段卸载时，以K56为卸载刚度，沿56进行卸载，若卸载后反向加载未过峰值点，则沿路径6C′D′进行加载，若反向加载超过峰值点，则沿路径65′D′进行加载；在C′D′段的卸载情况与CD段相似。

文献[5]利用1 Hz的GPS观测数据获取了FAIR测站记录到的2002年的Denali地震波形，并与地震仪积分结果比较，两者具有较好的一致性。此后，越来越多的学者将高频GPS技术广泛应用于其他地震的研究，并取得了大量的研究成果[6-8]。由GPS获得的同震静态永久变形信息(static offset)可用于约束震源参数如震级、震源深度等的反演，地震波信息可作为地震仪、强震仪等观测数据的重要补充，用于反演地震断层破裂过程。但是，震后的高频GPS观测资料包含了由地震引起的静态永久变形和地震波的详细信息。当需要永久性变形信息反演震源参数时，地震波信息会影响准确的获得永久性变形，反之，当需要地震波信息反演断层破裂过程时，永久性变形信息也会对地震波信息的准确性产生影响。因此，实现两者的有效分离是准确获取两种信息进而反演地震参数的关键因素。以往为了从GPS观测资料中获取同震静态永久变形，通常采用震前和震后一段时间的资料进行计算，并需要剔除地震发生或以后一段时间的数据以避免地震波位移影响。文献[9]利用卡尔曼滤波方法估计了2003年日本Tokachi-Oki地震主震和震后静态永久性变形效应，但是在应用该方法时，为了避免动态位移对滤波结果的影响，在估计中必须剔除地震发生时的观测数据，这样不仅增加了数据处理的工作量，而且也会丢失地震波信息。文献[10]利用小波分解重构方法提取高频GPS数据静态永久性变形信息，但是在利用小波分层技术时，无法完整保存地震波到时及振幅等信息，而到时和振幅信息对于准确确定震中和震级等参数具有重要意义。

为充分发挥高频GPS资料的优势，本文旨在研究在快速获取高频GPS数据静态永久变形信息的同时，如何有效分离并完整保存高频GPS资料中记录的地震波信息。本文选用平滑先验信息方法(SPM)[11-14]来实现地震引起的高频GPS观测结果中静态形变和地震波信息的分离。该方法是一种高通滤波器，对具有周期特性的高频信号最为敏感，能够有效地识别具有周期特性的高频周期信息，从而实现周期信号与非周期信号的有效分离，并且不需要剔除数据。由于地震波具有一定周期或准周期性而表现为高频信息，永久性变形则不具有周期性，正是针对高频GPS资料中地震波与永久性变形的这一特点，本文将SPM方法引入高频GPS数据后处理中，分别获取了静态永久性变形信息和地震波信息。

1 数据处理与方法

1.1 数据及地震位置

本文采用的1 Hz的GPS观测数据由美国地学研究组织UNAVCO网站提供(http:∥www.unavco.org/data/data.html)。本文选取了近几年的3次地震，地震的详细信息见表1。本文依照两个原则，每个地震对应选取了两个观测站：①测站到震中的距离；②在地震发生前后，观测站的数据记录是否完整。测站分布与震中位置分布情况如图1所示。

由图1和图2可知，全国大部分地区Cd污染较为严重，而Cr和Ni污染较轻，这与表4表5评价结果一致。总体而言，除 Zn之外，非工业区和工业区的重金属污染严重区域有相同的分布规律，即人口较多、工业集中、经济发达的中东部地区的降尘中重金属含量普遍高于人口稀少、发展滞后的西部地区。不同之处在于，工业区中重金属的污染情况更为严重且区域多集中在中西部重工业基地以及有色金属矿区和冶炼区。

表1 地震详细信息(USGS)

Tab.1 The detail information of earthquakes (USGS)

地震名称时间/UTC震中位置震源深度/km震级所选测站震中距/km2010⁃04⁃04EIMayor⁃Cucapah地震22：40：4232．286°N，115．295°W10．07．2P496P74458．863．62012⁃09⁃05哥斯达黎加地震14：42：0710．085°N，85．315°W35．07．6QSECSAJU27．643．32011⁃03⁃11东日本大地震主震05：46：2438．297°N，142．373°E29．09．1余震06：15：4036．281°N，141．111°E42．67．9余震06：25：5038．058°N，144．590°E18．67．7MIZUUSUD142．4428．9

1.2 高频GPS数据处理

本文采用Bernese 5.2软件进行动态精密单点定位，卫星轨道及钟差采用欧洲定轨中心(CODE)提供的15 min间隔的精密星历与5 s间隔的精密卫星钟差[15]，并且采用IERS推荐的最新的天线相位中心模型、大气潮汐及海潮模型，卫星截止高度角设置为3°。同时，估计了载波相位模糊度、卫星钟差及接收机钟差。电离层误差通过采用双频电离层无关组合消除，天顶对流层延迟采用WET_GMF模型估计。计算得到的动态位移解仍然受多路径效应等周期性噪声的影响，本文采用改进的恒星日滤波(modified sidereal filtering,MSF)方法[16]进行处理。对P496和P744测站选用震前2 d(P496和P744测站震前只有2 d的数据)，QSEC、SAJU、MIZU和USUD测站选用震前4 d相同恒星日段的时间序列进行加权平均，得到背景噪声。图2以MIZU测站N方向为例，对比了分别采用震前2 d和震前4 d的观测数据计算得到的背景噪声。其中绿色曲线分别对应2011-03-07(年积日066)到2011-03-10(年积日069)动态位移时间序列，红色曲线为滑动平均后的结果，消除了高频噪声。蓝色曲线即为所得的背景噪声。两组背景噪声的相关系数达到了90.4%，从图中也可以看出，在震前4 d该台站周围环境稳定，没有受到明显的噪声事件干扰。从图2结果来看，由于背景噪声只与观测台站质量和周围环境的噪声影响有关，在没有突然噪声的干扰下，选用震前2～4 d的数据处理的结果是基本相同的。从地震当天的动态位移时间序列中扣除背景噪声，扣除后的坐标时间序列基本上不再含有多路径效应等周期性噪声的影响。

1.3 信号分离

这在世界文物考古发掘史上。都是绝无仅有的。这个奇观说明，地处发掘地的汉江流域古老文明，历经悠悠，一脉长流，五千年从未断绝。

（1）水基淬火剂淬火异常现象可能是由于淬火冷却介质老化和浓度改变引起的，因此应该及时检测浓度变化并做出调整。

对于任意时间序列z，均可表示为z=[z1 z2… zn]T。其中，n为采样点个数。z可看作近似稳态的周期分量和非周期趋势项之和

z=zstat+znonstat

(1)

式中，zstat表示近似稳态的周期分量(即地震波)；znonstat表示非周期趋势分量(即静态永久变形)，可表示为

znonstat=Hθ+ε

(2)

(1) 利用GPS静态天解解算地震引起的静态永久变形，不仅需要几天时间才可获得稳定解，而且解算结果包含了主震、余震及震后断层余滑位移，仅靠GPS静态天解结果是难以区分的。而高频GPS数据可以清晰地记录地震引起的静态永久变形和高频地震波信息，采用SPM方法，无须剔除地震时段的数据，因而在完整保留地震波信号的同时，还可以迅速提取静态永久变形信息。分离得到的静态永久变形信息可用于震源参数如震级、震源深度等的反演，而高频地震波信息可以结合地震仪、强震仪等观测资料反演断层破裂过程，为震后快速反应及灾后评估提供保障。

(3)

式中，λ是正则化参数;Dd为第d次导数操作的离散逼近。显然，式(3)是在普通最小二乘方法的基础上增加了一个条件，属于条件极值问题。极值条件的目的是使式(3)的解在满足取得极小值的同时，使得取值也尽可能小。这样就可以在推估线性项Hθ的过程中，加入一些已知的先验信息。式(3)的解可表示为

(4)

(5)

式中，为非周期分量的估计值。观测矩阵H的选择可以包含时间序列z的一些已知信息，例如使用一些高斯型函数或三角函数，这里选取单位阵I作为观测矩阵，即H=I∈Rn×n，详细参见文献[13]。式(3)中正则化部分可以理解为求解正则化矩阵Dd零空间里的解，二阶差分矩阵的零空间包含了所有的一阶曲线，因此D2是用来估计非周期项的一个不错的选择[11,13]。二阶差分矩阵D2∈R(n-2)×n形式如下

(6)

有了这些特殊选择后，该方法称为平滑先验信息方法(smoothness priors method,SPM)[11-14]。这样，周期信号可以表示为

(7)

1.3.2 频率响应及平滑参数的选取

以上现象表明，即使我们可以分析出较早的字形，这个字形却也未必是它的本字。连绵词其本质是以音表义，有些表面看来意义与字形相关，事实却是受到了汉字形义相关的构形规律的影响。比如，表示“辽阔，不清”貌，有从“目”的“眇莽”“眇漠”“眇茫”，有从“水”表义的“渺茫”“渺莽”“渺漭”“淼茫”等词，事实上，这些连绵词只是记录表示“辽阔、迷濛”意义的语音。

台肇地区注水系统节能降耗工作以注配间降压节电为目标，通过单井、管线、注配间“两点一线”开展治理工作，按照地面设备改造、降压增注措施、地质方案优化相结合的思路，深化注水系统能耗挖潜。鉴于注水井可以按照配注完成情况分为两类：多数的正常井，注入压力高低有别；少数的欠注井，注入压力接近泵压。针对前者，注重从注配间地面设备改造入手，通过采取分压注水、安装变频器、管线改造等措施，力图将注入压力接近的井接入统一柱塞泵；针对后者，注重利用采油工程和油藏工程的手段从单井入手，通过利用酸化解堵、化学洗井和方案优化等措施，力图将高压欠注井降压。

式(7)可以写为其中相当于一个时变高通滤波器，对于每一个离散时间点，L的频率响应(以其行的傅氏变换形式获得)如图3(a)所示，本文给出了前半部分的频率响应，后半部分是完全相同的。可以看到，滤波器基本为一常数，仅对起始和结束的信号做了不同处理。其滤波效果仅对数据的起始和末尾有扭曲现象，最大程度的避免截断效应的影响。对于信号z的第一个和最后一个元素，滤波效果是衰减的，因此避免了数据终点的扭曲。平滑参数λ对于频率响应的影响如图3(b)所示，截止频率随λ增大而减小，图中平滑参数λ分别为1、2、4、10、20、50、300，相应的3 dB截止频率分别是0.189、0.132、0.093、0.059、0.041、0.025 Hz和0.011 Hz(对应1 Hz采样频率)。除了平滑参数λ，频率响应当然还与信号z的采样率有关。因此在给定信号z的情况下，其频率响应的调整仅依赖于一个平滑参数λ，受人为因素影响较小，使用起来非常简单。

对于地震波来说，体波的频率范围约为0.02～100 Hz[17]，一般地震波首波都是高频信号。为了进一步验证，笔者对P744测站E-W方向分量作了功率谱分析，结果如图4所示。从图中可以看出，高频地震波的能量主要集中在0.05～0.40 Hz频段上，在一般理论地震波中的体波频率范围之内，因此本文假设频率大于0.04 Hz的信号是所需要保留的高频地震波信号，则所设计滤波器的截止频率为0.04 Hz。对应频率响应关系，本文选取λ=30。在一般的正则化方法中，正则化参数的选择是难点，而本文根据实际物理问题本身的特点(地震波传播频率固定，与震级大小、震中距远近无关)选取了截止频率，进而确定正则化参数。这也解决了应用正则化方法时平滑因子的选取往往会受到人为因素影响的问题。根据地震波的频率范围，同时结合FFT结果，笔者认为在以后使用SPM方法处理1 Hz的GPS数据时，可选取30作为滤波参数的经验值。

图3[11] L的时变频率响应 Fig.3[11] Time-varying frequency response of L

图4 测站P744在E-W方向动态位移序列功率谱图 Fig.4 The power spectrum of dynamic displacement sequence of station P744 in E-W direction

2 结 果

图5是3次地震前后相应测站高频GPS动态位移时间序列。从图中结果可以明显看出，高频GPS动态位移记录了地震导致的地壳永久变形和地震波，但两种信号是混叠在一起的。在利用GPS资料反演地震参数或地震破裂过程时，通常需要永久性变形信息或地震波信息作为约束条件。如果参数不够准确，势必会导致反演结果存在一定的偏差，因此，在实际应用中，需要将永久性变形信息和高频地震波信息精确分离。

图6至图8分别给出了3次地震前后相应测站高频GPS动态位移的SPM信号分离结果。从图中可以看出，地震引起的永久变形和地震波得到了有效分离。将SPM方法获取的静态永久变形与用地震前后静态天解计算的结果进行了比较，结果见表2和表3。从表中的结果可以看出，SPM分离结果与GPS静态天解结果有一些差异。产生该差异的原因有二。一是在地震发生后，震中距较近的GPS站产生较大的静态永久变形，震后一段时间由于断层余滑的影响，在静态天解的静态永久变形结果中包含了主震与震后余滑的共同效应。因此SPM结果与GPS静态天解结果具有差异。如2010-04-04 EI Mayor-Cucapah地震P496和P744测站南北分量，2012-09-05哥斯达黎加地震QSEC和SAJU测站东西分量及南北分量。二是在主震发生后，强余震也可能会进一步引起震中距较近的GPS站产生较大的静态永久变形。因此，在静态天解结果中包含了主、余震引起的静态永久变形的共同效应，使得SPM结果与GPS静态天解结果具有差异。例如2011-03-11东日本大地震在Mw 9.1级主震发生后约30 min又发生了Mw 7.9级强余震，此后约10 min，Mw 9.1震中东部外海又发生了一次Mw 7.7级强余震。由图8和表3可以看出，Mw 7.9级强余震发生后，分别引起MIZU和USUD有2.3 cm和2.0 cm的东西向永久性变形。由于Mw7.7级强余震距离陆地较远，两个测站都没有观测到明显的静态永久变形。

[11] TARVAINEN M P, RANTA-AHO P O, KARJALAINEN P A. An Advanced Detrending Method with Application to HRV Analysis[J]. IEEE Transactions on Biomedical Engineering, 2002, 49(2): 172-175.

图1 测站分布与震中位置 Fig.1 The locations of the stations and the epicenters

图2 MSF恒星日滤波背景噪声分析(以MIZU测站N方向为例) Fig.2 The background noise analysis of MSF method (taking the north component of MIZU station as an example)

注：(a)、(b)对应2010-04-04 EI Mayor-Cucapah地震；(c)、(d)对应2012-09-05哥斯达黎加地震；(e)、(f)对应2011-03-11东日本大地震 图5 各测站高频GPS动态位移时间序列 Fig.5 Time series of high-rate GPS dynamic displacement of each station

图6 2010-04-04 EI Mayor-Cucapah地震P496和P744测站SPM滤波的结果 Fig.6 The results of SPM filtering at P496 and P744 station

图7 2012-09-05哥斯达黎加地震QSEC和SAJU测站SPM滤波的结果 Fig.7 The results of SPM filtering at QSEC and SAJU station

图8 2011-03-11东日本大地震MIZU和USUD测站SPM信号分离的结果 Fig.8 The results of SPM method at MIZU and USUD station

如图8(b)和图8(d)给出的2011-03-11东日本大地震经过SPM信号分离得到的MIZU和USUD测站高频地震波信号。其中黑色虚线对应的时间为Mw9.1级主震P波到时，紫色虚线对应的时间为Mw 7.9级强余震P波到时。从图中可知，MIZU测站在Mw 9.1级主震发生后约40 s记录到了地震波信号，在Mw 7.9和Mw 7.7级余震发生后并没有记录到地震波信号；USUD测站在Mw 9.1级主震发生后约105 s记录到了地震波信号，在Mw 7.9级余震发生后约50 s记录到了地震波信号，而在Mw7.7级余震发生后没有记录到地震波信号。

表2 不同方法所得静态永久变形及中误差或标准偏差

Tab.2 Static offset obtained by different methods and the corresponding mean square errors cm

测站名称SPM方法解算1Hz静态永久变形GPS静态天解解算永久变形UNAVCO网站结果东西估值/中误差南北估值/中误差东西估值/中误差南北估值/中误差东西/南北P4962．55/0．18-16．42/0．222．44/0．20-18．01/0．102．40/-18．11P7441．10/0．13-6．79/0．171．33/0．11-8．05/0．041．01/-8．08QSEC-1．14/0．14-32．97/0．14-5．70/0．69-39．92/0．59-5．91/-39．94SAJU-31．29/0．23-19．30/0．12-34．61/0．44-26．00/0．53—

表3 不同方法所得静态永久变形及中误差或标准偏差

Tab.3 Static offset obtained by different methods and the corresponding mean square errors cm

测站名称SPM方法解算1Hz静态永久变形Mw9．1主震后Mw7．9余震后GPS静态天解解算永久变形文献[18]东西估值/中误差南北估值/中误差东西估值/中误差南北估值/中误差东西/南北东西/南北MIZU208．32/0．35-114．41/0．25210．63/0．25-115．26/0．21230．06/-123．36204．9/-125．0USUD17．67/0．174．25/0．2819．69/0．234．37/0．2422．45/5．3722．7/5．3

从以上的分析结果可以得出，利用SPM方法，在完整、清楚地获取地震波信息的同时，还可以清晰地分离主震引起的静态永久变形、余震引起的静态永久变形和震后断层余滑量。

3 结 论

式中，H∈Rn×m为观测矩阵；θ∈Rm为回归参数；m是回归参数的个数；ε为观测误差。求解回归参数的估计值时，常用最小二乘准则来估计，但由于静态永久变形信号为非稳态信号，因此，本文引入了正则化估计方法

(2) 相比于小波等其他信号分离的方法，SPM滤波方法最大优点是其频率响应的调整仅依赖于一个平滑参数λ，使用起来非常方便简单。当利用SPM方法分离高频GPS动态位移中的静态永久变形和高频地震波信息时，不必对其做频谱分析，λ=30可作为经验值直接使用。

(3) 本文选取的3次地震均为大地震，从分离效果来看，能够实现快速、准确地分离，与前人结果也比较一致。对于小地震引起的静态永久变形和地震波信息，只要高频地震波信号能量在噪声水平以上，也可采用SPM方法快速分离。

时钟拨回到43年前。1975年初，北京大学教师王选带领团队开始研究“748工程”中的子项目“汉字精密照排系统”，用一系列颠覆性的技术创新和发明，解决了汉字字形信息的数字化储存、信息压缩、复原和输出等重大难题，成功研制第四代汉字激光照排系统。自此，中国印刷技术告别“铅与火”，走向“光与电”。

参考文献：

[1] 刘启元, 陈九辉, 李顺成, 等. 新疆伽师强震群区三维地壳上地幔S波速度结构及其地震成因的探讨[J]. 地球物理学报, 2000, 43(3): 356-365.

LIU Qiyuan, CHEN Jiuhui, LI Shuncheng, et al. Passive Seismic Experiment in Xinjiang-Jiashi Strong Earthquake Region and Discussion on Its Seismic Genesis[J]. Chinese Journal of Geophysics, 2000, 43(3): 356-365.

[2] 陈九辉, 刘启元, 李顺成, 等. 青藏高原东北缘-鄂尔多斯地块地壳上地幔S波速度结构[J]. 地球物理学报, 2005, 48(2): 333-342.

CHEN Jiuhui, LIU Qiyuan, LI Shuncheng, et al. Crust and Upper Mantle S-wave Velocity Structure Across Northeastern Tibetan Plateau and Ordos Block[J]. Chinese Journal of Geophysics, 2005, 48(2): 333-342.

1.3.1 SPM方法

[3] 刘瑞丰, 党京平, 陈培善. 利用速度型数字地震仪记录测定面波震级[J]. 地震地磁观测与研究, 1996, 17(2): 1-4.

LIU Ruifeng, DANG Jingping, CHEN Peishan. Determination or Surface Wave Magnitude by Using Velocity Type Digital Seimograph Recordings[J]. Seismological and Geomagnetic Observation and Research, 1996, 17(2): 1-4.

[18] 杨少敏, 聂兆生, 贾志革, 等. GPS解算的日本Mw 9.0级地震的远场同震地表位移[J]. 武汉大学学报(信息科学版), 2011, 36(11): 1336-1339.

[5] LARSON K M, BODIN P, GOMBERG J. Using 1-Hz GPS Data to Measure Deformations Caused by the Denali Fault Earthquake[J]. Science, 2003, 300(5624): 1421-1424.

[6] JI Chen, LARSON K M, TAN Ying, et al. Slip History of the 2003 San Simeon Earthquake Constrained by Combining 1 Hz GPS, Strong Motion, and Teleseismic Data[J]. Geophysical Research Letters, 2004, 31(17): L17608.

[17] LAY T, WALLACE T C. Modern Global Seismology[M]. San Diego, California: Academic Press, 1995.

[8] 李志才, 张鹏, 金双根, 等. 基于GPS观测数据的汶川地震断层形变反演分析[J]. 测绘学报, 2009, 38(2): 108-113, 119.

LI Zhicai, ZHANG Peng, JIN Shuanggen, et al. Wenchuan Earthquake Deformation Fault Inversion and Analysis Based on GPS Observations[J]. Acta Geodaetica et Cartographica Sinica, 2009, 38(2): 108-113, 119.

本实验选用的移动终端APP软件统一为沪江开心词场，因为根据实验前对学生的课堂调查，大部分同学在平时的学习中经常使用此款软件，经过指导教师和多位其他教师商讨，此款软件设计合理，趣味性强，收词全面，师生应用体验一致较好。此款软件在应用软件下载排行上名列第一，用户体验也是最高，5颗星。

[9] LARSON K M, MIYAZAKI S. Resolving Static Offsets from High-rate GPS Data: the 2003 Tokachi-Oki Earthquake[J]. Earth, Planets and Space, 2008, 60(8): 801-808.

[10] 郭爱智, 王勇, 苏晓庆, 等. 利用小波分解重构提取高频GPS数据静态永久性变形[J]. 武汉大学学报(信息科学版), 2013, 38(10): 1192-1195.

因此，在译文中，为了取得语言，语法与词句的通顺，可以做出适当改变，考虑译文风格，接受者语境，使得接受者达到与原文读者一样的反应，产生归化的作用。当译语原文达不到形式上和意义上的完全对等时，形式要完全让位于意义[3]。

GUO Aizhi, WANG Yong, SU Xiaoqing, et al. Resolving Static Offset from High-rate GPS Data by Wavelet Decomposition-reconstruction Algorithm[J]. Geomatics and Information Science of Wuhan University, 2013, 38(10): 1192-1195.

从图6—图8中可以看到，利用SPM方法不仅能够准确获得主、余震引起的静态永久变形，同时，得到的高频地震波信号(图6(b)、6(d)，图7(b)、7(d)，图8(b)、8(d))比未分离前(图5(a)—(f))更加清晰。

巴哈赛车转向节结构优化设 计 ………………………………………… 王新建，张 蕊，耿 杰，钟海涛（42）

[12] GERSCH W. Smoothness Priors[M]∥BRILLINGER D, CAINES P, GEWEKE J, et al. New Directions in Time Series Analysis. New York: Springer, 1993.

滩地地块最高出水点高程为510.0 m，引水进水口水面高程为528.5 m，地形水头为18.5 m。管道水力损失为9.5 m，输水剩余水头为9 m，大于最不利位置给水栓所需的水头2.0 m，管道满足灌溉供水要求。

[13] KARJALAINEN P A. Regularization and Bayesian Methods for Evoked Potential Estimation[D]. Kuopio: University of Kuopio, 1997.

[14] 詹金刚, 王勇, 史红岭, 等. 应用平滑先验信息方法移除GRACE数据中相关误差[J]. 地球物理学报, 2015, 58(4): 1135-1144.

ZHAN Jingang, WANG Yong, SHI Hongling, et al. Removing Correlative Errors in GRACE Data by the Smoothness Priors Method[J]. Chinese Journal of Geophysics, 2015, 58(4): 1135-1144.

[15] BOCK H, DACH R, JGGI A, et al. High-rate GPS Clock Corrections from CODE: Support of 1 Hz Applications[J]. Journal of Geodesy, 2009, 83(11): 1083-1094.

微信平台因其近乎垄断的市场地位，规模庞大的用户群体，强大的管理、推广、统计功能，受到越来越多企事业单位的青睐。微信平台作为图书馆进行阅读推广活动的重要媒介，有着举足轻重的地位，许多基于微信平台的线上活动竞相开展，同时微信平台还承担着各种线下活动的宣传推广。近年来，基于微信平台的阅读推广研究也逐渐增加，利用中国知网CNKI期刊检索中的高级检索功能，以“微信”和“阅读推广”两个词作为主题词进行检索，检出期刊文章共404篇，其中核心期刊67篇，国内最早开始研究“微信阅读推广”的文章出现于2013年，并呈逐年上升趋势（检索时间：2018年7月17日）。

消费者关于再制造产品的认知程度、环保认同度、价格差异接受度将直接影响其对再制造产品的支付意愿，从而引起市场需求的变化。对此，主要考虑以下问题：认知差异将如何影响再制造产品需求以及基于三者差异，支付意愿如何改变，对需求有何影响。根据研究需要，做如下假设：

[16] CHOI K, BILICH A, LARSON K M, et al. Modified Sidereal Filtering: Implications for High-rate GPS Positioning[J]. Geophysical Research Letters, 2004, 31(22): L22608.

[7] ANZIDEI M, BOSCHI E, CANNELLI V, et al. Coseismic Deformation of the Destructive April 6, 2009 L’Aquila Earthquake (Central Italy) from GPS Data[J]. Geophysical Research Letters, 2009, 36(17): L17307.

[4] LARSON K M. GPS Seismology[J]. Journal of Geodesy, 2009, 83(3-4): 227-233.

YANG Shaomin, NIE Zhaosheng, JIA Zhige, et al. Far-field Coseismic Surface Displacement Caused by the Mw 9.0 Tohoku Earthquake[J]. Geomatics and Information Science of Wuhan University, 2011, 36(11): 1336-1339.

[19] 许才军, 刘洋, 温扬茂. 利用GPS资料反演汶川Mw 7.9级地震滑动分布[J]. 测绘学报, 2009, 38(3): 195-201, 215.

XU Caijun, LIU Yang, WEN Yangmao. Mw 7.9 Wenchuan Earthquake Slip Distribution Inversion from GPS Measurements[J]. Acta Geodaetica et Cartographica Sinica, 2009, 38(3): 195-201, 215.

[20] 吴继忠, 吴文坛. 利用高频GPS进行地震动态变形分析及地震定位[J]. 大地测量与地球动力学, 2012, 32(2): 20-23.

WU Jizhong, WU Wentan. Kinematic Analysis of Coseismic Deformation and Seismic Location Using High-rate GPS[J]. Journal of Geodesy and Geodynamics, 2012, 32(2): 20-23.

[21] JIANG Zaisen, WANG Min, WANG Yanzhao, et al. GPS Constrained Coseismic Source and Slip Distribution of the 2013 Mw 6.6 Lushan, China, Earthquake and Its Tectonic Implications[J]. Geophysical Research Letters, 2014, 41(2): 407-413.

[22] 谭凯, 乔学军, 杨少敏, 等. 汶川地震GPS形变约束的破裂分段特征及滑移[J]. 测绘学报, 2011, 40(6): 703-709.

TAN Kai, QIAO Xuejun, YANG Shaomin, et al. Rupture Characteristic and Slip Constrained by GPS Coseismic Deformation Induced by the Wenchuan Earthquake[J]. Acta Geodaetica et Cartographica Sinica, 2011, 40(6): 703-709.

[23] WANG Qi, ZHANG Peizhen, FREYMUELLER J T, et al. Present-day Crustal Deformation in China Constrained by Global Positioning System Measurements[J]. Science, 2001, 294(5542): 574-577.

[24] 李萌, 黄丁发, 严丽, 等. 汶川地震前后四川盆地CORS站运动特性分析[J]. 测绘学报, 2014, 43(6): 582-589. DOI: 10.13485/j.cnki.11-2089.2014.0107.

LI Meng, HUANG Dingfa, YAN Li, et al. Characteristics of Position Time Series at CORS Stations in Sichuan Basin before and after Wenchuan Earthquake[J]. Acta Geodaetica et Cartographica Sinica, 2014, 43(6): 582-589. DOI: 10.13485/j.cnki.11-2089.2014.0107.

[25] 王敏, 李强, 王凡, 等. 全球定位系统测定的2011年日本宫城Mw 9.0级地震远场同震位移[J]. 科学通报, 2011, 56(20): 1593-1596.

WANG Min, LI Qiang, WANG Fan, et al. Far-field Coseismic Displacements Associated with the 2011 Tohoku-Oki Earthquake in Japan Observed by Global Positioning System[J]. Chinese Science Bulletin, 2011, 56(20): 1593-1596.

[26] 徐韶光, 熊永良, 廖华, 等. 汶川地震同震形变的静态和动态分析[J]. 大地测量与地球动力学, 2010, 30(3): 27-30, 54.

XU Shaoguang, XIONG Yongliang, LIAO Hua, et al. Static and Kinematic Analysis of Coseismic Deformation of Wenchuan Ms 8.0 Earthquake[J]. Journal of Geodesy and Geodynamics, 2010, 30(3): 27-30, 54.