非等间距GM(1,1)模型的总体最小二乘算法及其病态问题

在非等间距GM(1,1)模型中，系数矩阵中有无误差的常数项和有误差的随机项，并且系数矩阵与观测向量误差同源，即系数矩阵与观测向量中有相同的元素存在，这些相同元素应该有相同的改正数，为此本文推导了一种适合非等间距GM(1,1)模型求解的总体最小二乘算法。同时，考虑到非等间距GM(1,1)模型中存在病态问题时影响总体最小二乘计算结果的稳定性，提出对系数矩阵常数列乘以某一常数的方法，以改善病态问题。     

鄂尔多斯地块周缘垂直形变场负位错反演与应变积累研究

利用鄂尔多斯地块及其周缘1970～2014年的垂直形变速率场资料，借助负位错反演研究该区域长期应变积累。结果表明，地块东北缘山西断陷带中北段年均能量积累增量、剪应力强度都较高，西南缘六盘山断裂与渭河断裂西段次之；山西断陷带中南段至晋陕交界处年均剪应力强度较高且显示一定程度的能量积累；西秦岭构造区尤其西秦岭北缘断裂西段、晋冀蒙交界区也反映一定程度的能量积累特性。     

阿尔金断裂带中段现今构造形变模式的InSAR研究

利用2007～2010年间14景ALOS PALSAR数据及SBAS InSAR技术，获取阿尔金断裂带中段91°E附近现今地壳形变速率场，并反演该地区断层的滑动速率和闭锁深度。结果表明，阿尔金断裂中段地区的形变速率自北向南呈3个线性梯度变化区，分别为阿尔金山东段8～12 mm/a、索尔库里盆地6～7 mm/a、阿尔金断裂带以南约0 mm/a。3个速率梯度变化区主要集中在喀腊达坂断裂和阿尔金主断裂上；拟合的断层就位于金雁山南缘、喀腊达坂断裂南邻，走滑速率从西（7.1 mm/a）向东（14.0 mm/a）逐渐增大，闭锁深度自西(4.5 km)向东(10.6 km)逐渐趋深。结合前人研究推测，金雁山（阿尔金山链东部）与索尔库里拉分盆地组成的复合破裂构造模式，是转换断层运动时应力和应变调整的主要驱动机制。     

岛屿海岸线自动化提取及其时空变迁研究——以大连市獐子岛为例

海岸带是受人类活动和全球海平面上升影响的敏感地带,海岸线的提取和监测是海岸带生态系统研究和社会管理的重要内容。本文在遥感和地理信息系统的支持下,以修正的归一化水体指数(Modified Normalized Difference Water Index,MNDWI)为基础,结合遥感影像处理和直方图均衡化等技术,实现了大连市獐子岛1985—2016年海岸线的自动化提取。结果表明:(1)通过与三位专家目视解译的成果比对,本文提取海岸线的精度能满足后续研究的要求(相对误差分别为0.045%,0.032%和0.023%);(2)近30年来,獐子岛海岸线总体呈现蚀退趋势,岸线长度与岛屿面积分别呈现变短和变小的趋势,獐子岛(主岛)和大耗岛的岸线蚀退速率最大,褡裢岛次之,小耗岛最小;在人类活动较为密集的区域,海岸线呈现出较为强烈的增长趋势,海水养殖和圈海建坝是岸线增长的主要驱动力;(3)獐子岛海岸线具有显著的分形性质,分形维数随时间呈现增大的趋势,獐子岛(主岛)的分形维数最大,褡裢岛的分形维数最小。     

杭州湾中部实测波浪特性分析

利用杭州湾中部一年实测波浪资料,分析该地区的波参数统计变化特性,采用最小二乘法拟合得出相关参数之间的关系式,并统计分析了频谱特性。结果表明:杭州湾中部以小浪和轻浪为主,在各月分布较为均匀,年平均周期为2.97 s;最大一次波浪过程为冷空气影响所致,影响强度大于台风莫拉克;常浪向分布在东北至东南向,其中又以小浪出现较多,对应的波周期以2~4 s为主;强浪向主要分布在西北和东北方向,对应的波周期主要分布在3~5 s;多数特征波参数之间相关关系较好;由谱分析得知风浪占多数,且以单峰为主。研究结果可为相关工程和理论研究提供基础资料。     

利用短期静态PPP结果反演海潮负荷位移

