DreamBlooms/HansPub-15000-to-45000 huggingface dataset:
system: [SYSTEM]

user: 2. 库鲁斯台草原概况 2.1. 地理概况 如图1，库鲁斯台草原位于新疆塔额盆地腹地，地处塔尔巴哈台山、乌尔喀沙山和巴尔鲁克山三座断块山地包围的陷落盆地腹部中心地带，北以塔额公路为界，南以托里–裕民公路为界，西至中国与哈萨克斯坦国界，东接额敏县杰勒阿尕什乡沿加依尔山延伸至托里县加尔巴斯牧场。草原地势平坦，海拔400~600 m，额敏河贯穿中央，东西长76 km，南北宽36 km，土地总面积25.93万m 2 ，可利用草场面积21万hm 2 。库鲁斯台草原是全国第二大连片平原草场，也是新疆绿洲经济的重要组成部分，在塔城盆地及周边发挥着“肾脏功能”作用，生态价值远远超过经济价值。重视库鲁斯台草原区域生态环境的保护和治理，不仅对塔城盆地及周边地区生态文明建设有着极为重要的意义，而且对于整个新疆社会经济稳定发展具有十分重要的影响。 图1. 库鲁斯台地理位置图 2.2. 库鲁斯台草原1961~2016年温度和降水变化 以塔城、额敏、裕民、托里四站1961~2016年气象资料(四站平均)作为表征库鲁斯台草原气候的资料。1961~2016年：年平均气温6.7℃；年平均累积降水279 mm。 由图2可以看出：库鲁斯台草原近56a气温增长了2.37℃，平均升温速率为0.474℃·(10a) −1 。 滑动平均库鲁斯台草原1961~2016年降水量得到图3：60年代至80年代库鲁斯台年降水呈下降趋势，20世纪80年代中期至低点后持续呈上升趋势，存在不规律波动。 图2. 库鲁斯台草原1961~2016年湿度变化 图3. 库鲁斯台草原1961~2016年降水量年变化趋势
system: [SYSTEM]

user: 3. 资料与方法 3.1. 资料 本文使用由新疆信息中心质控过库鲁斯台周边四个国家站点(塔城、额敏、托里、裕民) 1961~2016年温度及降水资料，该资料是无缺值连续资料。由图1可见库鲁斯台草原位于新疆塔额盆地腹地(E82˚36'~84˚00， N46˚09'~46˚44')，以四站平均资料作为库鲁斯台草原气象资料。构建库鲁斯台草原1961~2016年12个月尺度SPI和SPEI序列。 3.2. 方法 采用具有多时间尺度特征的标准化降水指数(SPI)与在SPI基础上考虑蒸散影响的标准化降水蒸散指数(SPEI)作为气候干湿状况的表征。不同时间尺度干旱指数的第一个月对应的是当前月，对n个月的时间尺度，从时间上自当月向前延续n−1个月。SPI指数是表征某时段降水量出现概率多少的指标之一，该指标适合月尺度以上相对于当地气候状况的干旱监测与评估。 SPI指数的构建过程如下： 假设某一时段的降水量为x，则其Γ分布的概率密度函数为： f ( x ) = 1 β λ Γ ( x ) x γ − 1 l − x l β       ( x > 0 ) (1) 式中：β为尺度参数；γ为形状参数；x为降水量；β和γ可用极大似然估计方法求得： λ ^ = 1 + 1 + 4 A / 3 4 A (2) β ^ = x ¯ / y ^ (3) 式中： A = lg x ¯ − 1 n ∑ i = 1 n lg x i (4) 式中：i x 为降水资料样本， x ¯ 降水量气候平均值，n为计算序列的长度。确定概率密度函数中的参数后，对于某一年的降水量 x 0 ，可求出随机变量x小于 x 0 事件的概率为： F ( x < x 0 ) ∫ 0 ∞ f ( x ) d x (5) 利用数值积分可以计算用(1)和(5)式后的事件概率近似估计值。降水量为0时的事件概率为： F ( x = 0 ) = m / n (6) 式中：m为降水量为0的样本数；n为总样本数。对Γ分布概率进行正太标准化处理，即将式(5)、(6)求得的概率值代入标准化正态分布函数，即： F ( x < x 0 ) = 1 2 π ∫ 0 ∞ e − z 2 / 2 d x (7) 对式(7)进行近似求解可得： S P I = S t − ( c 2 t + c 1 ) t + c 0 ( ( d 3 t + d 2 ) t + d 1 ) t + 1.0 (8) 式中： t = ln 1 F 2 ，当 F > 0. 5 时，   S = 1 ； F ≤ 0.5 时， S = − 1 。 c 0 = 2.515517 ， c 1 = 0.802853 ， c 2 = 0.010328 ； d 1 = 1.432788 ， d 2 = 0.189269 ， d 3 = 0.001308 。 SPEI指数的计算与SPI指数类似，计算时首先以某一时段降水量与潜在蒸散量的差值替代降水量，然后进行正态标准化，潜在蒸散采用Penman-Monteith公式(Allen等，1998)计算。 根据国家气象干旱等级标准，划分为5个等级，并确定相应的SPEI/SPI界限值，各级干旱指标见表1。 Table 1 Drought classification of SPEI and SP 等级类型 SPEI和SPI 1 无旱 (−0.5， 0) 2 轻旱 (−1， −0.5] 3 中旱 (−1.5， −1] 4 重旱 (−2， −1.5] 5 特旱 (−∞， −2] 表1. SPEI和SPI干旱等级划分 SPI和SPEI是一种标准化的干旱指数，负值表示有干旱事件发生，正值则表示气候相对湿润，绝对值越大，表明干旱或湿润程度也越严重。SPEI具有多时间尺度的优点，1~6个月时间尺度的SPEI反映的是地表土壤水分变化，可以表征干旱对农业的影响。而12个月时间尺度SPEI与水文变化密切相关，反映的是某一地区长期的水分盈缺状态。相关研究指出，12个月及以上尺度的SPEI在新疆地区有较好的适用性，因此，本文以12个月时间尺度SPEI作为研究旱变化的指标。
system: [SYSTEM]

user: 4. 结果分析 本文采用SPI和SPEI作为干旱指标，对库鲁斯台草原干旱演变特征进行对比分析，得到同一时间段上两种指数统计出的干旱事件发生次数及发生年份基本统一，但同一降水变化序列，两种指数在时间变化上存有明显差异。将两种指数时空变化趋势与降水量和平均气温的时空变化趋势结果进行对比，以考察温度变化对库鲁斯台草原干湿变化的影响。 图4. 库鲁斯台草原1961~2016年12个月尺度的SPI 图5. 库鲁斯台草原1961~2016年12个月尺度的SPEI 计算库鲁斯台草原56a 12个月尺度的SPI和SPEI，并进行线性回归平滑处理，得到其时间变化趋势图(图4、图5)。由图4可知，近56a库鲁斯台草原SPI整体呈现波动上升的趋势，线性趋势拟合也表明这种上升趋势是极其显著的，而SPEI整体呈波动下降的趋势。两种指数在以12个月为尺度的时间趋变化上均存在着较大差异，SPEI呈变干趋势，SPI则变湿趋势，与SPI相比，SPEI考虑了热量因子对潜在蒸散发的贡献，能较好地对气候变暖背景下研究区干旱进行刻画。 SPEI和SPEI序列曲线并分别于1962、1965、1974、1997、2000年及2008年下降到−1以下区域，意味着这6个年份发生了较为严重的干旱事件。文献 [ 5 ] [ 6 ] 及气象报道同样指出，新疆于1997、2006年和2008年发生普遍干旱灾害，特别是2008年被认为是有气象记录以来第二个严重干旱年。SPEI指数很好的探测到这些干旱事件的发生，充分显示了12个月尺度SPEI可以有效地反映库鲁斯台草原地区干旱发生及持续时间
system: [SYSTEM]

user: 5. 结论 1) 基于降水和蒸散的SPEI可以更灵活的反映月季干旱变化特征，而SPI指数没有温度变化对蒸散的影响，不利于表征月季干旱变化特征。 2) 近56a中，两种指数判定的同一时间尺度下的干旱事件次数和发生时间基本统一。 3) 两种指数在以12个月为尺度的时间趋变化上均存在着差异，SPEI呈变干趋势，SPI呈变湿趋势，与SPI相比，SPEI考虑了热量因子对潜在蒸散发的贡献，能较好地对气候变暖背景下研究区干旱进行刻画。
system: [SYSTEM]

