河北省崇礼幅1∶50 000矿产地质图数据库

根据《固体矿产地质调查技术要求(1∶50 000)》(DD 2019-02)和行业其他规范要求,完成河北省崇礼幅1∶50 000专项地质填图379.93 km²、矿点检查11处,采集水系沉积物样品1622件、重力测量物理点1623个,测试锆石U-Pb年龄样3组、流体包裹体样80余件,系统编制了河北省崇礼幅1∶50 000矿产地质图,并建立了数据库。突出了对燕山期岩浆岩和相关脉岩、碱性杂岩体和谷嘴子岩组、矿化蚀变等相关成矿信息的表达,厘定了崇礼幅中太古代、新太古代的变质深成侵入体,划分了海西期、印支期、燕山期的岩浆岩建造类型,划分了张家口组火山岩的7个建造类型、谷嘴子岩组变质岩的3个建造类型,系统总结了该区金矿的成矿地质体、成矿构造及成矿结构面、成矿作用特征标志。该数据库为该区金矿的找矿工作提供了基础数据支撑,对金矿目标层地质信息和构造建造信息的提取提供了便利条件,对今后矿产勘查部署具有一定的参考价值。     

基于不等式约束的三维电阻率探测混合反演方法

三维电阻率探测的线性反演和非线性反演中均存在着多解性的固有难题.电阻率线性反演方法的效率较高,但反演结果对初始模型的依赖性较强,易陷入局部极小;而非线性反演方法不依赖初始模型,但搜索效率极低,尚未见到关于三维电阻率非线性反演的文献.针对上述问题,融合线性与非线性反演方法的互补优势,提出了最小二乘法（线性方法）与改进遗传算法（非线性方法）相结合的混合反演方法的概念和思想.首先,提出了将介质电阻率变化范围作为不等式约束引入反演方程的思路,以实现压制多解性、提高可靠性的目标.提出了宽松不等式约束和基于钻孔推断的局部严格不等式约束的获取及定义方法.在此基础上,分别提出了基于不等式约束的最小二乘线性反演方法和遗传算法非线性反演方法.其次,对于遗传算法在变异搜索方向控制、初始群体产生等方面进行了改进,优化了其搜索方向和初始群体多样性.然后,提出了混合反演方法及其实现方案,利用改进遗传算法进行第一阶段反演,发挥其对初始模型的依赖程度低的优势,搜索到最优解附近的空间,输出当前最优个体;利用最小二乘法进行第二阶段反演,将遗传算法得到的当前最优个体作为初始模型,在最优解附近空间执行高效率的局部线性搜索,最终实现地电结构的三维成像.最后,开展了合成数据与实际工程算例验证,与传统最小二乘方法进行了对比,发现混合反演方法在压制多解性、摆脱初始模型依赖和提高反演效果方面有较好效果.

     

用天然地震探测青藏高原中部地壳、上地幔结构

从西藏南部的定日、嘎拉至青海铜铁山的天然地震探测剖面,实际路线长约２０００ｋｍ,布设了约１１０台便携地震仪,记录了数百次远震和近震事件,采用多种方法进行了资料处理与解释．依据ＳＫＳ,ＰＫＳ,ｐｓ等横波分裂特征计算的青藏高原中部上地幔的地震各向异性表明：研究区各构造单元内的地震各向异性有明显变化,发现上地幔各向异性快速波的偏振方向与造山带的走向不完全一致．在雅江缝合线、崩错－嘉黎、唐古拉山口、昆仑山口几条断裂带处南、北各向异性出现显著的差异,而金沙江缝合线和班公－怒江缝合线的南、北则没有明显的各向异性变化．由Ｐ波走时残差,利用层析技术反演了４００ｋｍ深度内的速度图像,可以看出近地表１００ｋｍ范围内速度的不均匀变化与地表划分的构造单元很吻合,进一步佐证了青藏高原是由不同时期的微板块拼合而成的认识．在青藏高原中部１５０ｋｍ深度以下发现了多处低速区．在金沙江缝合带下方约２００ｋｍ深度处有一长２５０ｋｍ以上、延伸１５０ｋｍ的低速体,推测可能是一地幔柱．利用ＰＳ转换波划分的界面,显示出青藏高原北部具有低速层和高速层交替出现的地壳结构．     

初夏西北太平洋副高东西变动对中国南部降水东西差异的影响

