衰减夹层GPR模拟的时频域全局反射误差

全局反射误差分析是深入研究探地雷达（GPR）吸收边界条件吸收效率的有力工具.基于常规完全匹配层（PML）的标准交错网格有限差分算法必须满足严格的CFL条件限制,即在单位时间步长内,不容许电磁波传播的距离超过单元网格尺寸.为了提高主区域所有网格节点的计算效率,并有效地吸收传播后期出现的低频隐失波,提出基于非分裂递归卷积完全匹配层（UCPML）的旋转交错网格（RSG）GPR正演算法.该算法采用不同的网格交错策略,并在边界条件中引入了吸收低频隐失波的自由可变因子,使得该算法允许选取较大的时间步长,提高了计算效率,并且实现了对低频隐失波的高效吸收.本文首先给出了RSG差分格式,推导了基于UCPML的RSG差分更新方程,探讨了数值色散的稳定性条件,然后以绕射现象严重的衰减夹层数值模拟为例,分别从波场快照、单道波记录、时间域反射误差（TDRE）、频率域反射误差（FDRE）四个方面分析了UCPML与常规PML的全局反射误差,说明了基于UCPML和RSG的GPR正演算法能更有效地吸收低频隐失波.

     

利用地震折射和反射波资料研究银川盆地浅部结构和隐伏断裂

银川盆地是华北克拉通西部构造活动较为强烈的一个新生代断陷盆地.为了研究银川盆地的地壳浅部结构和活动断裂特征,我们利用2014年在银川盆地完成的深地震反射剖面数据,采用初至波层析成像方法得到了银川盆地高精度的基底P波速度结构和构造形态;考虑到仅根据速度结构剖面还难以确定断裂的准确位置、断层上断点埋深、断层的近地表构造组合样式等特征,研究中还采用浅层地震反射波勘探方法对银川盆地内的隐伏断裂和1739年平罗8.0级地震的地表破裂带浅部结构进行了高分辨率成像.研究结果表明：银川盆地与两侧地块的浅层P波速度结构和沉积盖层厚度差异较大,银川盆地总体呈现出明显的低速结构特征,盆地基底面起伏变化较大,基底最深处位于芦花台断裂和银川断裂之间的银川市下方,其深度约为7000~7200 m;贺兰山隆起区显示为明显的高速特征,地表出露中-古生代基岩地层,缺失新生代地层;鄂尔多斯地块西缘的浅层P波速度明显高于银川盆地,基底埋深相对较浅,推测其新生界地层厚度小于2500 m.浅层地震反射剖面揭示的地层反射界面形态和断裂的浅部构造特征非常清楚,黄河断裂、贺兰山东麓断裂、银川断裂和芦花台断裂不仅是错断盆地基底的断裂,而且还是第四纪以来的隐伏活动断裂,这些断裂的交替活动形成了"堑中堑"的盆地结构,并对银川盆地的形成、盆地内的新生代地层厚度和第四纪沉降中心具有重要的控制作用;在近地表这些断裂表现为由2~3条断层组成的"Y字形"断裂构造,且主断裂的最新活动可追踪至晚更新世末期或全新世,是构造继承性活动的结果.本文的研究结果不仅可为进一步分析银川盆地的基底结构、隐伏断裂特征和活动构造研究等提供新的地震学证据,而且还可为该区城市规划中避让活动断层提供科学依据.     

四川盆地深部地壳结构——深地震反射剖面探测

四川盆地位于扬子地块的西北部,被褶皱构造带所围绕,受周缘构造带的侧向挤压作用,盆地卷入了多期次和多边界的构造变形,为开展盆山耦合作用及多边界、多期次构造叠加与复合关系的研究提供了不可多得的理想野外实验室.为揭示四川盆地地壳结构,本文通过对3条不同时间采集的深地震反射剖面数据进行拼接联线处理,获得跨越四川盆地的330 km深地震反射偏移成果剖面,揭示了四川盆地地壳上地幔细结构：沉积层从西北向东南逐渐变薄,在龙门山前沉积层厚度超过15 km,在华蓥山下沉积层减薄到~8 km,且褶皱变形形成华蓥山薄皮褶皱冲断带;莫霍面出现在13~15 s（双程走时）,埋深约40~45 km;并发现从下地壳延伸至地幔的东南向的倾斜反射,从13 s向下延伸至18 s,结合四川盆地及其周边地区其他地球物理和地球化学花岗岩同位素年龄等资料,我们认为这些倾斜反射层是扬子克拉通地台西北缘发生的新元古代俯冲的遗迹.

     

临南洼陷带滩坝砂岩储层分布

临南洼陷带沙四段滨浅湖滩坝相成因的砂体广泛发育,近年来已成为该区十分重要的勘探后备阵地。但该区火成岩发育,地震资料品质较差,滩坝砂的分布发育不清,不利于勘探部署。在深化研究沙四上亚段古地貌、沉积特征的基础上,运用频谱分解等技术提高地震资料分辩率,综合利用了正演及属性分析等技术,对滩坝砂岩的发育范围进行了预测,取得了良好的效果,为该层系油藏的勘探提供了可靠的依据。     