user: 摘要 通过对大学生预防公共突发事件各个因素分析，结合个人、学校、社会三个方面，选取18个指标建立公共突发事件评价指标模型，利用层次分析法求解，计算出各个指标对总目标的权重，且进行排序。从建模结果可以看出：权重值排前四的分别是公共交通、社会服务、社会治安和出行不便。在预防公共突发事件时，根据建模的结果可以更有目的地实施对策，这对大学生有效预防公共突发事件具有重要参考意义。
system: [SYSTEM]

user: 1. 引言 近年来，随着我国经济快速发展，经贸活动、人员流动频繁，我国的社会问题日益凸显，且个别问题容易被放大，损失了政府的公信力，公共突发事件发生的频率大增。公共突发事件具有社会影响剧烈、爆发频率高、波及范围广等特点，给个人、学校和社会带来不良的影响。大学生预防公共突发事件是目前社会公共突发事件研究热点之一，国内外学者一直关注该领域的研究动态。高校在应对公共突发事件时行政干预较多，缺少完备的法律法规作为约束规范。高校的法律法规制度不够完善，容易遭受社会各界的质疑，有些规章制度甚至出现相互抵触的现象 [ 1 ]。为了提高大学生应对公共突发事件的能力，提供合理可行的应对措施，笔者通过查找文献、结合实际情况从个人、学校、社会三个方面，选取18个指标，建立大学生预防公共突发事件指标体系模型，利用层次分析法求解，计算各个指标的总目标权重值，并进行排序，从而确定各个指标对研究大学生预防公共突发事件的重要性，进行定性与定量研究，根据建模结果提出一些合理的建议。
system: [SYSTEM]

user: 3. 公共突发事件的结论及应用对策 3.1. 公共突发事件研究的结论 公共突发事件的影响指标建模表明：指标权重排名第一的是公共交通，广义的公共交通是指公共运输，其包括民航、铁路、公路、水运等交通方式；狭义的公共交通是指城市范围内定线运营的公共汽车及轨道交通、公共自行车、渡轮、索道等交通方式。在今年新冠疫情期间，国家采取严格管控交通运输 Table 5 The total order of each indicator elemen 三级指标(C i ) B 1 0.1884 B 2 0.0810 B 3 0.7306 相对一级指标 各个权重值 层次总排序 C 1 0.0902 0.0170 12 C 2 0.2215 0.0417 6 C 3 0.0592 0.0112 14 C 4 0.1379 0.0260 9 C 5 0.0368 0.0069 16 C 6 0.4544 0.0856 4 C 7 0.4148 0.0336 8 C 8 0.1532 0.0124 13 C 9 0.0896 0.0073 15 C 10 0.0336 0.0027 18 C 11 0.0540 0.0044 17 C 12 0.2548 0.0206 11 C 13 0.4148 0.3031 1 C 14 0.1532 0.1119 3 C 15 0.2547 0.1861 2 C 16 0.0896 0.0655 5 C 17 0.0541 0.0395 7 C 18 0.0336 0.0245 10 表5. 各个指标元素的总排序 方面的措施，武汉封城，铁路、高速路临时关闭，禁止通行，以防人员流动造成大面积传染。经济水平较高的城市通常会因公共交通管控影响较大，而无法及时采取交通管控措施，会出现不及时抑制疫情传播而导致确诊人数显著增加的结果 [ 4 ]。为了更好的防控疫情同时也要有必要的物资运输，国家也出台了很多关于货车的通行方案。铁路部门主动与各级地方政府、企业联系，根据物资运输需求制定运输组织方案，协调帮助重点企业尽快恢复正常生产。民航部门利用通用航空和全货机运输物资，既保障了紧急物资的及时配送，又最大限度地减少了接触感染 [ 5 ]。指标权重排第二的是社会服务，目前对公共突发事件的对策也符合社会服务的数学建模结果。根据数学建模结果的指标排序以及当前实际，以后应对公共突发事件可根据数学模型结果权重前五的指标进行制定与修改，另外目前属于数字化、智能化时代，在做好前面这几项工作后，应加强完善机械化，做到机械智能稳定有效，精益求精，这对公共突发事件的防控的提升很有帮助。指标权重排第三的是社会治安，在这次的新冠疫情时，设立许多的卫生检查点，疫情检查点，建设疫情防控医院投入大量的医护技术人员参与防控疫情；在应急法制建设上，完善相应的应急法律法规，使各个单位机构能有相应的法律地位，确保各项防控工作有序进行。在处理社会治安方面，把握好处理应急措施与常规问题的尺度，确保最大限度保障人民的健康和财产；对可能引发疫情再度扩散的事情，做好充分预案准备，随时将损失降到最低，尽全力去维护社会稳定。处理疫情引起的社会治安案件时，杜绝笼统处理，要分类处理；做好预防和监控滥用职权强制措施的现象，避免矛盾的激化，保证社会的安定。 在预防公共突发事件中，做好预防公共突发事件是大学生体现为祖国担当的一种责任。为了更好地预防公共突发事件，大学生应当做好本分的事情，学会科学预防公共突发事件知识、运用科学预防公共突发事件知识，也可以宣传预防公共突发事件的知识，从而维护高校校园的稳定，也有利于提高我国公民应对公共突发事件的能力，提高我国预防公共突发事件的管理水平。 3.2. 公共突发事件应对策略 (1) 公共交通在指标权重中排第1。通过地方交通方式以及时间，选择便利的交通出行，且做到预防公共突发所带来的一些影响，比如人与人之间保持安全距离，戴口罩等。目前国家注重公共交通，在各公共交通场所标有提示语，有些还有体温检测、核酸检测等。在高校里，大学生错峰上下课，错峰进入公共场所，按时打卡和上报当日的身体健康状况等。 (2) 社会服务在指标权重中排第2。社会服务在日常生活中必不可少的。在预防公共突发事件时，需要建立一些网络机制，做好可以解决问题的准备。例如，面对疫情危机，公共图书馆要发挥专业优势，多途径开展面向社会公众的信息素养教育等专业活动，用文化和知识的力量减少社会公众的恐慌情绪，助力全社会疫情防控 [ 6 ]。 (3) 社会治安在指标权重中排第3。疫情下的社会治安状况是一个动态的过程，随着疫情的变化而发生改变，公安机关应对这种变化，必须构建一套系统完备、科学有效的应对机制，确保整个应对过程精准有力 [ 7 ]。社会治安是我国社会稳定的基础，同时也是人民人身安全的保障。社会治安是有效预防公共突发事件的重要保证。在高校里，必须建立预防公共突发事件的治安制度，有效地在高校里实施。大学生应当严格遵守学校的规章制度。
system: [SYSTEM]

