基于重力异常识别断裂的位错窗口相关法及其在大兴安岭与邻区的应用.pdf

大庆石油地质与开发 Petroleum Geology Oilfield Development in Daqing2023 年 6 月 第 42 卷第 3 期June，2023Vol.42 No.3DOI：10.19597/J.ISSN.1000-3754.202202024基于重力异常识别断裂的位错窗口相关法及其在大兴安岭与邻区的应用李希元1，2 胡望水1，2（1.长江大学地球科学学院，湖北 武汉430100；2.长江大学油气资源与勘探技术教育部重点实验室，湖北 武汉430100）摘要：断裂的划分及位置的确定对区域地质研究具有重要的作用。通过分析重力的方向导数、小子域滤波等传统划分断裂的方法，提出突出重力线状异常或梯度带拐点的划分断裂、确定断裂位置的位错窗口相关法。结果表明：与方向导数、小子域滤波等方法相比，位错窗口相关法能够有效提高强背景重力异常中断裂线状异常与梯度带的信息；以单一断层模型、剪切断层模型、正常埋深地质体与控凹边界断层模型、较深埋深地质体与控凹边界断层模型验证了位错窗口相关法具有准确识别断裂重力异常信息、确定断裂位置的作用；应用位错窗口相关法能够较好地指示大兴安岭及邻区的德尔布干断裂、大兴安岭莫霍面陡变带、嫩江断裂、孙吴双辽断裂和赤峰开源断裂的位置。研究成果为深化大兴安岭与邻区深部地质规律的认识提供了新的依据。关键词：位错窗口相关法；重力异常；线状异常；梯度带；断裂识别；大兴安岭与邻区中图分类号：P631 文献标识码：A 文章编号：1000-3754（2023）03-0122-10Dislocation window correlation method based on gravity anomaly identifying fault and application in Greater Hinggan Mountains and its adjacent areasLI Xiyuan1，2，HU Wangshui1，2（1.School of Earth Science，Yangtze University，Wuhan 430100，China；2.MOE Key Laboratory of Exploration Technologies for Oil and Gas Resources，Yangtze University，Wuhan 430100，China）Abstract：Division and location of faults play an important role in regional geological research.Based on analysis of gravity directional derivative and small domain filtering，a gravity dislocation window correlation method is proposed to classify faults with prominent gravity linear anomalies or inflection points in gradient zones and to determine location of faults.The results show that dislocation window correlation method can effectively improve the information of fracture linear anomaly and gradient belts in strong background gravity anomaly compared with traditional methods such as directional derivative and small domain filtering.Single fault model，shear fault model，normal buried geological body and sag-controlling boundary fault model，deeper buried geological body and sag-controlling boundary fault model verify that dislocation window correlation method can accurately identify gravity anomaly information and determine fault location.The location of Delbugan fault adjacent to Greater Hinggan Mountains and adjacent areas，Moho surface abrupt change zone of Greater Hinggan Mountains，Nenjiang fault，Sunwu-Shuangliao fault and Chifeng-Kaiyuan fault are well indicated by using dislocation window correlation method.The re收稿日期：2022-02-24 改回日期：2022-06-10基金项目：国家重点研发计划项目“北方东部复合造山成矿系统深部结构与成矿过程”（2017YFC0601300）。第一作者：李希元，男，1993年生，在读博士，从事盆地构造研究。E-mail：通信作者：胡望水，男，1963年生，教授，博士生导师，从事盆地构造教学与研究。E-mail：第 42 卷 第 3 期李希元 等：基于重力异常识别断裂的位错窗口相关法及其在大兴安岭与邻区的应用search provides new basis for deepening understanding of deep geological laws in Greater Hinggan Mountains and its adjacent areas.Key words：dislocation window correlation method；gravity anomaly；linear anomaly；gradient belt；fault identification；Greater Hinggan Mountains and adjacent areas0引言断裂是地壳中广泛发育的一种构造形态及构造方式。在一定程度上，地球的构造演化史就是断裂发生、发展的演化史1。大多数大、中型断裂构造直接控制了矿产的分布与成矿类型24。断裂研究对于地球的构造演化、区域大地构造、油气及矿产资源分布规律的分析都具有十分重要的意义510。断裂使地质体在三维空间发生位移和错断，形成地层密度界面的陡变带。断裂的规模越大，界面间的密度差异越显著，界面陡变带规模越大，重力异常梯级带越明显。重力线状异常或梯度带是断裂最基本的重力异常特征。资源勘探实践表明重力异常的识别是研究区域性断裂十分有效的手段1112。利用断裂在重力异常上大多表现为线状异常或梯度带的特征，前人研究了应用重力划分断裂与确定断裂位置的方法，最具代表性的是方向导数13、水平梯度矢量模1415、小子域滤波163 种，但这些方法均存在各自的不足。方向导数方法会使异常发生畸变，在突出线状与梯度带的同时，也使得不具一定走向的异常变成类似断裂的线状与梯度带异常；水平梯度矢量模方法克服了方向导数方法的不足，但由于导数的高通放大特性，增加了大量非线状异常信息的干扰；小子域滤波能够较好地突出断裂重力线状异常与梯度带的信息，不过要想取得良好的地质效果，必须要进行多次迭代，在迭代的过程中存在边缘损失。此外，3种方法都不能有效地提高与增强隐藏在重力中的断裂线状异常与梯度带信息。本文通过分析断裂在重力异常上表现的特征，提出利用重力异常划分断裂与确定断裂位置的方法位错窗口相关法，通过对比位错窗口相关法与方向导数、水平梯度矢量模、小子域滤波3种方法在利用重力异常划分断裂与确定断裂位置上的优势，并应用位错窗口相关法较好地指示了大兴安岭及邻区的德尔布干断裂、大兴安岭莫霍面陡变带、嫩江断裂、孙吴双辽断裂和赤峰开源断裂的位置，展示了位错窗口相关法在断裂构造的重力异常研究方面的良好前景。1断裂在重力异常上的特征基于断裂在重力异常上的 6 种特征1进行断裂的划分。（1）存在明显差异的异常分界线（图1（a）。重力异常具有明显差异的分界线，沿着一定走向分布的正、负异常之间存在断裂。由于断裂导致断裂两盘的地层发生差异性的沉降或抬升，使得同一深度的两盘相邻地层存在较大的密度差形成的线状重力异常。（2）异常等值线同形扭曲（图 1（b）。在地质体被断层错断的位置产生扭曲的重力异常。（3）存在串珠状异常（图 1（c）。反映沿断裂有岩浆侵入或密度明显不同于断裂内其他物质的较大块体。（4）沿着一定方向延展的梯度带（图 1（d）。由于两侧重力场特征不同，反映深大的断裂。断裂两侧发生较大规模的升降，使得断裂两盘所处的异常D oc umentPrPocsCD oo Pol:C yaMD oa umgYwPoD o Bk?r?D oe PoC?、?D of?l:CraM?Po图1断裂在重力异常上的特征Fig.1 Features of faults in gravity anomaly1232023 年大庆石油地质与开发差异大，下降盘为宽缓的低值异常，靠近断裂的重力异常呈现线状异常特征，在接近上升盘的部位，异常水平梯度较大。（5）异常线突然变窄、变宽（图 1（e）。平面投影呈现闭合的地质体在断裂两侧发生升降，产生重力异常突然变窄、变宽的特征。（6）梯度带等值线的扭曲（图 1（f）。表现为梯度带等值线扭曲、间距发生变化或等值线突然错动，反映大规模断裂被后期次级断裂所切割。从断裂在重力异常上所表现的特征可以看出，重力线状异常或梯度带是断裂存在的基本特征。断裂能够产生以上重力异常特征，但这些特征并不是断裂特有的，不同密度的岩性界面也能产生类似的特征，因此，在应用重力异常识别断裂时，应结合地质和物探等资料进行综合判别。2位错窗口相关法2.1原理基于重力异常的位错窗口相关法的基本原理是利用断裂在重力异常上表现为线状异常或梯度带的特征，沿着2个大小适当的相邻窗口在同一线状异常或梯度带所取得的重力异常数据进行位错窗口相关的计算，在线状异常或梯度带异常的拐点处会出现2个相邻窗口对应异常相关系数的最大值。重力线状异常或梯度带在异常的拐点处的邻域内，可以用1个重力异常g(x，y)=Ax+By+C的平面来代替，系数A、B、C能够很好地拟合任意方向的线状异常或梯度带。图 2 所示的计算点“O”处的主窗口重力异常的计算公式为gOW(xoi，yoi)=Axoi+Byoi+C（1）式 中：gOW 主 窗 口 重 力 异 常，mGal；xoi、yoi 主 窗 口 第i点 的 横 坐 标、纵 坐 标，i=1，2，N；A、B、C平面方程系数。在主窗口的邻域取与主窗口等尺寸的副窗口的重力异常计算公式为gMW(xmi，ymi)=Axmi+Bymi+C（2）式 中：gMW 副 窗 口 重 力 异 常，mGal；xmi、ymi 副 窗 口 第i点 的 横 坐 标、纵 坐 标，i=1，2，N。主、副窗口数据点的坐标具有xmi=xoi+x和ymi=yoi+y的对应关系。因此gMW(xmi，ymi)=Axoi+Byoi+Ax+By+C（3）式中 x、y实际重力异常网格间距。设：gO(i)=Axoi+Byoi+C(i=1，2，N)，gM(i)=Axoi+Byoi+Ax+By+C(i=1，2，N)，则重力线状异常或梯度带拐点处主、副窗口所对应数据的相关系数为R=k=1NgO(k)gM(k)-Ng Og M()k=1Ng2O(k)-Ng 2O()k=1Ng2M(k)-Ng 2M（4）式中：R主、副窗口重力异常的相关系数；g O、g M主、副窗口重力异常的平均值，mGal；n主、副窗口点数；N主、副窗口重力异常点数，N=n n。将主、副窗口的N个重力异常数据代入式（4）中得到主、副窗口重力异常的相关系数，当R=1.0时，为在具有线状重力异常或梯度带的拐点邻域内取得的最大相关系数，通常情况下，窗口宽度n取3或5时能获得理想的处理效果。在实际资料处理中，线状异常或梯度带在拐点邻域内也不可能完全在拟合的平面上，因此，线性相关系数(R)不可能为 1.0，但确实在线状异常或梯度带异常的方向上能够取得相关系数的极大值，这个极大值的延展方向就代表了线状异常或是梯度带异常所代表断裂位置的拐点的延展方向。利用该理论方法，不仅能够窄化异常突出断裂图2位错窗口相关法原理示意Fig.2 Schematic diagram of principle of dislocation window correlation method124第 42 卷 第 3 期李希元 等：基于重力异常识别断裂的位错窗口相关法及其在大兴安岭与邻区的应用信息，而且还能够有效地增强隐藏在重力异常中的断裂线状异常信息。2.2计算步骤在实际资料处理时，计算位错窗口相关可分为4个步骤。（1）以重力网格计算重力异常的位错窗口相关（图 2）。在计算网格节点 O 处，获取以 O 为中心、窗口宽度为 3 或 5 的主窗口重力异常，再分别以M1、M2、M3、M4、M5、M6、M7、M8 为副窗口中心点依次获取与主窗口相同数量的重力异常并分别计算相关系数(R)，取绝对值。（2）求取8个位错窗口与计算网格点主窗口所对应重力异常相关系数的最大值，作为 O 点的位错窗口相关的结果。这8个位错窗口覆盖了网格点O 的 N、S、E、W、NW、SE、NE、SW 共 8 个方向，便于追踪重力拐点的延展方向。（3）移动到下 1 个网格点，重复步骤（1）和（2），直至重力异常覆盖区域的所有网格点，完成重力异常位错窗口相关的计算。（4）将计算结果绘制成可用于划分断裂的图件，彩色浮雕图的效果较好。2.3模型验证为展示位错窗口相关法突出断裂重力线状异常特征的作用，设计了单一断层、剪切断层以及具有地质体与控凹边界断层的 3 种模型，利用棱柱体17为计算单元分别正演了模型的理论重力异常，对模型的重力异常进行了位错窗口相关计算。2.3.1单一断层模型单一断层为断层的基本样式（图 3（a），断层 沿 走 向 产 生 标 准 的 重 力 线 状 异 常 带（图 3（b），经过位错窗口相关处理清楚地显示了断层位于理论异常的拐点位置（图 3（c），且表示断裂位置的位错窗口相关具有显著的窄化特征，使断裂位置显示得更加准确可靠。2.3.2剪切断层模型剪切断层是另一种常见的断裂构造组合（图 4（a），常形成挤压应力下的共轭断裂。依据所设计的模型正演重力异常（图 4（b），并计算了其位错窗口相关结果（图4（c）。位错窗口相关结果简洁明了地显示了重力异常中剪切断层的信息，较好地再现了剪切断层的特征。D oc umentProcssCoclD or rocs:CyD oa rocsMagYlwBk?222426283032图3单一断层模型的重力位错窗口相关效果Fig.3 Effect of gravity dislocation window correlation in single fault model图4剪切断层模型的重力位错窗口相关效果Fig.4 Effect of gravity dislocation window correlation in shear fault model1252023 年大庆石油地质与开发2.3.3含有地质体与控凹边界断层的模型2.3.3.1正常埋深地质体与控凹边界断层模型在沉积盆地中，常见具有边界断层的断陷或凹陷。为分析位错窗口相关在含有地质体复杂构造中的应用情况，设计了正常埋深条件下含有地质体与控凹边界断层的模型（图 5（a），正演该结构模型的重力异常（图 5（b），并对重力异常进行了重力位错窗口相关处理（图5（c）。从图 5（c）可以看到，弯曲的边界断层、横切凹陷的断层及局部地质体的边界清晰，重力异常的位错窗口相关对于重力梯度异常具有良好的刻画能力，不仅能够突出断层产生的重力线状异常，而且还能够对具有密集带特征的地质体边界重力异常进行表征。2.3.3.2较深埋深地质体与控凹边界断层模型埋藏较深的地质体或深部断层，在地表产生的重力异常不仅弱，而且也失去地质体在浅部具有的明显线状异常或梯度带的特征（图 6（a）。当含有地质体与控凹边界断层的模型（图 5（a）埋深很大时，重力异常所刻画的地质体的边界与断层所特有的异常消失，给识别地质体的边界与断层带来困难。但经过对图 6（a）所示的较深埋深地质体重力异常进行重力位错窗口相关处理后（图 6（b），地质边界与断层又被较为清晰地刻画出来，说明重力位错窗口相关法具有从重力异常中增强地质边界与断层的作用。2.4效果对比2.4.1位错窗口相关与方向导数前人常采用重力异常的方向导数方法对地质体的边界、断裂进行识别与划分。方向导数方法能够突出与求导方向相垂直方向上的线状构造。图 7（a）与图 7（b）分别为正常埋深地质体与控凹边界断层模型（图 5（b）重力异常沿x与y方向的导数，只突出了与求导方向相垂直方向的构造，不能如图 5（c）那样把全部存在的构造特征整体刻画出来。DocumentPr oc DocuesClrr os Docu:yanMgtYr ot wuBk?90786654423018图5正常埋深地质体与控凹边界断层模型的重力位错窗口相关效果Fig.5 Effect of gravity dislocation window correlation between normal buried geological body and sag-controlling boundary fault modelD oc umentProcsCl:yD oa umentPMagYwen BlaktP8369554127图6较深埋深地质体与控凹边界断层模型的重力位错窗口相关效果Fig.6 Effect of gravity dislocation window correlation between deep buried geological body and sag-controlling boundary fault model126第 42 卷 第 3 期李希元 等：基于重力异常识别断裂的位错窗口相关法及其在大兴安岭与邻区的应用2.4.2位错窗口相关与小子域滤波小子域滤波是增强重力异常中的线状异常的一种常用的非线性滤波方法，可以利用重力异常划分断裂的同时增强重力线状异常。通过对正常深埋地质体与控凹边界断层模型（图 6（b）重力异常的小子域滤波（图 8（a）与同一异常的位错窗口相关（8（b）相比可知，小子域滤波方法只是将重力异常的密集带变得更加清晰，不能像位错窗口相关那样反映出整个模型各地质体清晰的地质边界。对较深埋深地质体产生的异常（图 6（a）分别进行位错窗口相关处理（图 9（a）及小子域滤波处理（图 9（b）可知，较深埋深地质体在地表产生的重力异常应用位错窗口相关能较为清图7正常埋深地质体与控凹边界断层模型重力异常的方向导数处理效果对比Fig.7 Comparison of processing effect of directional derivative of gravity anomaly between normal buried geological body and sag-controlling boundary fault modelD ocumentProsCl:D ocumeya MgCl:ecumentea Y YY YY YY YY YY YY YlolwolBoBwcko2图8正常埋深地质体与控凹边界断层模型的位错窗口相关与小子域滤波处理效果对比Fig.8 Comparison of processing effect between filtering in small domain and dislocation window correlation of normal buried geological body and sag-controlling boundary fault modelD ocumentProsCl:yaD ocumentMagYw:y B BB BB BB BB BB Blclk4l332c39ment/Galm图9较深埋深地质体与控凹边界断层模型位错窗口相关与小子域滤波处理效果对比Fig.9 Comparison of processing effect between filtering in small domain and dislocation window correlation of deep buried geological body and sag-controlling boundary fault model1272023 年大庆石油地质与开发晰地反映出地质体的轮廓与边界；而小子域滤波处理后，除了能够显示出特别强的边界的特征外，其余的边界则变得十分模糊，展示了位错窗口相关增强与识别被背景掩盖下弱重力线状异常的能力。3应用实例3.1重力异常特征图 10 是大兴安岭与邻区重力异常平面等值线图。从图 10可以看出，位于图中央呈反“S”形的大兴安岭重力梯度带从北向南纵贯全区。梯度带东侧为松辽盆地，西侧为大兴安岭，由于受艾利地壳均衡的作用，松辽盆地的上地幔为隆起区，大兴安岭的上地幔为深凹陷区。正是这样的深部地质结构使得梯度带东侧呈现区域正异常，西侧呈现区域负异常的重力异常特点。此外，在重力梯度带东西两大正、负区域重力异常之上，还有众多走向各异的局部异常，这些局部异常为地壳浅部不同埋深、不同性质的地质体产生。3.2方法应用重力异常是由地表到地下深处密度不均匀体共同作用引起的。为研究大兴安岭与邻区的深部断裂构造，对重力异常进行了小波多尺度分解1819，获得不同深度的重力异常。本文对大兴安岭与邻区的重力异常进行了6阶分解，经对各阶分解结果平均径向功率谱深度估计2021，分解的5阶小波细节重力异常（图11）的功率谱估计平均深度为 32.3 km（图 12），该深度的重力场应是地壳及更深地质体产生的重力异常。在 32.3 km 深度的重力异常体现了大的构造界面起伏及与深大断裂引起的块段效应产生规模宏大、宽缓重力异常的特征。图 11 中的大规模重力异常梯度带，要么是构造界面的陡变带，要么是深大断裂所处的位置。为确定深部断裂或地质构造陡变带的位置，对重力异常5阶小波细节进行了位错窗口相关处理。处理结果（图 13）非常规律地显示了断裂构造的特征。图10大兴安岭与邻区重力异常Fig.10 Gravity anomaly in Greater Hinggan Mountains and its adjacent areas图11大兴安岭与邻区深源重力异常Fig.11 Deep source gravity anomaly in Greater Hinggan Mountains and its adjacent areas128第 42 卷 第 3 期李希元 等：基于重力异常识别断裂的位错窗口相关法及其在大兴安岭与邻区的应用3.3效果分析利用大兴安岭与邻区的重力异常，采用滑动趋势分析获取了反映莫霍面重力异常，应用帕克密度界面三维反演方法2224得到了莫霍面深度（图14），很好地反映了研究区内地幔起伏的情况，即大兴安岭反“S”形重力梯级带的东西两侧表现为幔凹与幔隆。从莫霍面深度与重力位错窗口相关叠合图（图 15）可以看到，反映深部陡边带或是深大断裂的重力位错窗口相关狭窄的线性区与莫霍面的起伏及形态有良好的对应关系，说明这些深部断裂大多切割到上地幔，在一定程度上可以认为，图 13 所示的深大断裂多为超壳深断裂。大兴安岭与邻区的几条重要深大断裂均在重力位错窗口相关图上且均呈现连续分布的狭窄线状异常，也是控制大兴安岭与邻区深部构造的主要大断裂。具体表现为：（1）在松辽盆地内的一条规模宏大的孙吴双辽断裂被清楚地刻画出来，并对莫霍面有一定的控制作用，间接说明该断裂是一条岩石圈断裂（图 15）；（2）具有长期发育历史、被多期断裂构造改造且规模宏大的中朝地台北缘赤峰-开原断裂不仅控制了莫霍面的扭曲与深部构造形态，而且在重力位错窗口相关图上也清楚地表明该断裂曾受多期、不同方向的断裂构造的切割与改造，把整个大断裂肢解得支离破碎，呈现断续的狭窄位错窗口相关线状形态，即使这样也能从位错窗口相关图上看出该深大断裂的延展特征并能指示出被后期断裂改造后呈现的断裂形态特征（图 15）；（3）在深源重力位错窗口相关图（图 13）上大多表现为图14大兴安岭与邻区莫霍面深度Fig.14 Moho depth in Greater Hinggan Mountains and its adjacent areas图13大兴安岭与邻区深源重力位错窗口相关处理效果Fig.13 Processing effect of deep source gravity dislocation window correlation in Greater Hinggan Mountains and its adjacent areas图12大兴安岭与邻区深源重力异常对数功率谱Fig.12 Logarithmic power spectrum of deep source gravity anomaly in Greater Hinggan Mountains and its adjacent areas1292023 年大庆石油地质与开发NE、NEE 向及近NS向的线状构造，NW 向的线状构造在大兴安岭及其邻区，尤其是松辽盆地所属区域少量发育或不发育，这在一定程度上说明，NW向断裂构造不是深部的主要控制断裂，而是形成在地壳深度范围中的壳断裂、区域性断裂。在图 15 中较好地指示了德尔布干断裂、大兴岭莫霍面陡变带、嫩江断裂、孙吴双辽断裂以及赤峰开源断裂的位置，为深化大兴安岭及其邻区深部地质规律的认识提供了新的依据。4结论（1）重力异常位错窗口相关法是一种新型的突出断裂线状异常的非线性滤波方法，在突出与增强重力异常线状特征方面具有显著优势。（2）重力异常位错窗口相关法在识别深源断裂相关重力线状异常方面具有较高的分辨能力。（3）重力异常位错窗口相关法不仅能够识别断裂的线状异常信息，而且对于具有线性分布的异常紧密带也具有极强的分辨能力，能够很好地应用重力异常确定地质体的边界。（4）德尔布干断裂、大兴安岭莫霍面陡变带、嫩江断裂、孙吴-双辽断裂以及赤峰-开源断裂在大兴安岭与邻区深源重力位错窗口相关结果中得到了较好地体现，表明了重力异常的位错窗口相关法具有较好的实用价值。参考文献：1 李成立，陈树民，姜传金.重磁勘探基础理论与资料处理解释方法M.北京：科学出版社，2015.LI Chengli，CHEN Shumin，JIANG Chuanjin.Basic theory of gravity and magnetic exploration and data processing and interpretation methodM.Beijing：Science Press，2015.2 张太富.断裂构造对黔东北铅锌矿的成矿作用J.资源信息与工程，2016，31（2）：25-26.ZHANG Taifu.Mineralization of lead-zinc deposits in northeast Guizhou province by fault structure J.Resources Information and Engineering，2016，31（2）：25-26.3 杨明桂，徐梅桂，胡青华，等.鄂皖赣巨型矿集区的构造复合成矿特征J.地学前缘，2016，23（4）：129-136.YANG Minggui，XU Meigui，HU Qinghua，et al.The structural composite metallogenic characteristics of Hubei-Anhui-Jiangxi giant ore concentration area J.Earth Science Frontiers，2016，23（4）：129-136.4 杨明桂，祝平俊，熊清华，等.新元古代早古生代华南裂谷系 的 格 局 及 其 演 化 J.地 质 学 报，2012，86（9）：1367-1375.YANG Minggui，ZHU Pingjun，XIONG Qinghua，et al.Framework and evolution of the Neoproterozoic-Early Paleozoic south-China rift system J.Acta Geologica Sinica，2012，86（9）：1367-1375.5 李成立，崔瑞华，刘益中.盆地基底岩性的综合地球物理预测方法：以松辽盆地滨北地区基底岩性预测为例 J.地球物理学报，2011，54（2）：491-498.LI Chengli，CUI Ruihua，LIU Yizhong.Comprehensive geophysical prediction method of basement lithology：Example of Binbei area，Songliao Basin J.Chinese Journal of Geophysics，2011，54（2）：491-498.6 崔瑞华，谷社峰，李成立.综合物探在滨北地区基底结构研究中 的 应 用 J.大 庆 石 油 地 质 与 开 发，2011，30（1）：170-174.CUI Ruihua，GU Shefeng，LI Chengli.Application of integrated geophysical technique in basement structure research in Binbei area J.Petroleum Geology&Oilfield Development in Daqing，2011，30（1）：170-174.7 闫百泉，孙杨，杨贵茹，等.松辽盆地深层断陷分布特征及其控制因素J.大庆石油地质与开发，2020，19（1）：1-9.YAN Baiquan，SUN Yang，YANG Guiru，et al.Distribution characteristics and their controlling factors of the deep fault 图15大兴安岭与邻区莫霍面深度与深源重力位错窗口相关叠合Fig.15 Superimposition between Moho depth and deep source gravity dislocation window correlation in Greater Hinggan Mountains and its adjacent areas130第 42 卷 第 3 期李希元 等：基于重力异常识别断裂的位错窗口相关法及其在大兴安岭与邻区的应用depression in Songliao Basin J.Petroleum Geology&Oilfield Development in Daqing，2020，19（1）：1-9.8 姜冠一.徐家围子断陷升平兴城地区深层火山岩岩性、岩相特征及喷发期次J.大庆石油地质与开发，2021，40（4）：9-25.JIANG Guanyi.Lithology，lithofacies characteristics and eruption stages of deep volcanic rocks in Shengping-Xingcheng area of Xujiaweizi rift J.Petroleum Geology&Oilfield Development in Daqing，2021，40（4）：9-25.9 潘磊，徐祖新，李让彬，等.川东南涪陵地区基底断裂特征与油气成藏J.特种油气藏，2020，27（4）：19-25.PAN Lei，XU Zuxin，LI Rangbin，et al.Basement fault characterization and hydrocarbon accumulation in Fuling of southeastern Sichuan J.Special Oil&Gas Reservoirs，2020，27（4）：19-25.10 马海陇，王震，邓光校，等.塔里木盆地和田河东地区断裂特征及其油气地质意义 J.断块油气田，2021，32（3）：329-334.MA Hailong，WANG Zhen，DENG Guangxiao，et al.Fault features in eastern Hetianhe area，Tarim Basin and its petroleum geological significance J.Fault-Block Oil&Gas Field，2021，32（3）：329-334.11 郇恒飞，郭灿文，孙守亮，等.重力数据处理中几种断裂构造识别方法的组合应用 J.地质与资源，2020，29（5）：467-475.HUAN Hengfei，GUO Canwen，SUN Shouliang，et al.Combined application of several methods for fault structure identification in gravty data processingJ.Geology&Resources，2020，29（5）：467-475.12 雷晓东，戚帮申，关伟，等.北京平原区断裂构造重力异常识别研究J.地球物理学报，2021，64（4）：1253-1265.LEI Xiaodong，QI Bangshen，GUAN Wei，et al.Research on the faults identification based on gravity anomaly in Beijing plainJ.Chinese Journal of Geophysics，2021，64（4）：1253-1265.13 舒晴，朱晓颖，周坚鑫，等.矩形棱柱体重力梯度张量异常正演计算公式J 物探与化探，2015，39（6）：1217-1222.SHU Qing，ZHU Xiaoying，ZHOU Jianxin，et al.Forward modeling expressions for right rectangular prism with constant densityJ.Geophysical and Geochemical Exploration，2015，39（6）：1217-1222.14 孙鹏飞，吴燕冈利用重磁水平和垂直二阶导数确定东北地区梯度带J 吉林大学学报（地球科学版），2007，37（增刊1）：27-31.SUN Pengfei，WU Yan gang.Make use of vertically and horizontal differential coefficient of the gravitation and magnetic force to ensure the grads of the northeast China J.Journal of Jilin University（Earth Science Edition），2007，37（