卫星重力测量技术[1]日趋成熟，如CHAMP、GRACE以及GOCE卫星[2-5]。卫星重力代表的是地球重力场的中长波部分[6]，即地球深部信息。进行重力场与其场源深度关系的研究，能更好地推进卫星重力数据的实际应用。

通过研究小波基和小波阶数，可以对重力场进行小波多尺度分析[7-11]，再结合功率谱埋深估计技术，就能对重力场与场源深度的关系作出推断。本文使用恰当的小波基和小波阶数分解卫星重力场，并根据其结果计算南海北部陆缘地区的莫霍面深度，结果表现良好。

潮州麦秆画作品中所蕴含的历史文化内涵，是潮州麦秆画区别于其它地区麦秆画作品的重要因素，潮州麦秆画精良的制作技艺与深刻的文化内涵，使潮州麦秆画有别具一格的艺术表现力与艺术感染力。在体现历史文化内涵方面，风景题材的麦秆画作品是其中的翘楚，与传统的花鸟麦秆画作品不同，风景题材的作品更具地域特色，通过对潮州特定景观进行取景创作而成的潮州麦秆画作品，具有独特的潮州历史文化内涵。

1 重力场场源模型设置

卫星重力是不同深度物质引起的叠加重力场。因此，在设置场源时，模型应设置为多层，且场源必须包括地球表面和深部物质。本文以两层分布的5个质量块体来模拟地球局部重力场，也符合多尺度重力场的设定。

1.1 模型重力场原理

采用常用的直立六面体组合作为重力场场源模型，其计算公式为[12]：

(1)

使用最优小波分解参数(小波基db19，分解阶数4)计算南海北部陆缘地区布格重力异常的4阶小波逼近(图15)。布格重力异常一般表征了莫霍面的起伏形态，其小波分解的低频场(即区域布格重力异常)与莫霍面起伏更为吻合[20]。根据功率谱波数图(图16)，南海北部陆缘地区莫霍面深度约为28.3 km，这与北部陆架地区莫霍面深度在26～30 km之间的结论[20]相吻合。

基因检测结果分析:MEN1基因是多发内分泌腺肿瘤1型(Multiple endocrine neoplasia 1)的致病基因,为常染色体显性遗传方式(AD)。对于该类遗传方式,杂合变异可能导致发病。在受检者 MEN1基因发现c.24_31delGACGCTGT(编码区第24_31号核GACGCTGT缺失)的杂合核苷酸变异,该变异导致从第8号氨基酸开始的氨基酸合成发生改(p.8fs),为移码变异。该变异可能导致蛋白质功能受到影响。该变异的致病性尚未见文献报道。该变异不属于多态性变化,在人群中发生的频率极低。

1.2 模型参数

不同阶的小波低频重力场对应不同的地质体深度信息。根据位场频谱理论[17-19]，其场源埋深为：

图1 模型示意图 Fig.1 The sketch figure of model

表1 组合模型参数 Tab.1 The parameter of the combination model

模型剩余密度/g·cm-3ξ1/kmξ2/kmη1/kmη2/kmζ1/kmζ2/kmM0．2356535652040Ⅰ0．248524852Ⅱ0．248521317Ⅲ0．213178387Ⅳ0．2838783870．52．5

1.3 模型结果

本文研究区域为南海北部陆缘地区[20]，范围115°～120°E、18°～23°N。重力数据为分辨率1′×1′的布格重力异常(图14)，来自斯克里普斯海洋协会和美国国家海洋与大气局卫星测高实验室共同维护的全球卫星重力异常数据库。

2 小波多尺度分解算法

重力场是叠加场，为研究不同规模、形态和埋深的场源物质，必须对重力异常进行分离[13-15]，主要方法有滑动平均算法、导数计算法、小波算法等。小波多尺度分析是近年来研究热点，但小波基和分解阶数的选择仍然是一个难点[16]。

图2 模型重力异常 Fig.2 The map of model gravity anomaly

本文利用Biorthogonal、Coiflets、Daubechies、Discrete Meyer、Reverse Biorthogonal、Symlets 6种小波基对总场重力异常(图2(a))进行分解，阶数选取为1～20阶，并求取分解结果的高频场与局域重力异常的相关系数。图3～8只列出了差异较大的结果。

图3 Biorthogonal小波基系相关分析图 Fig.3 The correlation figure of Biorthogonal wavelets