user: 4. 结语 在处理公共突发事件中，目前还没有一套完全不变的应急预案，也没有绝对的应对措施。在应对过程中有着诸多不定性因素的影响，应急措施还需进一步完善，从模糊数学出发，建立公共突发事件模糊综合评价，尽量将不定性因素的影响降到最低。应急预案不能全面地应对公共突发事件，只有在不断改进与完善才能更好应对。在实践中找到缺陷，从而进行更改，不停地反思与改进。
system: [SYSTEM]

user: 摘要 1 h变温灾害是影响杭州萧山南北白对虾生长发育的重要气象灾害，本文采用杭州萧山区21个区域自动站2020年1~6月的逐日逐时平均气温资料以及两个典型区域自动站：对虾基地站(白对虾集中养殖区)和十七工段站(建站时间最长的站点)从建站至2020年6月的逐日逐时平均气温资料，利用Excel数据统计和Python、Origin、ArcGis进行数据分析，从月、日两个时间尺度分析1 h变温灾害的时间变化规律和在萧山地区不同区域的空间分布特点，探讨其在不同区域的差异性并对其变化规律进行分析。得到以下几点主要结论：1) 2020年2~6月中，变温轻度灾害在4月受灾程度最重，6月受灾程度最轻。十七工段站在2~6月的受灾程度均为最重；外二十工段站受灾程度相对较轻。2) 2020年2~6月中，变温中度灾害在2月受灾程度最重，5月的受灾程度最轻。十七工段站在2~6月的受灾程度均为最重，其次为南沙村站；东沙村站、群围村站和长兴村站受灾程度相对较轻。3) 变温重度灾害在2020年2~6月中，3月的灾害发生略高于其他月份，4月和5月无重度灾害发生；在龙虎村站重度灾害发生略多于其余站点，绝大部分站点无重度灾害发生。4) 对虾基地站1 h变温轻度灾害在0~1.57次之间，最大受灾次数1.57次出现在4月08 h，且1 h变温轻度灾害在对虾基地站呈现出双峰分布的特点。1 h变温中度灾害在0~0.29次之间，且多发于午后，最大受灾次数0.29发生于3月16 h、4月17 h和7月的16 h和18 h。1 h变温重度灾害次数在0~0.14次之间，多发于7月。5) 十七工段站1 h变温轻度灾害在0~5.29次之间，最大受灾次数5.29次出现在4月07 h，且1 h变温轻度灾害在十七工段站呈现除双峰分布的特点。1 h变温中度灾害在0~1次之间，最大受灾次数1发生于11月16 h，且变温中度灾害在十七工段站也呈现除双峰分布特点。1 h变温重度灾害在0~0.41次之间，最大受灾次数发生在7月15 h，且重度灾害多发于7月。
system: [SYSTEM]

user: 1. 引言 南美白对虾是当今世界上三大养殖对虾中但产量最高的虾种 [ 1 ] [ 2 ]，南美白对虾肉质鲜美深受广大市民的喜爱 [ 3 ]，且适温范围广，是一种广温广盐虾，具有适应性强、生长速度快、抗病能力强三大特点，在我国的养殖业有重要地位 [ 4 ] [ 5 ] [ 6 ]，带动了我国水产养殖业的发展，促进沿海地区的经济发展，同时带动了冷藏加工、育苗、饵料生产等相关产业链的发展。 对于虾苗的放养，为了确保养成能获得良好的效果，除了要求虾苗的数量和质量要完全符合养成生产的要求外，由于虾苗个体小，身体脆弱，对环境条件的变化十分敏感，当外界的环境条件变化时 [ 7 ] [ 8 ]，对白对虾的生长发育，具有极大的影响，当环境变化剧烈，白对虾容易产生应激反应，使其自生的机体功能和免疫力下降，容易诱发疾病，使得白对虾产量急剧下降，造成养殖业的经济损失。变温是引起白对虾病害的重要因素，白对虾为变温动物，对环境的变化十分敏感，自身温度随环境变化而变化，当外界的气温变化过大、变化速度过快，白对虾无法快速调节，从而使白对虾产量急剧下降。所以研究变温灾害的时空分布对预防白对虾减产、减少养殖业的损失有重要意义 [ 9 ] [ 10 ] [ 11 ] [ 12 ] [ 13 ]。 萧山区是我国白对虾养殖的主要地区之一 [ 14 ] [ 15 ] [ 16 ]，拥有白对虾养殖规模7万余亩，南美白对虾的养殖推动了杭州经济的发展，白对虾主要是淡水养殖，盐度对其生长也有显著影响(黄凯等，2004)。白对虾的淡水养殖近年来表现出些许劣势，白对虾的生长受水质和气候、气象灾害影响较大，人工养殖密度也是一个影响因素。在萧山区变温灾害时有发生，大大影响了本地的白对虾养殖，造成了严重的经济损失。所以本文重点就萧山区变温灾害的时间变化规律以及空间分布特征做出分析。
system: [SYSTEM]

user: 2. 资料与方法 2.1. 研究区概况 萧山区位于浙江省北部，钱塘江南岸，属杭州市，地理坐标在29˚50'54''N~30˚23'47''N，120˚04'22'' E~120˚43'46'' E。萧山区地处浙东低山丘陵区北部、浙北平原区南部。地势南高北低，自西南向东北倾斜，中部略呈低洼。地貌分区特征较为明显：南部为低山丘陵地区，间有小块河谷平原；中部和北部为平原，中部间有丘陵，全区平原约占66%，山地占17%，水面占17%。 萧山区位于北亚热带季风气候区南缘。总的气候特征为：冬夏长、春秋短，四季分明；光照充足，雨量充沛，温暖湿润；冷空气易进难出，灾害性天气较多；光、温、水的地域差异明显。年平均气温为16.1℃。年平均地面温度为18.3℃。 2.2. 白对虾养殖概况 南美白对虾养殖一般一年两茬，分为大棚养殖和外塘养殖两种模式。大棚养殖第一茬放苗时间一般在3月中旬到4月初，5月5到20日出塘，养殖周期为45到56天；第二茬放苗时间一般在6月中旬到7月份初，9月中旬到10月初出塘，养殖周期为3个月左右。外塘养殖即南美白对虾幼苗在外塘露天养殖，第一茬放苗时间一般在3月底到4月中旬，5月20日到6月5日出塘，繁殖周期约为56天。第二次幼苗时间通常是从6月底到7月中旬。9月25日到10月中旬出塘，养殖周期为3到4个月。 2.3. 资料概况 本文观测资料来源于萧山区气象局，根据21个区域自动站2020年1-6月的逐日逐时平均气温资料以及两个典型区域自动站：对虾基地站(白对虾集中养殖区)和十七工段站(建站时间最长的站点)从建站至2020年6月的逐日逐时平均气温资料，运用基本数理方法计算出1 h变温数据。本文研究选取的21个站点为：十七工段站、外六工段站、外二十工段站、东江闸站、丰乐村站、新围大队站、临江海港站、对虾基地站、舒兰农业站、新北桥社区站、南沙村站、党山村站、龙虎村站、幸福村站、河庄村站、群围村站、长兴村站、益民村站、三联村站、南翔村站和东沙村站。图1为21个区域自动站在杭州萧山的分布。 图1. 选取21个站点分布图 2.4. 研究方法 1) 根据白对虾生长发育确定了一小时变温的灾害指标如下： 一小时变温指标 轻度灾害(1 h变温3.0℃~4.9℃) 中度灾害(1 h变温5.0℃~6.9℃) 重度灾害(1 h变温≥7.0℃) 2) 对资料进行筛选，并根据南美白对虾的生长周期截取所需日期，并根据需要对缺测数据进行补充或者删除。 3) 利用编程计算每年、每月、每天特殊天气的小时占比。 4) 根据21站点的温度变化情况，分析不同站点的1 h变温灾害发生情况。 5) 结合前人关于南美白对虾的研究，总结不同程度1 h变温灾害对南美白对虾影响的程度。
system: [SYSTEM]