采用香港11个GPS测站的观测资料进行1 h、2 h、3 h和4h静态PPP解算，获得4组PPP坐标序列，利用调和分析求取11个测站处8个主要分潮的负荷位移参数（振幅和相位），将其与海潮模型计算的负荷位移参数进行对比，并比较分析PPP反演值与海潮模型值改正海潮负荷信号的效果。结果表明，垂直和水平方向上，不同PPP结果反演8个分潮的负荷位移分别具有约5 mm和7 mm的差异；PPP反演8个分潮垂向负荷位移优于全球海潮模型，但水平方向上的反演效果稍弱。     

基于重力测点和段差变化的前郭5.8级震群研究与三维反演

利用东北地区近几年流动重力观测数据，分析前郭5.8级震群周边重力点值的时序变化特征，并基于重力段差变化对东北地区整体和局部重力变化引入能够反映异常显著性程度的指标量G和C值，最后对前郭5.8级震群前的重力场动态变化作三维密度反演。结果表明：1）前郭5.8级震群发震构造两侧的局部重力变化具有较明显的差异性，其变化趋势可作为识别构造活化或解耦运动的标志；2）显著性指标量G和C值能够为地震重力前兆的定量描述提供新思路；3）震前重力变化的反演结果显示地震发生在质量运移的过渡区域，垂向反演结果对本次地震的震源深度有一定的揭示意义。     

平衡剖面的制作流程及其地质意义

平衡剖面技术是地质思维和计算机技术的结晶,使对断层构造的研究提高到定量阶段,其依据是在垂直构造走向的剖面上,地层长度和面积(2D)或体积(3D)是均衡的。在此原理基础上利用数学手段对盆地的构造发育史进行正演和反演模拟,直观地再现地下构造的原始几何形态,迅速提供地震剖面的构造解释方案,并对解释结果进行检验(不平衡的剖面其解释一般有问题),为深刻认识构造发育史、分析油气运移及聚集规律提供依据,提高了工作效率。其结果也为盆地模拟、油藏模拟、定量计算构造伸缩量等地质研究打下了坚实的基础[1]。     

用双三次样条函数和GPS资料反演现今中国大陆构造形变场

将中国大陆现今构造变动视为一种连续的地壳变形,利用双三次样条函数模拟了近期GPS测定的大陆内部及周边地区412个测站速率,反演大陆地区自洽的构造变动速度场和应变率场.模拟结果显示:印度板块与欧亚板块的碰撞、挤压是构成中国大陆内部岩石层水平形变的主要驱动力.印度板块在东喜马拉雅构造结深深插入青藏高原,造成地壳大规模的缩短和抬升.青藏高原东南部的喜马拉雅带、拉萨和羌塘地块以及青藏高原东南边的川滇地区,内部构造活动强烈,其内部的构造变形包含地壳碎片的冲断、褶皱和侧向逃逸.大陆地壳(或岩石圈)的增厚,尤其是喜马拉雅山脉南北向的快速缩短和青藏高原东西向的缓慢拉张,大约吸收了印欧板块会聚量的85%,西藏中东地区东西向的拉张速率达到了(16±2.0)mm/a,且顺时针方向扭转明显.印度板块相对欧亚板块运动的欧拉极为(29.7°N, 19.3°E, 0.392°/Ma);华南地块相对于欧亚大陆向东(102°±7.4°)南的运动速率是(11±1.54)mm/a,华南块体相对欧亚板块运动的欧拉极为(62.25°N, 126.56°E, 0.141°/Ma);塔里木地块相对较稳定,其西部运动速度高于东部运动速度,作顺时针方向旋转.总体上讲,中国大陆运动方向为北偏东呈辐射状,从西部近南北方向的运动转向东部地区东南方向的运动,绕东喜马拉雅构造结有一顺时针方向的旋转.横穿喜马拉雅构造带及青藏内部的南北向压缩速率为(19±2.0)mm/a,横穿西天山构造带的南北向压缩平均速率为(13±1.5)mm/a,横穿东天山构造带的南北向压缩平均速率为(6.0±1.4)mm/a.阿尔金断裂带的左旋走滑速率为(6±1.2)mm/a.     

用实测波浪数据确定JONSWAP谱参数

本文给出应用实测波浪数据确定JONSWAP谱参数的方法，若Hs,Tz,Tc分别为实的有效波高，平均上跨零线周期，平均峰－峰周期并令y=Tz/Tc，本文导出求解谱参数γ值的方程f（γ，y）＝0。根据这个方程的解并应用最小二乘法，γ作为y的函数的表达式可写为（对于2.54≤γ≤15.34及1.6≤y≤1.79）。