图4 Coiflets小波基系相关分析图 Fig.4 The correlation figure of Coiflets wavelets

图5 Daubechies小波基系相关分析图 Fig.5 The correlation figure of Daubechies wavelets

图6 Discrete Meyer小波基系相关分析图 Fig.6 The correlation figure of Discrete Meyer wavelets

图7 Reverse Biorthogonal小波基系相关分析图 Fig.7 The correlation figure of Reverse Biorthogonal wavelets

图8 Symlets小波基系相关分析图 Fig.8 The correlation figure of Symlets wavelets

在每种小波基系中选取相关系数最接近1的小波基进行对比分析(图9)。由图9可知，小波分解阶数为4时，所有小波基的计算结果最接近模型局域重力异常，且其中coif3、db19、rbio2.8、sym7的计算结果极为相近。通过实际测试，最终选定db19为最优小波基，最优小波分解阶数为4。

2)根据功率谱波数图估算出的深度具有统计意义，代表区域平均深度。

图9 最优小波基系相关分析图 Fig.9 The correlation figure of optimization wavelets

图10 小波分解低频结果图 Fig.10 The map of wavelet low frequency

图11 小波分解高频结果图 Fig.11 The map of wavelet high frequency

3 深度估计算法

本文组合模型由1个大长方体和4个小长方体构成(图1)，区域重力异常由埋藏较深的大长方体产生，局域重力异常由埋藏较浅的小长方体产生。该组合具有复杂的排列分布(左右对称、中心不对称)，且有大异常、小异常，以此模拟真实重力场的复杂程度。组合模型参数见表1，M为大长方体(平均埋深30 km)；Ⅰ、Ⅱ、Ⅲ、Ⅳ为小长方体(平均埋深1.5 km)；剩余密度为0.2 g/cm3。

1.2.2.3 以方案的多版本呈现模式提高目视提醒效应以保证皮肤护理方案执行到位的可行性 将可行性皮肤护理方案制作成多个版本，护士随身携带便携式卡片版本，电子版本置于护士站工作电脑桌面滚动播出，B2纸质版本张贴于病区走廊目视可见处，通过上述呈现方式时刻警示与指引护理人员规范全面地落实皮肤护理方案。

(2)

式中，与Δkj是与kj的变化率，j=1,2,…,J为分解阶数，hj为埋深，kj为波数，为功率谱。

根据功率谱波数图(图12、13)，由式(2)计算低频和高频场的平均埋深分别为29.4 km、1.9 km，与预设模型的平均埋深基本一致。

图12 低频功率谱波数图 Fig.12 The power spectrum figure of low frequency

图13 高频功率谱波数图 Fig.13 The power spectrum figure of high frequency

4 卫星重力多尺度分析

根据表1中的参数，按照式(1)计算模型重力场(图2)，作为研究输入数据。

图14 研究区布格重力异常图 Fig.14 Bouguer gravity anomaly in study area

式中，P为牛顿万有引力常数；m、n为观测面横纵向点数；t为直立六面体个数；ξ、η、ζ为直立六面体的8个角坐标(如图1)为待计算点坐标。

图15 研究区4阶小波逼近图 Fig.15 4 orders wavelet approximations in study area

图16 研究区4阶小波逼近功率谱波数图 Fig.16 The power spectrum figure of 4 orders wavelet approximations in study area

为进一步研究功率谱深度估计算法的准确性，从南海北部陆缘地区裁剪出珠江口盆地和台西南盆地的4阶小波逼近布格重力异常(图17、18)，分别进行深度计算，结果为23.9 km、20.8 km(图19、20)。

图17 珠江口4阶小波逼近图 Fig.17 4 orders wavelet approximations in Zhujiangkou

图18 台西南4阶小波逼近图 Fig.18 4 orders wavelet approximations in Taixinan

图19 珠江口4阶小波逼近功率谱波数图 Fig.19 The power spectrum figure of 4 orders wavelet approximations in Zhujiangkou

图20 台西南4阶小波逼近功率谱波数图 Fig.20 The power spectrum figure of 4 orders wavelet approximations in Taixinan

1)布格重力异常区域幅值趋势与莫霍面深度分布相关，使用小波多尺度分解提取布格重力异常区域信息很有必要。

表2 莫霍面深度结果 Tab.2 The depth of the Moho