user: 5. 结论 1) 2020年2~6月中，变温轻度灾害在4月受灾程度最重，6月受灾程度最轻。在2020年，十七工段站在2~6月的受灾程度均为21个区域自动站中最重；外二十工段站在21个区域自动站中的受灾程度相对较轻。 2) 2020年2~6月中，变温中度灾害在2月受灾程度最重，5月的受灾程度最轻。在2020年，十七工段站在2~6月的受灾程度均为21个区域自动站中最重，其次为南沙村站；东沙村站、群围村站和长兴村站在21个区域自动站中的受灾程度相对较轻。 3) 变温重度灾害在2020年2~6月中，3月的灾害发生略高于其他月份，4月和5月无重度灾害发生；在龙虎村站重度灾害发生略多于其余站点，绝大部分站点无重度灾害发生。 4) 对虾基地站1 h变温轻度灾害在0~1.57次之间，最大受灾次数1.57次出现在4月08 h，且1 h变温轻度灾害在对虾基地站呈现出双峰分布的特点，在上午和午后各有一个峰值区，上午峰值在08~09 h附近，下午峰值在16~17 h附近，除4、5、11月以外，2~11月中，午后峰值高于上午。1 h变温中度灾害在0~0.29次之间，且多发于午后，最大受灾次数0.29发生于3月16 h、4月17 h和7月的16和18 h。1 h变温重度灾害次数在0~0.14次之间，多发于7月。 5) 十七工段站1 h变温轻度灾害在0~5.29次之间，最大受灾次数5.29次出现在4月07 h，且1 h变温轻度灾害在十七工段站呈现除双峰分布的特点，在上午和午后各有一个峰值区，上午峰值在06~08 h附近，午后峰值在15~17 h附近，除11月午后峰值略高于上午外，2~11月中，表现为上午峰值高于午后。 1 h变温中度灾害在0~1次之间，最大受灾次数1发生于11月16 h，且变温中度灾害在十七工段站也呈现除双峰分布特点，上午和午后各有一个峰值区，上午在07 h附近，午后在17 h附近，在上半年中(2~6月)，表现为上午峰值多大于午后，在下半年(7~11月)则相反。1 h变温重度灾害在0~0.41次之间，最大受灾次数发生在7月15 h，且重度灾害多发于7月。
system: [SYSTEM]

user: 摘要 将二维EMD方法应用于鲁西地块重力数据分解，获取反映研究区不同尺度深部地幔隆陷及其与金矿化空间分布关系的IMF 4 和IMF 3 重力分量图像。IMF 4 低通滤波重力分量图像显示鲁西地块存在1处一级区域性地幔隆起，临沂幔隆(I UY )和两处一级区域性地幔凹陷，平邑幔陷(I DP )和莱芜幔陷(I DL )。几乎所有金矿床及其相关的中生代火成岩体(侵入岩和火山岩)皆沿这种区域性幔隆或幔陷的边缘分布。IMF 3 带通滤波重力分量图像，揭示在一级区域性地幔隆陷的基础上发育的二级区域性地幔隆陷特征。在整个研究区这种二级区域性幔隆和幔陷呈带状相间分布。在研究区东部，临沂–沂南–沂水幔隆西侧分布有相同展布方向和近似规模的孟良崮–棱背岭–沂山–唐吾幔陷，在其南东侧发育临沭–莒南幔陷。在研究区中西部，分布4对呈北西向分布的地幔隆陷。自南至北，依次为张庄–滕州幔陷；向城–石井–水泉–南辛幔隆；罗庄–郑城–四海山幔陷；临沂–费县–平邑–泗水幔隆；蒙山–宫里–泰安幔陷；蒙阴–新泰–徂徕–柳埠–济南幔隆；雁翎关–跺庄幔陷；辛庄–腰关–阎家峪幔隆。研究区构造线的分布与二级幔隆和幔陷的展布方向具有高度一致性。金矿床及其相关的中生代岩浆岩(侵入岩和火山岩)绝大多数分布于二级幔隆区或其边缘。业已表明，中生代太平洋板块向欧亚板块俯冲导致的郯庐断裂大规模的左行平移引起鲁西二级区域性地幔隆陷，二级地幔隆陷导致的拉张和挤压作用引起的构造岩浆活动构成中生代金–铜金成矿的有利地质背景，发育于二级地幔隆起及其边缘的构造岩浆岩带是找寻Au(Ag)、Cu等矿床的远景地段。
system: [SYSTEM]

user: 1. 引言 重磁地球物理数据具有深穿透性，且通常含有地壳深部关键地质结构和矿化信息 [ 1 ] [ 2 ]。而这些关键信息由于传统信息提取数学模型(如傅氏变换等)的局限性而没有得到有效地提取和应用，非线性技术为提取深部地质矿化信息，进一步精确构建深部地质结构和地质体模型提供了新的机遇 [ 3 ]。 地球科学的复杂性主要涉及地质过程的多期多阶段性及其结果的叠加性，记录这一复杂过程和叠加现象的数据集往往具有非线性结构和非平稳特征，使得诸如地质统计学和傅立叶变换等常用数据处理方法，严格意义上，并不适用于处理非线性和非平稳数据 [ 4 ] - [ 9 ]。为此，Huang等 [ 10 ] [ 11 ] [ 12 ] [ 13 ] 建立一种普适的信号分析技术，称为Hilbert-Huang变换(Hilbert-Huang Transform， HHT)。该方法包括两部分算法，即，经验模分解(Empirical mode decomposition， EMD)和Hilbert谱分析，二者有机结合定量刻画非线性和非平稳过程。该方法不依赖于核函数的先验选择，而是分解信号到起源于连续性极值的本征振荡模。该方法的创立与发展为基于非平稳数据研究致矿地质异常事件与过程开辟一个新领域。 经验模分解(EMD)可将复杂数据集分解为一系列有限个数的本征模态函数(intrinsic mode function， IMF)，这一系列的IMF可以通过筛分过程(sifting process)得到，通过筛分过程还能够从原始数据中得到局部的相对高通滤波振动。该方法适用于地震数据处理 [ 14 ]，噪声滤波 [ 15 ]，航空重力数据一维处理 [ 16 ]。Nunes等 [ 17 ] [ 18 ] 发展二维EMD (Bi-EMD，简称BEMD)方法，使之适用于处理二维图像；沈滨等 [ 19 ] 将BEMD方法应用于纹理分割及图像瞬时频率估计。目前，该方法已被广泛应用于提取区域和局部重磁致矿异常和地球化学元素组合致矿异常信息提取，取得显著效果 [ 3 ] [ 20 ] [ 21 ] [ 22 ]。近年来，许多学者应用多维分型模型(multi-fractal mode)提取重磁异常亦取得明显进展 [ 23 ] [ 24 ] [ 25 ] [ 26 ]。 本文尝试应用BEMD方法，对鲁西地块1:20万重力数据进行分解，结合其地质矿化特征，试图揭示鲁西地块深部地质结构及其与金矿化分布的关系，为金矿深部成矿预测与评价提供科学依据。
system: [SYSTEM]