非对称圆弧片状声源在井旁地层中产生的声场

提出了一种可用于随钻三维反射声波成像测井的圆弧片状声源,推导了该声源声学性质在波数-频率域内的数学描述,并利用实轴积分的方法对该声源在井旁地层中产生的声场进行求解,考察了声源的线度、频率等因素对该声源向地层中辐射的声场的影响.研究结果表明,特定尺寸和特定频率下的圆弧片状声源产生的纵波场的水平和垂直指向性图中仅具有一个明显的主瓣,且主瓣三分贝角宽较窄,方位分辨率较高,指向性良好,适用于随钻三维反射成像测井;SV波场和SH波场的指向性图中存在两个或者多个角瓣,SV波场在垂直于井轴方向上辐射的能量为0,另外在井中发射和接收的反射SH波场和SV波场会相互干扰,所以难以利用SV波场和SH波场进行三维反射声波成像;声源的线度和频率等因素对其辐射声场影响较大,考虑到激励效率和方位分辨率等因素,对于本文描述的井孔模型,选择圆周角在75°和90°之间, 主频为12 kHz左右的声源是合适的.     

基于模糊集分析的土壤质量指标高光谱反演

通过对横山县84个土样350~2 500 nm波段的光谱曲线分析,进行土壤质量12项指标的光谱反演。为获得有效光谱反演指标,利用单相关分析方法计算土壤质量指标与光谱反射率的相关系数并绘制相关系数曲线,根据极大相关性选择最佳波段作为光谱反演指标;剔除异常样本后,利用模糊识别理论建立土壤质量指标反演模型,通过优化得到模型的最佳参数。土壤质量12项指标光谱反演模型的平均检验误差均小于15%,说明模型的精度较高。     

反射地震走时层析成像中的大型稀疏矩阵压缩存储和求解

反射地震走时层析成像是一种精度较高的速度求取方法,最终可归结为线性方程组的求解.方程组具有很大的维数,常规解法需要很大的存储量和计算量.本文考虑到当投影函数取为走时残差,图象函数取为慢度残差时,灵敏度矩阵中的元素表示射线经过网格的长度的特殊物理意义,采用行索引的压缩存储方式,在射线追踪正演模拟过程中直接压缩存储灵敏度矩阵,在层析反演过程中利用压缩后的矩阵进行求解,大大降低了存储量和计算量.     

沂沭断裂带深反射地震资料处理

地震时间剖面可以对地质目标进行精细刻画,但由于变质岩、侵入岩地区存在地质体不规则、地质界线间波阻抗差异不明显、断裂构造倾角高陡等特征,导致地震原始资料具有波场复杂、干扰波发育、资料信噪比低等特点.为获得客观反应实际地质情况的沂沭断裂带深反射地震剖面,通过对研究区原始地震资料细致分析,针对各种干扰波的不同特性特征,分别采...     

帕米尔东北侧地壳结构研究

1998年在帕米尔东北侧伽师及其周边地区完成了两条深地震宽角反射/折射剖面. 结果表明,西昆仑、塔里木和天山在地壳速度结构、构造特征上显示出较大差异. 塔里木块体具有稳定地块的地壳结构特征,地壳平均速度较高(6.5km/s). 向南进入西昆仑,地壳明显增厚,厚度可达0km左右,且地壳平均速度偏低(6.0-6.2km/s),偏低的地壳平均速度主要来源于相对低速度的下地壳结构,反映了西昆仑褶皱系下地壳介质的特征. 向北进入天山后,地壳同样明显增厚,但增厚的程度低于西昆仑下,约为50-55km. 天山地壳同样具有明显低的平均速度(6.2km/s),显示了天山地壳相对"软"的特征,但天山地壳偏低的平均速度来源于广泛分布于中地壳的低速度层和速度偏低的下地壳. 在印度块体向北强烈推挤的作用下,该区地壳遭受强烈的不均匀变形,塔里木块体向南插入西昆仑下,向北插入天山下,形成了该区强烈地震频繁发生的深部构造环境.     

标定大陆科学钻探孔区地震反射体

在中国大陆科学钻探孔区(江苏省东海县南部)进行了系统的地球物理调查,包括二维地震测网和专门的地震剖面,大地电磁法和位场方法等.地震调查表明,在超高压变质岩出露区上地壳充满了反射体,包括倾斜反射体与上拱的反射弧. 本文介绍大陆科学钻探先行研究中地震调查的成果.根据大陆科学钻探预先导孔的岩芯和测井资料、井旁VSP和数值模拟结果证实,高波速的榴辉岩体、破碎断裂带和大型韧性剪切带都可引起倾斜的地震反射,而上拱的弧形反射体则是由近似直立的榴辉岩体和其中的破碎带的综合反映.由于地壳深部广泛分布着经受变质的岩石,上述研究结果对标定地壳中的反射地震信号具有一定意义.