夏日哈木矿床是青海省近几年在东昆仑地区发现的首例岩浆熔离型铜镍硫化物矿床,矿床规模达超大型,在其发现及评价过程中,化探综合找矿方法起到了最为主要的作用。在系统评价1:20万水系沉积物异常→1:5万水系沉积物地球化学特征→1:1万土壤异常特征→地质及工程验证的方法组合及应用效果基础上,研究发现:1:20万水系沉积物测量可大致圈定找矿靶区,但相应地质工作未发现与铜镍矿有关的基性—超基性岩体,引起异常的原因尚不能定性;1:5万水系沉积物测量进一步分解了1:20万水系沉积物异常,相应地质工作发现了与铜镍矿有关的基性—超基性岩体,但其含矿性需进一步确认;1:1万土壤测量进一步缩小了找矿靶区,配套完成的1:1万地质草测工作,基本查明了引起异常的原因,并详细圈定了含矿基性—超基性杂岩体,为进一步工程验证并发现矿(化)体提供了直接依据。通过研究确认其方法组合是达到最经济成本、最快发现矿体的有效方法组合,可为成矿带相似类型矿床的评价起示范作用。
针对南海北部陆坡T研究区的地质情况和天然气水合物的赋存特征,在缺少钻井资料的情况下,利用T研究区准三维地震数据和叠加速度谱资料,建立三维层速度模型,并将其转化为低频波阻抗背景模型。提取虚拟井进行精细标定,通过稀疏脉冲反演方法,最终获得研究区绝对波阻抗体和精细速度体。同时,综合利用地震、反演结果分析研究区水合物识别的最佳门槛值,并以此为依据对水合物进行定量解释。研究结果表明:① 研究区内水合物的速度在2 050~2 150 m/s之间,纵波阻抗在3.9×106 kg·m-3·m·s-1以上,与围岩速度、阻抗反差较大;② 水合物分布范围较广,在研究区中西部及东南部尤为发育,初步估算其厚度最大可达48 m左右。
鄂尔多斯盆地苏里格气田奥陶系马家沟组风化壳储层厚度薄,非均值性强,气水关系复杂,其储层发育受沉积、构造、成岩等多因素共同作用,勘探难度极大。以桃7区块马五13 段为例,结合风化壳储层发育特点及储层段地球物理响应特征,应用现代体属性分析技术对古岩溶风化壳储层进行了针对性研究,其中包括应用构造导向滤波技术提高原始地震资料信噪比,并在此基础上应用体曲率和蚂蚁追踪技术综合预测微断裂;应用地层切片及体融合技术精细刻画古岩溶沟槽。预测结果表明,现代体属性分析技术可以更充分地利用常规地震资料进行风化壳储层精细描述,减少预测多解性。多种技术方法通过实钻井验证,吻合度较高,为研究区下一步的油气勘探提供了良好的依据。
通过对雪峰山西侧地区震旦系—下古生界海相层系开展1:5万高精度重力、磁法、大地电磁测深综合物探剖面测量,结合钻井、地质资料及物性统计特征,对获取的重、磁、电资料进行联合正、反演解释,解译了断裂构造及地层空间展布特征,定量解释了地层的分布、埋深及厚度,预测了油气有利储集构造带。研究结果表明,在雪峰山西侧的复杂地质条件和人文、工业严重干扰下,重、磁、电方法组合能够最大程度地减少地质解译的多解性,较准确地查明了该区震旦系—下古生界海相层系分布及褶皱、断裂、局部构造,为该区油气远景评价发挥了重要作用。
宁东地区是宁夏规模最大的含煤区,已发现9个石炭—二叠纪和侏罗纪大煤田。笔者对宁东地区1:20万区域重力资料进行进一步处理,根据宁东含煤区石炭—二叠纪煤田与侏罗纪煤田2种类型,采用2.5D正演建立煤田重力场模型,并说明采用位场分离方法能够分离与识别含煤区;利用小波分析方法对布格重力进行分析,结果发现:小波一、二、三阶细节与已知的煤田对应较好,而滑动平均等传统的位场分离方法分离的局部异常则无明显对应关系。利用重力小波一、二、三阶细节预测了宁东地区石炭—二叠系和侏罗系含煤区的远景,充分发挥了1:20万区域重力资料在煤田勘探中的作用。
怀化地区金矿成矿研究程度低,为提高该地区研究级别,首次采用了基于MRAS的证据权重法对该地区金矿进行成矿预测。在充分收集、整理该地区各种地质、矿产、构造、物探、化探和遥感资料的基础上,建立怀化地区金矿多元地学空间数据库。在系统总结怀化地区金矿成矿规律的基础上,重建了该地区的金矿找矿模型,提取并构置了新元古代地层、断裂、航磁-60~0 nT、Au化探异常、Cu化探异常、Pb化探异常、Sb化探异常、Zn化探异常、岩浆岩、遥感推断环形构造、遥感推断线性构造、重力5~20 mGal和30~55 mGal、Au矿物重砂异常和其他矿物重砂异常等14个预测要素作为证据因子。预测网格单元划分采用2.5 km×2.5 km并进行了后验概率计算,根据后验概率计算值的分布特征将怀化地区的预测区块划分为A、B、C三个等级,并圈定了10个找矿靶区。根据找矿靶区的面积、已发现矿床(点)的数量和后验概率平均值等对找矿靶区进行了分级,最终将找矿靶区分为Ⅰ级找矿靶区4处,Ⅱ级找矿靶区3处,Ⅲ级找矿靶区3处。
对赛坝沟金矿床Ⅳ-3号矿体原生晕元素进行理论性分析后得到原生晕轴向分带特征,其前缘晕、近矿晕、尾晕分布混杂,呈现出明显的逆向分带特征,其中轴向分带序列最后出现的Au和Ni呈现良好的相关性以及前缘元素Hg的出现,指示矿体向下的延伸空间较大。进一步对其轴向地球化学参数进行对比,矿体由浅至深,地球化学参数具有"升—降—升"的变化特征,显示矿体在3 162 m中段往下还有很大的延伸性。通过对该矿体原生晕轴向分带序列及地球化学参数的研究对比,建立矿体的叠加晕理想模型,为矿区深部及其周边小型矿床的评价预测提供地球化学依据。
对赞比亚西北省卡森帕A矿区进行土壤地球化学测量,综合分析测区Cu、Co、Ni、Au、Ag、As、Pb、Zn等8个元素在土壤中分布的均匀程度、富集程度、相关关系、单元素分布特征与套合关系等,认为测区最佳找矿指示元素组合为Cu-Ni-As-Co-Zn,应用异常打分法圈定了一个Cu-Ni-As-Co-Zn综合异常Ht-1,经过野外地质调查,认为该异常为矿致异常;后期通过激电及高精度磁测,发现高磁、高极化、低电阻异常与化探异常吻合;经两个钻孔验证,见工业矿体,最终探获了一个中型铜多金属矿床。
花岗岩型铀矿床一般产在岩体内部或岩体周边,因此岩体发育情况的探测对铀矿勘查工作的进一步开展具有重要意义。简要介绍了可控源音频大地电磁法在卫境地区的实际应用,结合地质资料和岩石电阻率特征,利用可控源音频大地电磁测量资料,查明了花岗岩体的厚度以及花岗岩和变质岩接触带的深部发育特征,经钻孔验证,效果较好。综合分析认为:可控源音频大地电磁测深法有较好的垂向分辨力,能够较准确探测岩体的厚度、断裂发育的位置及向深部延伸情况、确定破碎带及裂隙带和脉体的空间展布特征,是花岗岩型铀矿勘查的有效技术手段。
湘赣边区是一个多期复合造山活动和构造叠加的复杂构造域,地质构造发育,也是华南地区有色金属、多金属及铀矿的主要产地,但受诸多因素的影响限制,地质构造的研究长期以来缺乏系统性、全局性。笔者利用最新实测的高精度航空物探(磁/放)数据资料,结合区重资料,系统的对湘赣相邻边区的地质构造进行了圈划,对深大断裂带的分布及控岩特征、钦杭结合带的地球物理特征及基本展布情况进行了论述。研究成果对本区进一步开展基础地质研究及矿产资源远景预测具有一定的参考价值。
乌尔吉地区地处中亚造山带东段,属大兴安岭南段铅锌银铜锡铁成矿带。区内二叠系大石寨组为一套变质火山岩夹碎屑沉积岩组合,火山岩具有典型的不对称双峰式特征,基性峰强。岩石地球化学研究表明,中基性、基性岩类属准铝质钙碱性系列和低钾(拉斑)系列,中性岩类总体属钾玄岩系列,中酸性岩属过铝质钙碱性岩系列,岩石为海西中晚期大陆边缘裂谷盆地的岛弧环境下的产物。大石寨组火山岩总体上富集Ag、Sn、Pb、Zn、Cu等成矿元素,是本区Pb、Ag、Sn、Zn等成矿元素的矿源层之一,成矿潜力较大。结合前人研究认为,大石寨组火山岩是本区寻找火山热液型锡、铅、锌、铜等多金属矿床的有利场所。
借助三维自然电位电源密度反演成像方法,分析南京中山植物园试验场地的野外电位观测资料,研究包气带的介质结构信息。研究结果证明,电源密度反演成像结果与探地雷达获得的地下介质结构特征具有很好的吻合性,这表明三维自然电位电源密度成像可以应用于探测包气带介质结构。这为研究地下空间结构提供了一种新的地球物理方法。
将等效导电平面快速算法与合成孔径相关叠加技术相结合,实现了航空瞬变电磁快速成像解释算法。为了进一步获得更好的成像效果与更优的分辨率,对合成孔径算法中的关键技术进行了研究。首先给出最佳孔径尺寸的求取方法,并通过多个模型进一步验证孔径尺寸对于成像效果的显著影响;接着利用大量数值模拟结果分析孔径大小与异常体埋深、围岩与异常体电阻率比值等因素之间的关系,总结规律利用数值拟合方法求取其函数表达式,并根据函数表达式直接确定最佳孔径尺寸的大小;最后,利用多个模型示例证实合成孔径成像方法能够有效圈定异常范围、提高横向与纵向的目标分辨率,并且该方法不受海量航空数据的影响,计算速度快、可用于实时成像。以上研究成果进一步提高了航空瞬变电磁法的解释精度,为丰富航空瞬变电磁解释方法、提高勘探效率提供了理论依据。
Radon变换是地震数据处理中广泛应用的一种技术方法,但实际应用中存在分辨率及运算效率低的问题。笔者在相移双曲线的基础上,推导了时不变双曲积分路径,提出了基于贝叶斯原理的稀疏双曲Radon变换,给出离散变换的参数选取准则,并且为保证不同反射层的反射波具有最佳收敛,在变换中引入了动态时窗。频域稀疏双曲Radon变换不仅具有较高的分辨率,而且可以在频域实现,计算效率大幅提高。模型数据测试和实际资料应用结果表明,频域稀疏双曲Radon变换可以有效去除规则干扰波、压制随机噪声,提高数据信噪比,且具有较好的保真度和保幅性。
复合电偶源瞬变电磁张量测量能够有效解决单一电偶源在复杂三维电性结构体的瞬变场响应失真的问题。采用复合电偶源瞬变电磁张量的坐标不变量P2研究均匀大地介质和一维典型层状介质,计算结果表明,坐标不变量P2参数不论在早期或晚期都能准确反映一维介质的电性规律。对于复杂的三维地质体,采用坐标不变量P2以及张量电阻率的极大值ρmax和极小值ρmin组成的旋转椭圆表征三维电性结构的特征,用低阻长方体模型进行正演模拟,计算结果表明:不同时间窗口坐标不变量P2的响应以及极大值ρmax和极小值ρmin的旋转椭圆,不但可以准确反映出地下三维结构的变化情况,而且能够直观地勾画出异常区域边界的变化范围,这为后期的资料解译提供了理论指导。
针对页岩脆性各向异性特征的研究相对薄弱,亟待发展与之相应的页岩储层脆性研究。为实现页岩脆性指数的表征,基于多尺度建模思想,引入各向异性元依次建立纳米级多孔黏土颗粒、亚微米级富有机质黏土颗粒和微米级干岩石骨架,并综合考虑吸附气、束缚水等影响因素,改进了常规基于DEM理论及Gassmann理论的建模流程,建立了页岩储层等效介质模型并实现了弹性参数的表征。然后,引入方位角的影响,依据Bond变换推导得到了杨氏模量、泊松比等脆性指数各向异性表达式,实现了以等效介质模型弹性模量为参数的脆性参数表征。最后,总结了页岩储层脆性的各向异性特征以及页岩储层脆性对矿物组分的敏感性,发现平行于储层走向且垂直于层理方向上的脆性是垂直于储层走向且与层理45°夹角时的2.5倍;同时,脆性对石英的敏感强度约为黏土的2倍。该特征研究可以为储层脆性预测及水力压裂方案设定提供参考依据。
通过搜集柴北缘1:20万区域地球化学数据,运用分形(多重分形)方法来确定荒漠戈壁景观区受风成沙影响的地球化学异常特征,并利用C-A和S-A模型分别探讨了景观地区Au的地球化学异常特征。研究表明,C-A模型在确定异常下限时有优势,但也存在局限性;S-A模型能够很好地排除背景场和干扰因素,缩小异常面积,突出弱小异常。基于该方法提取的异常元素区域能够较好地识别出绝大多数已知矿床,同时指示了值得作进一步异常查证的空白成矿元素异常区域,为勘查工作提供参考,值得开展下一步的找矿工作。
地电化学泡塑载体分析测试中一直缺乏有效的分析质量监控方法,影响了方法技术的应用推广。针对应用较广泛的泡塑预处理方法——灰化法,实验了在空白泡塑中加入标准溶液的方法,对分析质量进行监控。结果表明,使用这种方法应该可以解决Au的分析质量监控问题,实验的分析准确度甚至可满足多目标区域地球化学调查中土壤样品分析质量监控要求;而针对Cu、Ni、Pb、Zn等,在灰化过程中存在不同程度的元素损失,暂时不适合使用该方法。
在地震勘探中,强振幅、低频率的面波对于后续地震资料的处理和解释会产生重要影响。基于最小平方滤波的面波压制方法是一种新的面波压制方法,它与一般的频率域滤波方法不同。该方法的主要原理是利用扫描信号与滤波因子的褶积来预测面波,其中滤波因子的求取需要建立一个目标泛函。当目标泛函取得最小值时可从地震数据中求取滤波因子进而预测出面波,再将预测出的面波作为输入进行第二次面波预测,从地震记录中减去提取的面波,就得到去除面波后的地震数据。采用非线性的扫描信号提取面波,并对扫描信号起止频率进行了分析;通过对合成地震记录与煤田实际地震资料的处理,充分证明该方法能有效地压制面波。
在浅水区,由于水深较浅,缺少近偏移距的海底反射信号,采用常规SRME方法和预测反褶积不能有效地衰减海底相关多次波。笔者采用DWD+SRME组合法来衰减海底相关多次波,其中DWD方法用来衰减与海底相关的短周期多次波,SRME方法压制自由表面相关的长周期多次波,另外运用Radon变换去除层间多次波。结果表明,运用DWD+SRME组合法压制海底相关多次波效果较好,并且没有损害浅层的有效反射信号,采用Radon变换使层间多次波得到了良好的压制,深部有效反射信号逐步体现。因此,对于浅水区低信噪比地震资料,采用DWD+SRME组合方法和Radon域去多次,是比较可行的去多次波方法。
地震勘探勘查区往往位于偏远地区,已知地形资料较少。Google Earth(简称GE)采用的数字高程模型数据为SRTM,精度较高。笔者在批量提取GE高程数据的基础上,分别对山区和平原地区的GE高程数据的精度进行了分析。结果表明,无论是在山区还是平原地区,GE高程数据均与地表实际高程基本一致;在平原地区的精度高于山区;GE高程数据可以满足等高线生成和观测系统参数分析的需要。
为了适应地形起伏的实际地质情况,开展带地形的最小二乘二维反演研究。鉴于大地电磁(MT)反演的不适定问题,引入Tikhonov的正则化方法,从而获得关于总目标函数的方程,利用光滑约束最小二乘法求解总目标函数方程。由于正则化因子值与反演精度以及稳定性相关,采用主动约束平衡方法获取最优化的正则化因子,以确保反演精度和稳定性都达到最佳。与此同时,利用电磁场互易定理以节省反演迭代过程求解雅可比矩阵的计算时间。构建了若干地质构造模型进行试算,分别讨论TE、TM模式以及二者联合模式的反演结果,并与前人研究工作对比以说明本文方法的反演效果。
目前,可控源音频大地电磁法测深资料反演主要基于各向同性介质,在复杂的地电体中,由于岩石裂隙或矿物形状等原因使得在各个方向的电阻率不同。本文主要探讨更为符合实际地质体情况的方位非各向同性层状介质在正交水平电偶源激发条件下的可控源音频大地电磁数据的反演。运用Occam反演方法,实现了可控源音频大地电磁法方位非各向同性层状介质反演。构建理论地电模型,通过理论模型合成数据的反演算例检验了所实现的可控源音频大地电磁法方位非各向同性层状介质反演算法的有效性和稳定性。同时,不同数据的反演结果表明:方位非各向同性介质的反演应同时使用视电阻率ρsxy、ρsyx和相位φxy、φyx数据。
在海洋磁测数据日变改正过程中,常会因各种因素限制,导致日变数据缺失,这时就需要结合多台站数据推算工区位置处的日变数据。通过借鉴前人研究成果及考察相关分析结果,可以近似认为台站间日变数据具有一定的线性关系,并引入统计学中常用的回归分析法拟合出某一位置日变数据的推算公式,进而得到推算数据。选取南海区域实测数据验证可知该方法具有一定的可靠性,可以解决日变数据缺失等问题,从而提升海洋磁测数据的精度,是一种改进的磁测数据日变改正方法。
通过分析测定渤海湾100个站位表层沉积物中重金属元素的含量,探讨其分布特征及富集状况,并通过地累积指数法和潜在生态危害指数法对渤海湾西部岐口凹陷海域底质环境进行了评价。结果表明:渤海湾西部海域表层沉积物中As、Hg、Cu、Pb、Cd的平均含量分别为12.9×10-6、0.031×10-6、29.5×10-6、27.4×10-6、0.159×10-6。结合5种元素的地球化学图及地累积指数法分析表明,研究区表层沉积物整体上未受到As、Hg、Cu、Pb、Cd的污染,而在高值区Pb和Cu达到轻度污染的水平;潜在生态危害指数法的评价结果指明,重金属元素的潜在危害程度顺序是Hg >Cd >As >Pb >Cu,研究区重金属的潜在生态危害轻微,只有在高值区,Hg和Cd才达到中等程度的生态危害。
以海南省东部某典型地区市县级农业地质普查为例,介绍了基于MapGIS的Lable点法样点布置方法。该方法首先将海南省第二次土地利用调查的26个土地利用小类归并成水田、旱地、果园、其他园地、林地、建设用地、滩涂、交通用地、水面、其他土地共10大类,图斑总数由9 932个归并为4 257个;然后将归并后的图斑按照面积大小划分为代表性依次降低的A、B、C三个级别;样点布置过程中,A级图斑采用传统500 m×500 m网格化方法布置样点,B级图斑采用基于MapGIS的Lable点法布置样点,C级图斑不布置样点。软件自动布置样点后,经必要的人工调整,在174 km2的工作区共布置1 610个采样点,样点密度达9.25个/km2,符合农业地质调查相关规范要求。方法有效性验证显示:与传统单一的333 m×333 m网格化样点布置方法相比,该方法布置的样点在可采性、代表性、准确性、合理性四方面都有显著优势,值得推广应用。
当井—地直流电法实现了高清晰地质扫描,电法勘探技术的发展就到了一个全新的定量解释阶段,不仅拓展了地质勘查更广阔的领域,而且能够解决一些过去无法解决的地质勘查难题。以某热力站工程为例,在建筑施工阶段发生岩溶地质坍塌时,成功使用井—地直流电法实施了高精度地质勘查,在短时间内,低成本、高效率地完成了复杂岩溶地质勘查难题,探明了岩溶地质发育带的大小、延伸方向和埋藏状态,确定了岩溶地质特征,为解决岩溶地质灾难隐患,设计工程施工方案提供了可靠的地质资料依据。
将超高密度电阻率法用于某水库坝址区隐伏断层探测。介绍了测区地质及测线布置情况,采用汉克尔变换方法模拟了层状介质模型视电阻率响应,通过反演计算结合地质钻孔资料,确定了基岩埋深及区域断层的走向。断层破碎带在电阻率剖面中呈现明显斜立条带状低阻异常,其产状与地质及钻孔资料一致,验证了超高密度电法探测的有效性。
在研究可控源音频大地电磁测深与NLCG反演方法基本原理的基础上,分析了运用钻探、综合测井、地质、地震、高密度电法等先验信息的电性约束反演方法,并结合宝兰客专某隧道的实测CSAMT资料进行了电性约束反演前、后的对比分析,证明了约束反演结果在断层、风化层、岩性分界、地层厚度确定等方面的准确性更高。
在鱼卡煤田1171工作面巷道采掘过程中,前方出水量明显增大。为避免水害发生,在整个工作面开展了地面与井下相结合的瞬变电磁探测工作,利用地面瞬变电磁测量从平面上控制水害分布范围,利用井下探测方法对富水区域进行圈定,指导下一步的抽排水工作。对比地面与井下瞬变电磁结果,两者对于富水区域结果是相同,联合探测降低了物探资料多解性问题,大大提高了勘探精度。