user: 2. 二维EMD基本原理及算法实现 EMD方法 [ 10 ] 能够自适应地将一组数据分解为不同频率的组成——本征模态函数(intrinsic mode function， IMF)。通过筛分，可以得到从高通滤波到低通滤波的一系列的IMF分量，而每一个IMF分量包含以下性质：1) IMF分量有相同的过零点数和极值点个数或者差1；2) 在任一点，由极大值点构成的上包络和由极小值点够成的下包络的平均值趋于0。筛分过程要保证数据的极值点个数至少为2。EMD分解可以将数据分解为有限个的IMF分量和剩余分量，不同的IMF分量代表了不同的振动模式(oscillatory mode)，而剩余分量描述了数据的整体趋势。 2.1. 二维EMD的基本原理 在二维空间中，对一维EMD方法进行推广，实现了二维的EMD的分解(BEMD)。二维EMD方法过程与一维EMD类似，主要的区别在于二维EMD方法矩阵中极值点的确定和包络面的拟合比一维EMD方法更复杂(后面将详细说明)。令 O r i ( m ， n ) 为待分解的二维数据，通过二维筛分过程，可将二维数据分解为有限个的二维IMF分量(BIMF)，分别代表二维数据的不同频率(尺度)的结构特征，按照频率的由高到低，依次为 B 1 ( m ， n ) ， B 2 ( m ， n ) ， ⋯ ， B t ( m ， n ) ，则有： O r i ( m ， n ) = ∑ i = 1 t B i ( m ， n ) + R e s ( m ， n ) (1) 其中， B i ( m ， n ) 为第i个二维IMF分量， R e s ( m ， n ) 为剩余分量。在滤波过程中，与一维类似，可以设计不同的滤波器 S H P ( m ， n ) ， S B P ( m ， n ) ， S L P ( m ， n ) ，分别用于高通、带通和低通滤波。也可以选择性的选取某些(个)反映特定频率(尺度)结构特征的二维IMF分量作为滤波结果。 S H P ( m ， n ) = ∑ i = 1 k B i ( m ， n ) (2) S B P ( m ， n ) = ∑ i = k p B i ( m ， n ) (3) S L P ( m ， n ) = ∑ i = p t B i ( m ， n ) + R e s ( m ， n ) (4) 换句话说，我们能通过长周期分量相加获得低通滤波结果，通过短周期分量相加获得高通滤波结果，通过选择性周期分量相加获取带通滤波结果。这就是该方法能够用于信号分解和致矿异常信息提取的方法学基础。 此外需设定筛分过程的停止条件，这里主要通过限定SD的大小来得到，SD由两个相邻的筛分过程的结果得到 [ 13 ] [ 17 ] [ 18 ]，在二维中，设定： S D i j = ∑ m = 1 ∑ n = 1 | h j ( i − 1 ) ( m ， n ) − h j i ( m ， n ) | 2 h j ( i − 1 ) 2 ( m ， n ) < ε (5) SD的设定对BIMF有影响，如果SD的值较小，BIMF的个数将增加，筛选过程也更为细致，但运算时间会有所增加。但在二维中， ε 的取值大小还没有一个标准，往往通过经验来判断和设定，亦可借助局部均值矩阵 m e a n ( m ， n ) 的期望值、方差获取。 2.2. 算法实现 二维EMD的流程步骤与一维EMD类似，主要是将一维的极值求取和包络线的拟合变为二维的极值求取与二维包络面的拟合。二维EMD的算法流程(图1)如下： O r i ( m ， n ) 为待分解的二维数据， 1) 初始化：令 r 0 ( m ， n ) = O r i ( m ， n ) ，并 j = 1 (j为二维IMF的标号)； 2) 提取第j个IMF分量 ① 令 h 0 ( m ， n ) = r j − 1 ( m ， n ) ， i = 1 ， ② 查找二维数据 h i − 1 ( m ， n ) 的极大(小)值， ③ 计算上(下)包络 u p p e r i − 1 ( m ， n ) ( l o w e r i − 1 ( m ， n ) )， ④ 计算均值 m e a n i − 1 ( m ， n ) = ( u p p e r i − 1 ( m ， n ) + l o w e r i − 1 ( m ， n ) ) / 2 ， ⑤ 令 h i ( m ， n ) = h i − 1 ( m ， n ) − m e a n i − 1 ( m ， n ) ， i = i + 1 ， ⑥ 验证结束条件， ⑦ 循环执行步骤②~⑥，直到满足停止条件为止跳出，而此时的 h i ( m ， n ) 即为第j个二维IMF分量 s j ( m ， n ) ， 3) 剩余量 r j ( m ， n ) = r j − 1 ( m ， n ) − s j ( m ， n ) ， 4) 令 j = j + 1 ，循环步骤2-4，直到剩余量 r j ( m ， n ) 的极值数小于2为止。
system: [SYSTEM]