构造名称布格重力异常/mGal文献[20]/km计算结果/km南海北部陆缘-65．4～286．226～3028．3珠江口盆地-20～11023～2523．9台西南盆地-0．3～19017～2220．8

5 结 语

将以上深度计算结果与收集的资料进行对比(表2)，功率谱波数图(图16、19、20)得到的深度与文献[20]结论一致，表明功率谱的深度计算结果可以表征区域莫霍面平均深度。

使用最优小波分解参数对总场重力异常进行处理，其结果的高频和低频部分对应于局域、区域模型重力异常，幅值大小和异常形态与输入重力场极为接近(图10、11)。

3)深度估算结果表明，研究区布格重力异常的区域背景主要是由莫霍面产生的。

Hello Time 2 sec Max Age 20 sec Forward Delay 15 sec

742例患者共置入2 968枚螺钉，0级2 520枚，1级375枚，2级54枚，3级19枚，置钉优良率为97.5%。19枚3级置钉来源于17例患者，均表现为神经根刺激症状，根据神经根症状程度酌情使用脱水剂和地塞米松3 d，指导患者逐步行直腿抬高锻炼。10例患者随访观察症状均逐步缓解；7例患者症状明显，严重影响生活质量，行开放翻修手术调整螺钉位置，末次随访时恢复良好。

参考文献

[1] 曾华霖.重力场与重力勘探[M].北京：地质出版社，2005(Zeng Hualin. Gravity Field and Gravity Exploration[M]. Beijing: Geological Publishing House, 2005)

[2] 罗志才，吴云龙，钟波，等.GOCE卫星重力梯度测量数据的预处理[J].武汉大学学报：信息科学版，2009, 34(10)：1 163-1 167(Luo Zhicai，Wu Yunlong，Zhong Bo，et al. Preprocessing of the GOCE Satelite Gravity Gradiometry Data[J]. Geomatics and Information Science of Wuhan University, 2009, 34(10):1 163-1 167)

[3] 祝竺，周泽兵.用于卫星重力梯度测量的加速度计性能参数分析[J].大地测量与地球动力学，2012, 32(5)：148-153(Zhu Zhu，Zhou Zebing. Analysis of Accelerometer Performance Parameters for Satellite Gravity Gradiometry[J]. Journal of Geodesy and Geodynamics, 2012, 32(5)：148-153)

[4] 詹金刚，王勇，郝晓光.GRACE时变重力位系数误差的改进去相关算法[J].测绘学报，2011, 40(4)：442-446(Zhan Jingang，Wang Yong，Hao Xiaoguang. Improved Method for Removal of Correlated Errors in GRACE Data[J]. Acta Geodaetica et Cartographica Sinica，2011, 40(4)：442-446)

[5] 周新，邢乐林，邹正波，等.GRACE时变重力场的高斯平滑研究[J].大地测量与地球动力学，2008, 28(3)：41-45(Zhou Xin，Xing Lelin，Zou Zhengbo，et al. Study on Gaussian Smoothing of GRACE Temporal Gravity Variation[J]. Journal of Geodesy and Geodynamics, 2008, 28(3)：41-45)

[6] 王留朋，郭燕平，冯炜.GOCE重力卫星在军事上的应用前景分析[J].地理空间信息，2011, 9(1)：11-15(Wang Liupeng，Guo Yanping，Feng Wei. GOCE Satellite Gravity in Military Application Prospect Analysis[J]. Geospatial Information, 2011, 9(1)：11-15)

[7] 陈国雄，孙劲松，刘天佑.GRACE卫星时变重力场的小波多尺度分解——以2008年汶川Ms8.0大地震为例[J].武汉大学学报：信息科学版，2012, 37(6)：679-682(Chen Guoxiong，Sun Jinsong，Liu Tianyou. Wavelet Multi-Scal Decomposition of Time Variable Gravity Field Detected by GRACE Satellite: A Case from Wenchuan Ms8.0 Earthquake,2008[J]. Geomatics and Information Science of Wuhan University, 2012, 37(6)：679-682)

[8] 高德章，侯遵泽，唐建.东海及邻区重力异常小波多尺度分解[J].地球物理学报，2000, 40(6)：842-849(Gao Dezhang，Hou Zunze，Tang Jian. Multiscale Analysis of Gravity Anomalies on East China Sea Adjacent Regions[J]. Chinese Journal of Geophysics, 2000, 40(6)：842-849)