西北太平洋副热带高压（以下简称副高）是影响中国气候的大尺度环流系统,为了进一步了解副高对中国气候的影响,本文利用站点观测资料和大气环流再分析资料,通过资料诊断分析和数值模拟方法,探讨了6月副高东西变动对中国南部降水的影响,以及影响副高东西变动的前期海洋因子.结果表明副高东西变动对中国西南和华南地区降水的影响明显不同：副高偏东有利于降水西南偏多而华南偏少,偏西则降水变化刚好相反.其原因与副高东西变化引起的环流差异有关,华南降水与副高东（西）变动时西太平洋地区副高西北侧的东北（西南）风异常以及东亚中低纬度地区异常经向波列的变化直接有关,而西南降水异常不仅与副高东西变动在东南亚地区引起的纬向风异常有关,与青藏高原大地形动力作用对副高北侧异常纬向风的变化也有十分密切的联系.此外,副高东西变动时影响西南和华南地区的水汽来源不同,影响西南的水汽主要来源于赤道印度洋80°E附近越赤道气流,而影响华南的水汽主要来源于副高南侧偏东气流从西北太平洋地区输送的水汽.进一步分析发现前期冬春季热带西北太平洋和赤道西太平洋海温变化的偶极差异与后期初夏副高东西变动有密切联系,冬春季西北太平洋暖海温和赤道西太平洋冷海温变化有利于后期初夏副高偏西,相反则有利于副高偏东,数值模拟结果在一定程度证实了资料诊断分析结果.

     

养马岛附近海域藻华期间有色溶解有机质的生物可利用性研究

本研究利用吸收光谱和荧光激发-发射矩阵光谱-平行因子分析(EEMs-PARAFAC),研究了养马岛附近海域海水中有色溶解有机质(CDOM)的浓度、组成、来源和生物可利用性,并估算了浮游植物生长繁殖对CDOM及具有生物可利用性CDOM的贡献。结果表明,表、底层海水中CDOM浓度(以吸收系数a₃₅₀计)平均值分别为1.62±0.42 m^-1和1.30±0.47 m^-1,光谱斜率(S_275-295)平均值分别为0.022±0.003 nm^-1和0.023±0.003 nm^-1。利用PARAFAC模型识别出4种荧光组分,分别为陆源类腐殖酸C1、类色氨酸C2、类酪氨酸C3和微生物源类腐殖酸C4。荧光指数(FIX)、腐殖化指数(HIX)和生物指数(BIX)显示,CDOM受陆源输入和海洋自生源的综合影响。降解实验结果显示,表、底层海水中生物可利用性CDOM百分比(%△a₃₅₀)平均值分别为(23.36%±17.94%)和(8.93%±20.30%)。C1、C2和C4组分的荧光强度在培养之后降低,而C3组分的荧光强度上升。各荧光组分生物可利用性依次递减的顺序为:%△C1(23.75%±8.96%)＞%△C4(20.83%±11.71%)＞%△C2(11.67%±38.87%)＞%△C3(-29.61%±39.90%),显示培养之后CDOM的平均分子量和腐殖化程度降低。表层海水中a₃₅₀、%△a350与Chl a之间存在显著线性相关关系,据此可以估算出浮游植物生长繁殖对CDOM的贡献为36.9%,对具有生物可利用性CDOM的贡献为85.0%。     

利用磷灰石裂变径迹资料反演合肥盆地古地温和估计沉降率与剥蚀率

利用合肥盆地内６个磷灰石裂变径迹样品资料,反演模拟了该盆地自侏罗纪晚期以来各时代地层古地温变化,估算了沉降率与剥蚀率．模拟结果与其他地质资料推论一致,它揭示出该盆地南北两地存在不同的构造变化和受热史,反映了大别山构造发展对盆地南北两地区影响的差异．盆地南部靠近大别山地区的晚侏罗世地层在白垩纪早期埋藏温度曾大于１２０℃;早白垩世后期的构造抬升（剥蚀率约１３０ｍ／Ｍａ）使温度降至３０—４０℃;自白垩纪后期始,该地区一直处于动荡的但总体为持续抬升的构造环境中．盆地北部地区晚侏罗世与早白垩世早期地层在白垩纪期间埋藏温度曾达到和超过１００℃,但随后的大幅度构造抬升（剥蚀率约１９０ｍ／Ｍａ）使其温度降至３０—６０℃;第三纪早期,局部区域的裂陷（沉降率约６０ｍ／Ｍａ）使温度又升至８０℃左右．指出合肥盆地构造演化大体可分形成、隆升、局部裂陷和再隆升４个阶段．     

Detection of nighttime sea fog/stratus over the Huanghai Sea using MTSAT-1R IR data