user: 3. 二维EMD方法在鲁西地块重力致矿异常信息提取中的应用 前已述及，重磁信息的一个显著特点是具有“透视性”，它不仅能够反映浅部的地质现象，而且通过对重磁场的分解，还能够获取深部地质结构信息。重磁信息的另一个特点是具有多解性，这是因为通常我们获取的重磁数据测量的是不同规模、不同深度、不同形态和不同磁性(密度)地质体组合的叠加场。这就要求我们根据地质体场的性质和特点，借助于信息处理技术实现叠加场的分解，尽可能使场与地质体一一对应，并结合地质矿化信息等约束条件，获取目标信息，最终达到解决地质与找矿疑难问题之目的 [ 27 ] [ 28 ]。 图1. 经验模分解法流程图(根据Chen等，2017) [ 3 ] 3.1. 鲁西地块地质矿化特征 鲁西地区系指沂沭断裂以西的基岩隆起区。区内众多的金矿点和一些金矿床初步划分为三种成因类型：1) 变质热液金矿床；2) 与侵入岩有关的热液金矿床；3) 火山热液金矿床(图2)。与金矿成矿有关的主要因素如下： a) 断裂：以NE向沂沭断裂及NW向断裂为主的断裂系统。在沂沭断裂带内或其附近的中生代火山岩盆地，一些金矿点受其控制呈NE向分布。在沂沭断裂带以西，NW向断裂系统控制了前寒武纪变质岩系、古生代碳酸盐岩沉积盆地、中生代火山岩盆地和中生代侵入体以及金矿床、矿点的分布。因此，NE向沂沭断裂带及NW向断裂系统构成了研究区的一级控矿地质因素 [ 28 ]。 b) 地层：本区前寒武纪结晶基底由表壳岩(5%)和深成侵入岩(95%)组成。表壳岩主要包括中太古代(>3 000 Ma)沂水岩群和新太古代(2800 Ma)泰山岩群。沂水岩群分布于沂水断裂带，其主要岩性为斜长二辉麻粒岩、二辉角闪斜长片麻岩和斜长角闪岩以及含紫苏磁铁石英岩等一套深变质岩系；其原岩为超镁铁质–镁铁质熔岩、凝灰岩及泥质砂岩夹硅铁质岩石。泰山岩群分布于沂沭断裂带以西的鲁西隆起区，其主要岩性为斜长角闪岩，透闪片岩、角闪黑云变粒岩夹铁闪磁铁石英岩；其原岩为超镁铁质–镁铁质熔岩(科马提岩)、凝灰岩、泥质粉砂岩和硅铁质岩；泰山岩群是我国保存最好、发育最完整的典型新太古代绿岩带 [ 29 ]。 c) 中生代岩体：与岩体有关的金矿床是本区金的主要工业矿床类型。与金矿有关的侵入体皆为中生代浅成斑状杂岩体，成岩时代为190~125 Ma [ 30 ] [ 31 ] [ 32 ]，发育于隐伏基底(这里指被古生代沉积碳酸盐岩覆盖的前寒武结晶基底区)区。金矿床(矿点)通常环绕侵入体分布(图2)。这类金的工业矿床的形成归结于矿源(前寒武纪结晶基底)，热源(中生代侵入岩)，有利围岩(古生代沉积碳酸盐岩)和成矿构造等成矿要素的最佳匹配 [ 28 ]。 图2. 鲁西地块地质矿产图(据Chen et al， 2001修编) [ 28 ] 3.2. 鲁西地块重力致矿异常信息提取 论文研究所用重力数据来源于山东地质矿产局第物化探大队1980年代实施的1:20万地面综合物探调查中的重力测量结果。其测量网距为2 km × 2 km，工作总精度为±2.32 g.μ，控制面积约48，000 km 2 。鲁西地块出露地层密度参数由高至低依次为：泰山岩群(2.70~2.90 g/cm 3 ，平均值2.81 g/cm 3 )→寒武–奥陶系(2.64~2.73 g/cm 3 ，平均值2.70 g/cm 3 )→石炭–二叠系(2.58~2.67 g/cm 3 ，平均值2.61 g/cm 3 )→侏罗–白垩系(2.39~2.54 g/cm 3 ，平均值2.51 g/cm 3 )。侵入岩密度参数由高至低依次为：前寒武纪花岗质侵入岩(TTG岩系：2.59~2.84 g/cm 3 ，平均值2.75 g/cm 3 )→中生代花岗岩类(2.59~2.71 g/cm 3 ，平均值2.65 g/cm 3 ) [ 20 ]。 原始重力数据图像(图3)是不同深度、不同密度地质单元及其地质结构的综合反映。为揭示鲁西地块深部地质结构及其对金矿化系列形成与分布的控制，对其1：20万重力数据(图3)进行二维EMD分解，获得4个二维IMF分量，以及剩余量 R e s ( m ， n ) ， O r i ( m ， n ) = ∑ i = 1 4 B I M F i ( m ， n ) + R e s ( m ， n ) 。各个二维IMF分量代表二维数据的不同频率的结构特征。BIMF 1 、BIMF 2 、BIMF 3 和BIMF 4 按照频率从高到低产生的，对于相同的局部区域BIMF 1 带通滤波频率通常高于BIMF 2 的频率。在鲁西1:20万重力数据的处理过程中，设定SD = 0.02。 图3. 鲁西地块原始重力数据图像(数据来源山东地矿局物化探大队，1991)，地质图例同图3，AB和CD线段是对分解的不同的IMF重力分量进行正交性检验的一维重力剖面位置 3.3. IMF分量正交性检验 一维IMF满足正交性，虽然还不能在理论上进行证明，但实验结果证明一维IMF分量是正交的 [ 10 ]。假设各个一维IMF分量是正交的，则： X 2 ( t ) = ∑ i = 1 n + 1 I M F i ( t ) + 2 ∑ j = 1 n + 1 ∑ k = 1 n + 1 I M F j ( t ) I M F k ( t ) ，( j ≠ k ) (6) 中等式右边的第二项等于0 (或趋近于0)，其中 X ( t ) = ∑ i = 1 n + 1 C i ( t ) ， C i = I M F i ， i = 1 ， 2 ， ⋯ ， n ， C n + 1 = r e s (剩余分量)。在实际计算中，用IO来表示正交性指标， I O = a b s ( ∑ t = 0 T ( ∑ j = 1 n + 1 ∑ k = 1 n = 1 C j ( t ) C k ( t ) / X 2 ( t ) ) ) (7) 或者可以计算两个不同IMF分量之间的 I O f g 值： I O f g = a b s ( ∑ t C f C g C f 2 + C g 2 ) (8) 图4. 鲁西地块北东向(AB)一维重力数据分解的各个IMF分量 图5. 鲁西地块北西向(CD)一维重力数据分解的各个IMF分量 这里，T是总的时间区间，为获得IMF的最优分解，IO应该非常小，典型地，其值介于0.01和0.001之间是可以接受的。如果信号长度很短，IO的值会增大至0.05左右。线性分离系统需满足正交性，对于一些特殊的信号，Huang [ 10 ] 认为虽然分量之间不再满足正交性，但由EMD分解到的IMF分量依然是有意义的。分别详细研究图4中北东向(AB)和北西向(CD) 2个方向的二维IMF分量的一维数据结果(图4，图5)。 并对这两个方向的IMF分量进行正交性检验(表1)，结果表明两个方向的IMF分量近似满足正交性。 Table 1 Orthogonality assessment of IMF AB线段 CD线段 IO 12 0.0718 0.0228 IO 13 0.0089 0.0308 IO 14 0.0088 0.0571 IO 23 0.0163 0.0368 IO 24 0.0341 0.0379 IO 34 0.0283 0.0278 表1. IMF重力数据分量的正交性检验
system: [SYSTEM]