[9] 杨文采，施志群，侯遵泽.离散小波变化与重力异常多重分解[J].地球物理学报，2001, 44(4)：534-541(Yang Wencai，Shi Zhiqun，Hou Zunze. Discrete Wavelet Transform for Multiple Decomposition of Gravity Anomalies[J]. Chinese Journal of Geophysics, 2001, 44(4)：534-541)

[10] 秦前清，杨宗凯. 实用小波分析[M]. 西安：西安电子科技大学出版社， 1994(Qin Qianqing，Yang Zongkai. Practical Wavelet Analysis[M]. Xi’an: Xidian University Press, 1994)

[11] 宁津生，汪海洪，罗志才.小波分析在大地测量中的应用及其进展[J].武汉大学学报:信息科学版，2004, 29(8)：659-664(Ning Jinsheng，Wang Haihong，Luo Zhicai. Applications of Wavelet Analysis in Geodesy and Its Progress[J]. Geomatics and Information Science of Wuhan University, 2004, 29(8)：659-664)

[12] 重磁资料数据处理问题编写组.重磁资料数据处理问题[M]. 北京：地质出版社，1997(Editorial and Magnetic Data Processing Problem. Magnetic Data Processing Problem[M]. Beijing: Geological Publishing House, 1997)

[13] 赵勇.小波分析在重力异常分离中的应用研究[D]. 上海：同济大学，2003(Zhao Yong. The Application of Wavelet Analysis in Gravity Anomaly Separation Research[D]. Shanghai: Tongji University, 2003)

[14] 许闯. 多尺度重力反演方法及其在城市活动断层探测中的应用研究[D]. 武汉：武汉大学，2014(Xu Chuang. Study on Multi-Scale Gravity Inverse Method and Its Application to Detect Urban Active Faults[D]. Wuhan: Wuhan University, 2014)

[15] 侯遵泽，杨文采. 中国重力异常的小波变换与多尺度分析[J].地球物理学报， 1997, 40(1)：85-95(Hou Zunze，Yang Wencai. Wavelet Transform and Multi-Scale Analysis on Gravity Anomalies of China[J]. Chinese Journal of Geophysics, 1997, 40(1)：85-95)

[16] 李健，周云轩，许惠平. 重力场数据处理中小波母函数的选择[J]. 物探与化探， 2001，25(6): 410-417(Li Jian，Zhou Yunxuan，Xu Huiping. The Selection of Wavelet Generating Functions in Data-Processing of Gravity Field[J]. Geophysical and Geochemical Exploration, 2001, 25(6): 410-417)

[17] Spector A, Grant F S. Statistical Models for Interpreting Aeromagnetic Data[J]. Geophysics, 1970, 35(2): 293-302

[18] Syberg F J R. A Fourier Method for the Regional-Residual Problem of Potential FieIds[J].Geophysical Prospecting, 1972, 20(1): 47-75

[19] 邱宁，何展翔，昌彦君. 分析研究基于小波分析与谱分析提高重力异常的分辨能力[J]. 地球物理学进展，2007，22(1): 210-214(Qiu Ning，He Zhanxiang，Chang Yanjun. Ability of Improving Gravity Anomaly Resolution Based on Multiresolution Wavelet Analysis and Power Spectrum Analysis[J]. Progress in Geophysics, 2007，22(1): 210-214)

[20] 秦静欣，郝天珧，徐亚，等．南海及邻区莫霍面深度分布特征及其与各构造单元的关系[J]．地球物理学报，2011，54(12): 3 171-3 183(Qin Jingxin，Hao Tianyao，Xu Ya，et al. The Distribution Characteristics and the Relationship between the Tectonic Units of the Moho Depth in South China Sea and Adjacent Areas[J]. Chinese Journal of Geophysics, 2011, 54(12): 3 171-3 183)

在贫困地区培育新型职业农民的过程中，不仅需要落实相关的扶持政策，还应建立高效的衔接机制，有效实现贫困农民培训需求、培训内容与就业创业需求之间的对接。要将新型职业农民培育工作与土地流转、产业扶持、人才奖励激励、金融保险等多种扶持政策相衔接。要推动农技推广、农业科研机构和其他专家开展面向新型职业农民的后续服务，对接创业扶持项目，提升职业农民就业创业、发家致富的能力，切实将贫困地区的新型职业农民培育工作与精准扶贫战略有机结合起来，为贫困人口尽快实现精准脱贫提供有力保障。