A dual channel difference (DCD) method is applied to detect nighttime sea fog/stratus over the Huanghai Sea using the infrared (IR) data of shortwave (3.5–4.0 μm) and longwave (10.3–11.3 μm) channels from the Multi-functional Transport Satellite (MTSAT)-1R, i.e., shortwave minus longwave brightness temperature difference (SLTD). Twenty-four sea fog events over the Huanghai Sea during March to July of 2006 and 2007 are chosen to determine a suitable value of SLTD for nighttime sea fog/stratus detection, and ...     

GRAPES区域扰动预报模式动力框架设计及检验

设计了适用于四维变分同化系统的扰动预报模式GRAPES_PF。根据GRAPES的地形追随坐标非静力原始方程组,采用小扰动分离方法推导微分形式的线性扰动预报方程组,并利用与GRAPES非线性模式相似的数值求解方案求解线性扰动微分方程组。在设计扰动预报模式时采用了两个时间层半隐式半拉格朗日方案对动量方程、热力学方程、水汽方程和连续方程进行时间差分,空间差分方案的变量分布水平方向采用Arakawa C跳点网格,垂直方向采用Charney/Phillips跳层。利用代数消元法进一步推导得到只包含未来时刻扰动Exner气压的亥姆霍兹方程,进而通过广义共轭余差法(GCR)求解,在此基础上得到未来时刻扰动量的预报值。基于所开发的扰动模式开展了数值试验。首先在非线性模式中施加一个中尺度初始扰动高压,得到初始扰动在非线性模式中的后续演变,然后将相同的初始扰动作为扰动模式的初值进行时间积分,将扰动模式预报的结果与非线性模式的结果做了对比。结果表明,所开发的扰动模式GRAPES_PF较好地模拟了惯性重力内波的传播过程:初始高压扰动激发了一个迅速向外传播的惯性重力内波,...     

完全非线性波浪与二维固定结构物作用的IGN-BEM耦合分析模型

为高效准确地对完全非线性波浪与二维固定结构物的相互作用进行模拟分析,建立了二维完全非线性时域耦合模型。耦合模型将计算域划分为靠近结构物的内域和远离结构物的外域,每个区域均采用满足完全非线性自由水面边界条件的波浪模型进行求解。在内域使用Laplace方程描述流体运动并采用高阶边界元法（BEM）对其进行求解;而在没有结构物的外域,波浪运动的控制方程为Irrotational Green-Naghdi（IGN）方程并采用有限元法（FEM）对其进行求解。内域和外域通过一段重叠区域进行耦合,从而实现模型间变量的传递。首先利用耦合模型分别对规则波的传播、直墙前立波的生成以及相关物理模型试验进行模拟,数值结果与精确解和试验结果的良好吻合验证了耦合模型耦合方式的合理性以及处理非线性问题的准确性;然后使用耦合模型模拟分析了波浪与固定结构物间的相互作用,并将结果与线性解析解以及完全非线性BEM模型的结果进行了对比分析,进一步证明了耦合模型的正确性与高效性。     

帕米尔东北缘地壳结构的P波接收函数研究

利用位于新疆帕米尔东北缘地带12个固定数字地震台和天山动力学Ⅱ期10个流动宽频带数字地震台记录的高质量远震波形数据,应用接收函数H-κ叠加方法研究了帕米尔东北缘的地壳厚度-泊松比特征和部分台站下方的壳内界面深度.研究发现:(1) 帕米尔东北缘的Moho面起伏变化剧烈,其总体分布呈现东薄西厚、南厚北薄的特征,由塔里木盆地向天山延伸,地壳厚度约从45 km加深到55 km,从塔里木盆地向西昆仑山延伸,地壳厚度约从45 km加深到69 km;(2)沿着天山动力学Ⅱ期剖面,位于塔里木盆地北缘台站的壳内间断面的深度约为13~16 km,向北进入天山南麓加深到20 km左右,继续向北进入南天山山区壳内间断面不明显,可能暗示塔里木盆地基底向北俯冲,俯冲距离可能到达南天山的山前;(3)研究区地壳泊松比变化复杂(约从0.20到0.31),显示地壳物质组成的复杂性和显著的不均匀构造;(4)整个研究区的地壳厚度和泊松比之间没有明显的相关性,但天山动力学Ⅱ期剖面的结果表明,从塔里木盆地北缘到西南天山,地壳厚度和泊松比之间存在反相关关系,意味着天山地壳的增厚可能主要是通过以长英质岩石为主要组成成分的上地壳叠置而成;(5) 研究区全部地震台地壳厚度与海拔高程的线性回归方程表明地壳厚度与海拔的相关性相对较弱(相关系数为0.66),天山动力学Ⅱ期10个台站的地壳厚度与海拔具有很好的相关性(相关系数为0.85),可能表明沿该剖面地壳整体上处于相对均衡的状态.