user: 4. 结果与讨论 在上述分解的4个分量中，其中IMF 4 (图6)和IMF 3 (图7)具有较明确地质意义。 图6. 鲁西地块IMF 4 重力分量图像，地质矿产图例同图3 图7. 鲁西地块ＢIMF 3 重力分量图像，地质矿产图例见图3 4.1. IMF4和IMF3重力分量图像的地质意义 图6是低通滤波重力分量(IMF 4 )图像。该图表明，该重力分量值的高低分布与地表岩石的密度没有对应关系，反映了鲁西地块深部区域性地幔隆起和凹陷。在该图像上，沿沂沭断裂呈北东向分布的高重力分量异常带可能反映沿临沂—沂水分布一地幔隆(或称基底隆起)起，简称临沂幔隆(I UY )，重力IMF 4 分量值分布−40~80 μm。沂沭断裂带基本上沿该幔隆分布，其中火山岩盆地和幔隆两侧的金矿床亦呈北东向分布(图6)。 分布于沂沭断裂以西的IMF 4 低重力分量异常区反映了地幔凹陷特征。该区存在两处低重力分量异常区，分别对应两处地幔凹陷。其一是分布在平邑以南呈北西向分布的凹陷，简称平邑幔陷(I DP )，重力IMF 4 分量值分布−310~−190 μm，对应于该凹陷地表出露的地质单元是由太古代TTG和少量变质岩系构成的泰山岩群。在平邑凹陷北西侧分布有呈北西向展布的平邑–费县中生代火山岩盆地、燕山期鲁西杂岩体以及相关的大型归来庄金矿床等 [ 33 ] [ 34 ]，在其南东侧分布有苍山–龙宝山花岗岩带以及相关的金矿床。其二是分布于蒙阴–新泰断裂以北以莱芜为中心呈弧(扇)形分布的地幔凹陷，简称莱芜幔陷(I DL )；其重力IMF 4 分量值分布−310~−190 μm。对应地表出露的地质单元内弧为太古代TTG岩系 [ 35 ]，外弧为古生代碳酸盐岩。幔陷内自北向南发育三条呈弧形分布的中生代岩浆岩(侵入岩和火山岩)带，且不同程度的发育金矿化(图7)。此外，垂直于弧形幔陷的轴线发育放射性断裂。 总之，幔隆和幔陷从整体上控制研究区基本地质格局以及中生代与火成岩有关的金矿床和相关火成岩(侵入岩和火山岩)体的区域分布。 图7是带通滤波重力分量(IMF 3 )图像。该图表明，该重力分量值的高低分布与地表岩石的密度没有对应关系，可能反映在一级区域性地幔隆陷基础上发育二级区域性地幔隆陷 [ 3 ] [ 36 ]。在整个研究区这种二级区域性幔隆和幔陷呈带状相间分布。在研究区东部，临沂–沂南–沂水幔隆西侧分布有相同展布方向和近似规模的孟良崮–棱背岭–沂山–唐吾幔陷，在其南东侧发育临沭–莒南幔陷。在研究区中西部，存在4对呈北西向分布的地幔隆陷。自南至北依次为张庄–滕州幔陷；向城–石井–水泉–南辛幔隆；罗庄–郑城–四海山幔陷；临沂–费县–平邑–泗水幔隆；蒙山–宫里–泰安幔陷；蒙阴–新泰–徂徕–柳埠–济南幔隆；雁翎关–跺庄幔陷；辛庄–腰关–阎家峪幔隆。构造线的分布与幔隆和幔陷的展布方向具有高度一致性。绝大多数中生代火山岩及其相关的金矿床和铜金矿床分布于幔隆区，而绝大多数中生代侵入岩及其相关的金矿床和铜金矿床则分布于幔隆的边缘。 4.2. 金矿化及其相关火成岩形成的构造环境 一些学者 [ 37 ] [ 38 ] [ 39 ] [ 40 ] 的研究资料表明：中生代太平洋板块向欧亚板块俯冲，导致鲁西地区的整体隆起；伴随郯庐断裂大规模的左行平移，在其西侧形成了一系列NW-NNW向断裂和陆内断陷盆地，在以张应力为主的构造环境下，出现了大量火山喷发和浅成侵入活动，从而形成了一系列火山岩盆地，和浅成次火山斑状杂岩体。这种发育中生代火成岩的以张应力为主的构造环境恰是一级幔陷中的局部隆起所致。王先美等 [ 41 ] [ 42 ] 由出露于沂沭断裂带、鲁西地体、鲁东地体的中生代地层、岩浆岩，结合断裂活动年代学、区域地质等资料分析，将沂沭断裂带晚中生代构造演化划分出距今约160 Ma、130~110 Ma、90~80 Ma等3个关键时期，并分别与左行压剪、左行张剪、右行压剪构造活动相对应。晚中生代北北东走向的沂沭断裂带与北西向断裂系的共轭属性及其构造活动方式，是在华北地块周边不同板(地)块、不同时期的汇聚与离散等相对运动，以及地球内部不同圈层的物质与能量交换的背景下发生的。沂沭断裂带早期左行压剪与北西向断裂系的右行压剪构造活动的动力来自华北、华南地块的陆–陆碰撞；沂沭断裂带早白垩世左行张剪运动与北西向断裂系同期右行张剪活动是Izanagi板块俯冲与华北岩石圈拆沉引起上地壳快速伸展的联合效应；Kula板块晚白垩世向西的快速俯冲可能构成华北地块挤压隆升的主要动力，沂沭断裂带与北西向断裂系则相应表现出右行压剪与左行压剪构造活动。 某些学者 [ 43 ] [ 44 ] [ 45 ] [ 46 ] 根据同位素定年，微量元素，稳定同位素和稀土元素等分别对鲁西铜石、铜井、莱芜、龙宝山、莲子旺、粟园等杂岩体进行了较系统研究。其结果表明，其成岩时代分布112.5~189.8 Ma，众数值是110~126 Ma，这一时期与前述的沂沭断裂左行张剪相对应，亦与早白垩世中国东部岩浆、成矿等作用最为强烈的时期(120~130 Ma)相吻合 [ 47 ] [ 48 ] [ 49 ] [ 50 ] [ 51 ]，亦很可能对应岩石圈减薄的高峰时期 [ 52 ]。 林景仟等 [ 44 ] 通过对鲁西地区中生代火成活动的研究认为：鲁西地区的火成活动远与太平洋板块对欧亚板块的俯冲导致沂沭断裂的活动有密切联系。此时，在鲁西地块上北西向的断裂频频活动，该组断裂的强烈拉张最早发生在早侏罗世，发育于鲁西的南部，生成了高钾钙碱性–高钾碱性的岩石组合(如铜石杂岩体)，至早白垩世地壳拉张活动加剧，南部区域的北西向断裂切割深度大，从上地幔的深层生成了碱性岩浆，北部区域北西向断裂亦切割到上地幔，生成了钙碱性的辉长质岩浆，它与上部壳层产生的熔体发生了混染，形成了二辉闪长质岩浆(较高水条件下生成了角闪闪长质岩浆)，二长质岩浆是在混入了更多的酸性熔融体的情况下产生的。在此时期，近沂沭断裂带的地段，早白垩世较早阶段生成的辉长质岩浆发生混染生成闪长质岩浆，上部的壳层则产生了花岗质岩浆，后一组岩石多在沂沭断裂带近旁分布，它们亦沿北西向断裂活动，分布于较远离沂沭断裂带的部位，如滕州桑村岩体、蒙阴的虎头崖岩体等。 形成杂岩体的岩浆最初可能起源于上地幔，在其上升过程中不同程度地同化了太古代基底岩层，使部分Au、Cu、F、B等物质进入岩浆。随着岩浆的进一步侵位，冷凝结晶(局部伴有交代作用)，成矿组分(Au、Cu等)和矿化剂(F、B等)在富碱质的热液中进一步富集形成液态矿源，并在岩体内外有利的构造部位沉淀成矿，最终形成火成岩体(侵入岩和火山岩)和金与金–铜矿化 [ 53 ] [ 54 ] [ 55 ]。 鲁西在中生代除发育花岗质岩石及其相关的金(铜)矿床外，还发育基性火成岩 [ 56 ]。典型玄武岩发现于费县方城一带，K-Ar年龄为119~125 Ma [ 56 ] [ 57 ]。鲁西中生代火山岩的主体为钾玄质火山岩，蒙阴盆地钾玄质火山岩的 40 Ar- 39 Ar年龄变化于114.8~124.3 Ma之间 [ 58 ]。谭东娟和林景仟 [ 59 ] 测得济南辉长质杂岩体的 40 Ar- 39 Ar年龄为115 Ma，山东省地质矿产局(1991)测得邹平盆地侵入青山组火山岩的辉长质岩体 40 Ar- 39 Ar年龄为128~130 Ma [ 60 ]。它们皆属于燕山晚期岩浆活动的产物。 上述研究表明：中生代太平洋板块向欧亚板块俯冲导致的郯庐断裂大规模的左行平移引起鲁西地幔隆陷，地幔隆陷导致的拉张和挤压作用引起的构造岩浆活动是中生代金、铜成矿的有利背景，亦是该类矿床找矿的远景地段。
system: [SYSTEM]

user: 5. 结论 二维经验模分解(BEMD)方法被成功地应用于鲁西地块重力数据分解，获取反映鲁西地块不同尺度深部地幔隆陷及其与金矿化空间分布关系的IMF 4 和IMF 3 重力分量图像。 1) IMF 4 低通滤波重力分量图像显示鲁西地块存在1处一级区域性地幔隆起，临沂幔隆(I)和两处一级区域性地幔凹陷，平邑幔陷(I DP )和莱芜幔陷(I DL )。几乎所有金矿床及其相关的中生代火成岩(侵入岩和火山岩)皆沿幔隆或幔陷的边缘分布。此外，研究区北部近于平行分布的三个弧形花岗岩带的形成与分布受莱芜弧形幔陷的控制，且不同程度的发育金矿化。 2) IMF 3 带通滤波重力分量图像，揭示在一级区域性地幔隆陷的基础上发育的二级区域性地幔隆陷特征。在整个研究区这种二级区域性幔隆和幔陷呈带状相间分布。在研究区东部，临沂–沂南–沂水幔隆西侧分布有相同展布方向和近似规模的孟良崮–棱背岭–沂山–唐吾幔陷，在其南东侧发育临沭–莒南幔陷。在研究区中西部，分布4对呈北西向分布的地幔隆陷。自南至北，依次为张庄–滕州幔陷；向城–石井–水泉–南辛幔隆；罗庄–郑城–四海山幔陷；临沂–费县–平邑–泗水幔隆；蒙山–宫里–泰安幔陷；蒙阴–新泰–徂徕–柳埠–济南幔隆；雁翎关–跺庄幔陷；辛庄–腰关–阎家峪幔隆。研究区构造线的分布与二级幔隆和幔陷的展布方向具有高度一致性。金矿床及其相关的中生代岩浆岩(侵入岩和火山岩)绝大多数分布于二级幔隆区或其边缘。 3) 中生代太平洋板块向欧亚板块俯冲导致的沂沭断裂大规模的左行平移导致鲁西地块在一级地幔凹陷的基础上发育二级区域性地幔隆陷。张应力为主的构造环境中形成的中生代火成岩恰是二级幔隆所致；深部地质结构对金矿化的形成与分布的控制是通过二级(局部)地幔隆起对与金–铜矿化有关的中生代火成岩的控制而实现的。因此，发育于二级地幔隆起及其边缘的中酸性岩浆岩区是找寻金–铜矿床的远景地段。
system: [SYSTEM]

user: 摘要 阿尔金山首次发现恐龙足印化石，足印发育在侏罗系的粉砂岩岩壁上，为两类恐龙所遗留，一类为带足指的兽脚类恐龙，体长约在6 m左右，行动较快，属肉食类龙，另一类为不带足指的蜥脚类恐龙，长度在16.5 m以上，行动迟缓，具有草食性龙的特点。恐龙生活的环境为三角洲和湖泊的交汇带，这里动植物繁盛，水草肥美，侏罗纪潮湿的气候环境为恐龙提供了良好的生存和生活环境。恐龙足印的发现对于研究该地区侏罗系古环境、古地理具有重要意义。
system: [SYSTEM]

user: 1. 引言 恐龙活动痕迹主要表现为三种方式：一为遗迹(足迹化石)，二为遗物化石(恐龙蛋化石和粪便化石)，三为恐龙遗体的骨骼化石。这三种痕迹化石都直接反映了当时恐龙活动和生活的证据，也是认识和研究恐龙的主要实物资料。 中国恐龙化石研究记录可以追溯到解放前，国外科学家开展了部分研究工作，1929年Chardin等(Chardin et al.， 1929)对陕西神木发现的恐龙足迹化石进行了首次的报道。解放后国内专家对恐龙开展了大量的研究。杨钟健(Yang Zhongjian， 1960， 1966， 1979)等对陕西铜川、河北滦县、云南西双版纳傣族自治州等地的恐龙足印进行了研究和初步的总结 [ 1 ] [ 2 ] [ 3 ]。赵资奎对河北内乡的恐龙蛋化石类型和恐龙足印化石类型进行了研究 [ 4 ]。杨兴隆等对四川的恐龙化石足印进行了研究和总结 [ 5 ]。余心起等研究了安徽黄山的恐龙足印化石 [ 6 ]。甄朔南等对中国的恐龙足印化石进行了总结研究 [ 7 ]。蔡雄飞等研究了甘肃兰州民和盆地盐锅峡一带恐龙足迹 [ 8 ]。李日辉等研究了山东莱阳盆地莱阳群恐龙足迹化石 [ 9 ]。李大庆等(Li Daqin et al.， 2000)和宋杰己(Song Jieji， 1992)对民和盆地恐龙足迹进行了研究 [ 10 ]。王潞琳等(Wang Luling et al.， 2009)对内蒙巴彦卓儿地质公园内的恐龙化石进行了研究。胡松梅等(Hu Songmei et al.， 2011)对陕西商洛的足迹脚印进行了研究。2009年甘肃日报报道了甘肃省地矿局第三地质矿产勘察院李大庆等在甘肃酒泉地区发现恐龙化石，并进行了研究 [ 11 ]。另外在新疆鄯善地区发现大量恐龙足迹和世界上最大恐龙化石。这些研究表明，恐龙足迹和化石在我国广泛分布，并且研究越来越受到重视。 本次发现的恐龙足迹位于柴达木盆地阿尔金山西段吐拉盆地的格斯沟阳光煤矿区(图1)，坐标：东经88˚48'43''，北纬38˚10'10.2''，海拔4770米，该区山顶终年积雪，周边有小型冰川发育，恐龙足迹化石有三组二种类型，调研表明，这是在阿尔金山地区首次发现恐龙足迹化石。