2. 地质过程与矿产资源国家重点实验室, 武汉 430074;
3. 中国地质科学院, 北京 100037
2. Key Laboratory of Geo-detection, Ministry of Education, Wuhan 430074, China;
3. Chinese Academy of Geological Sciences, Beijing 100037, China
岩石圈与软流圈边界(LAB)的物理和化学起因及其控制板块运动的机制是长期困扰固体地球科学家们的问题之一[1-3].LAB最早是流变学上的定义.岩石圈意指地球上力学性质强的外壳, 可以在地质历史上长时间支撑差异应力.岩石圈之下为力学性质上弱的软流圈, 软流圈容许与均衡调整相关的物质流动.而今, 根据测定方法的不同岩石圈已经被定义成岩性岩石圈、地震学岩石圈、电性岩石圈、温度岩石圈、弹性岩石圈等[4].岩石学方法受限于包裹体的获取, 因而只能给出大范围的LAB信息[5].大地电磁学(MT)和地震学方法是测定LAB的主要手段[2-4, 6].一般地, 将MT反演所获得的上地幔第一个相对高导层的顶部视为LAB[6-7].欧洲大陆的对比研究表明[8]:根据地震波速度各向异性和接收函数成像确定的LAB与根据MT探测到的上地幔高导层所推断的LAB在显生宙的中欧地区较为一致, 而在其它地区差别很大, 说明至少存在两种可能的原因:①在不同的大地构造单元, LAB的地质属性不同, 因而在地球物理参数上并不具有同一性.也就是, LAB并不同时具有低速、高导以及速度强各向异性在LAB深度上截止等地球物理性质[9-10]; ②LAB横向变化大, 目前采用的地震学和MT方法不能精确确定其深度.相较于地震学方法, 利用MT确定LAB具有成本低的优势, 但根据单个测点资料确定的LAB精度有限.从MT具有空间平均效应和大量实测资料的结果来看, 它给出的结果一般是稳定的、在平均意义上是可靠的.
华南地区以东南沿海发育大规模火成岩和扬子板块内部发育巨厚元古界地层为主要特征.在华南地区业已开展的MT工作主要与油气勘探有关, 与LAB有关的研究工作则相对较少.李立[11]根据我国早期MT探测成果对中国大陆地壳上地慢电性特征进行了分析.徐常芳[12]探讨了中国大陆地壳上地慢电性结构和地震分布规律之间的可能关联.徐克定[13]通过多条大地电磁剖面, 对华南部分区域的大陆岩石圈层圈结构进行了研究.孙洁等[14]根据MT资料研究了华南地区岩浆活动与深部电性结构的关系.罗志琼等[15]根据375个MT测深点资料确定了中扬子地区的岩石圈厚度, 并探讨了其对该区油气勘探的意义.刘国兴等[16]给出了东南沿海地区岩石圈电性结构, 并讨论了其地球动力学意义.除MT外, 前人在华南地区利用地震学和地热学方法也开展了诸多研究工作[17-20], 但是不同方法定义的岩石圈不同, 各种方法所估计的LAB深度存在一定差异.
在SinoProbe-01-03课题资助下, 我们首次获得了扬子板块和华南造山带及海南岛4°×4°网格共12个测点最长周期为1万秒的MT资料, 并根据一维Occam反演[21]结果采用梯度法确定了LAB的深度, 计算了相应的岩石圈积分电导.为了验证结果的可靠性, 进行了必要的一维数值试验.本文主要报道上述结果, 并讨论所获得的LAB及岩石圈积分电导的地质意义.
2 根据一维MT反演结果确定LAB的模型数值试验前已述及, 一般将上地幔第一个高导层的顶界面视为LAB.对MT数据进行一维反演可以获得各向同性或者各向异性层状电性结构, 但根据反演结果一般难以准确读取上地幔第一个高导层的顶界面.由于LAB下方导电性迅速增强, 因此学者们将反演所获得模型的电阻率梯度最大处当做LAB的位置, 但仍缺乏理论模型的数值试验结果作为依据.进一步地, 当存在或缺乏壳内高导层、电性层存在方位各向异性等情形, 对采用梯度法确定上地幔高导层顶面的精度尚需进一步论证.为了研究华南地区长周期MT资料确定的LAB的可信度, 我们考虑华南地区典型的岩石圈一维地电结构, 进行了大量的数值试验.为简洁计, 此处给出2种模型(见表 1)的数值试验结果:
(1)存在壳内高导层和高导厚层软流圈的电各向同性层状模型, (2)电性层为方位各向异性模型.数值试验中对模型正演所得数据均加入了2%的高斯噪声, 初始模型均为具有1000 Ωm的30层对数域等厚层状模型.
对于电各向同性层状模型, 当同时存在壳内高导层和高导厚层软流圈时, 一维Occam法反演结果可以很好地恢复真实地电结构(见图 1a).此时无论采取极值法(把莫霍面以下第一个相对低阻层的上界面作为LAB)还是梯度法, 估计的LAB深度值均为127.225 km, 与理论模型的125 km十分接近.
对于方位各向异性层状模型, 采用了一维各向异性反演(图 1b, 1c)和各向同性反演(图 1d).各向异性反演采用Pek & Santos (2006)提出的方法[22], 输入数据含全部阻抗元素.各向同性反演采用Occam法, 输入数据为平均视电阻率.利用各向异性程序进行反演与真实模型非常接近, 根据电阻率最大梯度确定的LAB深度为127.2 km.根据笔者的经验, 各向异性反演对原始数据的质量要求很高, 实测数据往往很难满足.用平均视电阻率数据进行各向同性一维Occam反演(图 1d), 对所得模型用梯度法进行估计, LAB的深度在127 km附近, 与理论模型的125 km十分接近.这个数值试验表明:对于电性方位各向异性层状模型, 用平均视电阻率作为输入数据进行各向同性反演, 对反演结果应用梯度法估计LAB的深度具有较高的精度.
3 长周期MT资料的观测与数据处理在华南地区进行了4°×4°网格的MT数据采集, 计12个测点(表 2).每个测点观测时间不少于7天, 估计MT响应函数的频率范围为320 Hz~10000s.投入的采集仪器包括:加拿大凤凰公司生产宽频大地电磁仪MTU-Net、乌克兰LCISR研究所生产的长周期大地电磁仪LEMI-417M以及德国Metronix公司生产的全频带大地电磁仪GMS07e.采集方式有两种:(1) MTU-Net和LEMI-417M组合, (2) GMS07e系统(三根感应磁棒和FGS03三分量磁通门磁力仪组合).
采用人工时间序列挑样或功率谱筛选、Robust估计[23]、远参考处理[24]、阻抗张量分解[25-26]等技术进行数据处理, 最终获得视电阻率与相位曲线.利用WAL方法[27-28]判别了10s、100s、1000s及10000s周期地电结构的维数性(图 2), 计算了对应的平均Parkinson感应箭头[29](图 3).由于华南地区干扰因素比较多、地下结构复杂, 相位曲线也发生了畸变而且难以消除, 因此本文利用平均视电阻率(荆门测点例外, 由于ρyx质量差, 故只采用ρxy视电阻率曲线)进行一维Occam反演, 获得的一维地电模型及拟合情况分别见图 4-图 8.采用梯度法估计了每个测点的LAB深度(图 9a), 计算了岩石圈积分电导(LAB以上总纵向电导)(图 9b)和平均电阻率(表 2).
华南地区在10~104s周期范围内主要为三维电性结构(图 2), 不存在简单的一维和二维电性结构.在不同的深度范围内高导体的分布差异性明显(图 3).华南地区赣州、邵阳、施秉、云浮及揭阳测点的Parkinson感应箭头的幅值大, 其余测点均较小(图 3b).湖州、赣州、揭阳测点在10s周期的Parkinson感应箭头指向(SW 15°-25°)基本一致, 云浮测点Parkinson感应箭头指向发生约90°偏转(WNW方向), 它们均指示上地壳横向上由高阻体向高导体的转换.邵阳、施秉及云浮测点在1000s周期的Parkinson感应箭头指向基本一致(约SE 15°-25°), 指示岩石圈深部横向上由高阻体向高导体的转换.由于点距近400 km, 对Parkinson感应箭头的上述分析只对测点附近的电阻率横向变化有意义.
根据4°×4°网格长周期MT一维反演结果, 华南地区壳内和上地幔一般存在四个高导层:第一个高导层的深度在10~25 km范围内, 电阻率约10~60 Ωm, 为壳内高导层, 主要与发育壳内剪切带产生蛇纹岩化及含盐流体相关[30].第二个高导层在25~50 km之间, 电阻率约10 Ωm, 对应于电性Moho[31].除湖南邵阳和贵州施秉及广东揭阳测点外, 第三个高导层的深度在60~145 km范围, 电阻率低, 小于或等于10 Ωm, 对应于电性LAB; 第四个高导层的深度在200~350 km之间, 电阻率为10~30 Ωm, 为幔内高导层, 推测为幔内熔融体.沿北纬27°线LAB深度较大, 在120~230 km之间; 沿31°线测点的LAB深度次之, 在60~120 km之间.研究结果并不支持华南地区电性岩石圈东薄西厚的观点, 而是显示出岩石圈南北两侧上抬、中部下凹、东部受不均匀改造的趋势, 反映华南地区岩石圈稳定性较好, 晚中生代以来的构造伸展作用[32-33]对岩石圈的改造程度有限, 可能主要以不同形式的软流圈底辟[34]为主.
根据一维地电结构可以将华南地区岩石圈划分为五种类型:克拉通型(图 4)、构造边界型(图 5)、岩石圈中等改造型(图 6)、岩石圈强烈改造型(图 7)及造山带型(图 8).
华南克拉通型地电结构(图 4)表现为电性岩石圈厚度大、电性Moho清晰、岩石圈平均电阻率值高(表 2).湖南邵阳和贵州施秉测点实测视电阻率曲线类型与加拿大Slave克拉通基本相同(对比文献[31] 图 3), 估计的LAB分别为238.5 km和183.2 km, 岩石圈平均电阻率分别达到3223 Ωm和1489 Ωm.地热学研究[35]表明湖南邵阳和贵州施秉测点均位于平均热流值等于或低于50 mW·m-2区域, 在华南地区属于"冷"岩石圈.
构造边界型岩石圈电性结构以四川达州和彭州及湖北荆门测点为代表(图 5), 以无电性Moho为主要特征, 其它特征还包括岩石圈导电性好(平均电阻率小于60 Ωm)、中上地壳高导层发育、LAB以下电阻率低( < 20 Ωm).四川达州和彭州测点的电性结构指示地壳由弱的上地壳和强的下地壳组成, 上地壳与下地壳在力学性质上解耦, 佐证了四川地块发育双层变质基底且存在上部变质层滑脱[36].地热学[35]研究指示四川达州和彭州测点均位于热流值50~60 mW·m-2区域, 在正常热状态下, "湿"岩石圈是解释这一类电性结构的选项.因此, 达州和彭州位于四川地块与周边构造带的转换边界, 由发育岩石圈尺度的剪切带而形成高的导电性.湖北荆门测点位于太行山至武陵山重力梯度带上, 岩石圈厚度明显比其两侧小, 推测与所在构造带软流圈的上升有关.
岩石圈中等改造型电性结构以浙江湖州和广东云浮测点为代表(图 6), 主要特征为岩石圈导电性好(平均电阻率小于60 Ω.m)、上地壳电阻率远高于下地壳、LAB以下电阻率中等(20~50 Ωm).这种类型岩石圈的地壳由强的上地壳和弱的下地壳组成, 推测上地壳与下地壳在力学性质上解耦.地热学[35]研究指示浙江湖州和广东云浮测点均位于热流值60~70 mW·m-2区域, 暗示岩石圈遭受热侵蚀, 或有局部软流圈底辟.浙江湖州测点位于前人推测的大别-舟山岩石圈断裂带[34]和新安江壳幔韧性剪切带[19]端点处附近, 而广东云浮测点位于前人推测的赣湘桂弱岩石圈减薄带和华夏强岩石圈减薄带的转换带上[19], 故它们代表一类岩石圈中等改造型电性结构.
岩石圈强烈改造型电性结构以江西赣州、广东揭阳及福建霞浦测点为代表(图 7), 主要特征为岩石圈导电性中等(平均电阻率处于400~1200 Ωm之间)、中下地壳高导层不发育、电性Moho不明显而LAB清晰、岩石圈厚度偏大、LAB以下电阻率高(均值大于100 Ωm).上述3个测点的实测视电阻率曲线的尾支明显上升, 指示这种类型岩石圈对应的软流圈较薄.地热学[35]研究指示上述测点均位于平均热流值大于70 mW·m-2区域, 暗示岩石圈遭受强烈热侵蚀.值得注意的是上述测点位于前人所称的华夏岩石圈强减薄区域[19], 但长周期MT资料揭示其电性岩石圈厚度均超过130 km, 与中扬子正常岩石圈厚度相当.刘国兴等[16]完成的华南沿海地区MT二维反演结果也揭示江西赣州和会昌之间及潮州附近岩石圈厚度在100 km以上.如果用大洋板块俯冲时刮削下来的洋壳沉积物叠置在岩石圈下部导致岩石圈增厚机制[16]来解释福建霞浦和广东揭阳测点的电性结构, 则无法解释这两处偏高的岩石圈平均电阻率.另外, 江西赣州测点的电性结构不符合造山带"厚壳薄圈"特征(见后续分析), 因而不能用陆内造山机制[16]来解释.一种可能的机制是福建霞浦、江西赣州及广东揭阳测点在晚中生代末期以来构造伸展背景下发生了岩石圈底部的拆沉和垫托作用, 并伴随有大规模岩浆侵入事件, 造成改造后的岩石圈厚度增大而软流圈变薄, 同时由于壳内大量侵入岩和超镁铁质岩石的板底垫托使得岩石圈的平均电阻率偏高.东南沿海地区岩石圈遭受上述改造作用在时间和空间上极不均一, 某些区域的岩石圈拆沉作用可能还在进行之中, 如梅州-平和-长泰一带显示明显的软流圈物质上升[16], 而另一些区域已完成岩石圈的"新生", 因此简单的"减薄"机制难以解释东南沿海地区现今的岩石圈电性结构.
造山带型岩石圈电性结构以湖北英山测点为代表(图 8), 主要特征是厚壳薄幔且地壳电阻率高、电性Moho清晰.
5 结论(1) 利用大地电磁法探测LAB, 通过对不同的地电模型进行反演模拟发现:地表低阻层、地表方位各向异性对估计LAB没有影响; 当存在低阻软流圈时, 能较准确地判断LAB位置; 利用梯度法估计LAB位置具有较高的精度.
(2) 华南地区4°×4°长周期MT测点所获得的电性岩石圈并不是简单的东薄西厚, 而是显示出南北两侧上抬、中部下凹、东部受不均匀改造的趋势, 反映华南地区岩石圈稳定性较好, 晚中生代以来的构造伸展作用对岩石圈的改造程度有限, 而且在时间和空间上不均一.
(3) 本文将华南岩石圈电性结构划分成克拉通型、构造边界型、岩石圈中等改造型、岩石圈强烈改造型及造山带型, 分别描述了它们的主要特征, 对今后的长周期MT资料解释提供了参照物.华南地区经历了多期大规模构造运动, 岩石圈结构和构造复杂.加密测点和延长观测周期对于利用长周期MT资料研究华南地区岩石圈结构和构造、预测大型成矿带及油气田远景区均具有重要意义.
致谢感谢中国地质大学(武汉)地球物理与空间信息学院蒋礼、张秉政、周文豪、张西君、王重阳等在采集华南地区4°×4°网度高质量长周期MT数据中所付出的艰辛劳动.感谢Sinoprobe-01项目对本文的资助.
[1] | Fischer K M, Ford H A, Abt D L, et al. The lithosphere-asthenosphere boundary. Annu. Rev. Earth Planet. Sci. , 2010, 38: 551-575. DOI:10.1146/annurev-earth-040809-152438 |
[2] | Williams M L, Fischer K M, Freymueller J T, et al. Unlocking the Secrets of the North American Continent: An EarthScope Science Plan for 2010-2020, February, 2010, 78 pp. http://en.cnki.com.cn/Article_en/CJFDTOTAL-GJZT201011001.htm |
[3] | Lay T ed. Seismological Grand Challenges in Understanding Earth's Dynamic Systems. Report to the National Science Foundation, IRIS Consortium , 2009, 76. |
[4] | Eaton D W, Darbyshire F, Evans R L, et al. The elusive lithosphere-asthenosphere boundary (LAB) beneath cratons. Lithos , 2009, 109: 1-22. DOI:10.1016/j.lithos.2008.05.009 |
[5] | O'Reilly S Y, Griffin W L, et al. The continental lithosphere-asthenosphere boundary: Can we sample it?. Lithos , 2010, 120: 1-13. |
[6] | 魏文博, 金胜, 叶高峰, 等. 大陆岩石圈导电性的研究方法. 地学前缘 , 2003, 10(1): 15–22. Wei W B, Jin S, Ye G F, et al. Methods to study electrical conductivity of continental lithosphere. Earth Science Frontiers (in Chinese) , 2003, 10(1): 15-22. |
[7] | 袁学诚. 中国地球物理图集. 北京: 地质出版社, 1996 . Yuan X C. Atlas of geophysics in China (in Chinese). Beijing: Geological Publishing House, 1996 . |
[8] | Jones A G, Plomerova J, Korja T, et al. Europe from the bottom up: A statistical examination of the central and northern European lithosphere-asthenosphere boundary from comparing seismological and electromagnetic observations. Lithos , 2010, 120: 14-29. DOI:10.1016/j.lithos.2010.07.013 |
[9] | Ichiki M, Uyeshima M, Utada H, et al. Upper mantle conductivity structure of the back-arc region beneath northeastern China. Geophysical Research Letters , 2001, 28(19): 3773-3776. DOI:10.1029/2001GL012983 |
[10] | 赵国泽, 詹艳, 王立凤, 等. 鄂尔多斯断块地壳电性结构. 地震地质 , 2010, 32(3): 345–359. Zhao G Z, Zhan Y, Wang L F, et al. Electricstructure of the crust beneath the Ordos fault block. Seismology and Geology (in Chinese) , 2010, 32(3): 345-359. |
[11] | 李立. 中国大陆地壳上地慢电性特征. 地球物理学报 , 1996, 39(supp): 130–140. Li L. The geoelectrial characteristic of crust and upper mantle in the continental region of China. Chinese J. Geophys. (Acta Geophysica Sinica) (in Chinese) , 1996, 39(supp): 130-140. |
[12] | 徐常芳. 中国大陆地壳上地慢电性结构及地震分布规律(一). 地震学报 , 1996, 18(2): 254–261. Xu C F. The earthquake distribution and the resistivity structure in the Chinese mainland. Acta Seismologica Sinica (in Chinese) , 1996, 18(2): 254-261. |
[13] | 徐克定. 中国南方大陆岩石圈结构分析. 海相油气地质 , 2000, 5(3/4): 31–46. Xu K D. The structure analysis of South China mainland lithosphere. Marine Origin Petroleum Geology (in Chinese) , 2000, 5(3/4): 31-46. |
[14] | 孙洁, 晋光文, 江钊, 等. 用大地电磁测深法研究华南地区岩浆活动及深部结构. 地震地质 , 2001, 23(2): 328–329. Sun J, Jin G W, Jiang Z, et al. Magmatic activity and the deep structure of South China inferred from magnetotelluric sounding data. Seismology and Geology (in Chinese) , 2001, 23(2): 328-329. |
[15] | 罗志琼, 王家映, 邓忠凡, 等. 中扬子地区大地电磁与油气构造研究. 石油物探 , 1996, 35(3): 71–81. Luo Z Q, Wang J Y, Deng Z F, et al. The study of MT and hydrocarbon tectonic in central Yangzi area. Geophysical Prospecting for Petroleum (in Chinese) , 1996, 35(3): 71-81. |
[16] | 刘国兴, 韩凯, 韩江涛. 华南东南沿海地区岩石圈电性结构. 吉林大学学报(地球科学版) , 2012, 42(2): 536–544. Liu G X, Han K, Han J T. Lithosphere electrical structure in southeast coastal region, South China. Journal of Jilin University (Earth Science Edition) (in Chinese) , 2012, 42(2): 536-544. |
[17] | 安美建, 石耀霖. 中国大陆岩石圈厚度分布研究. 地学前缘 , 2006, 13(3): 23–30. An M J, Shi Y L. Review on lithospheric thickness research of the Chinese continent. Earth Science Frontiers (in Chinese) , 2006, 13(3): 23-30. |
[18] | 李廷栋. 中国岩石圈的基本特征. 地学前缘 , 2010, 17(3): 1–13. Li Y T. Principal characteristics of the lithosphere of China. Earth Science Frontiers (in Chinese) , 2010, 17(3): 1-13. |
[19] | 蔡学林, 朱介寿, 曹家敏, 等. 华南地区岩石圈三维结构类型与演化动力学. 大地构造与成矿学 , 2003, 27(4): 301–312. Cai X L, Zhu J S, Cao J M, et al. Three-dimensional tectonic types and evolutional dynamics of lithosphere of South China region. Geotectonica et Metallogenia (in Chinese) , 2003, 27(4): 301-312. |
[20] | 朱介寿, 蔡学林, 曹家敏, 等. 中国及相邻区域岩石圈结构及动力学意义. 中国地质 , 2006, 33(4): 793–803. Zhu J S, Cai X L, Cao J M, et al. Lithospheric structure and geodynamics in China and its adjacent areas. Geology in China (in Chinese) , 2006, 33(4): 793-803. |
[21] | Constable S C, Parker R L, Constable C G. Occam's inversion: a practical algorithm for generating smooth models from electromagnetic sounding data. Geophysics , 1987, 52: 289-300. DOI:10.1190/1.1442303 |
[22] | Pek J, Santos F A M. Magnetotelluric inversion for anisotropic conductivities in layered media. Physics of the Earth and Planetary Interiors , 2006, 158: 139-158. DOI:10.1016/j.pepi.2006.03.023 |
[23] | Egbert G D, Booker J R. Robust estimation of geomagnetic transfer functions. Geophysical Journal of the Royal Astronomical Society , 1986, 87(1): 173-194. DOI:10.1111/gji.1986.87.issue-1 |
[24] | Gamble T D, Goubau W M, Clarke J. Magnetotellurics with a remote magnetic reference. Geophysics , 1979, 44(1): 53-68. DOI:10.1190/1.1440923 |
[25] | Groom R W, Bailey R C. Decomposition of magnetotelluric impedance tensors in the presence of local three-dimensional galvanic distortion. Journal of Geophysical Research , 1989, 94(B2): 1913-1925. DOI:10.1029/JB094iB02p01913 |
[26] | McNeice G W, Jones A G. Multisite, multifrequency tensor decomposition of magnetotelluric data. Geophysics , 2001, 66(1): 158-173. DOI:10.1190/1.1444891 |
[27] | Weaver J T, Agarwal A K, Lilley F E M. Characterisation of the magneto-telluric tensor in terms of its invariants. Geophysical Journal International , 2000, 141: 321-336. DOI:10.1046/j.1365-246x.2000.00089.x |
[28] | Martí A, Queralt P, Ledo J. WALDIM: A code for the dimensionality analysis of magnetotelluric data using the rotational invariants of the magnetotelluric tensor. Computers & Geosciences , 2009, 35: 2295-2303. |
[29] | Parkinson W D. The influence of continents and oceans on geomagnetic variations. Geophys. J. R. astr. Soc. , 1962, 6: 441-449. DOI:10.1111/gji.1962.6.issue-4 |
[30] | 徐义贤. 中下地壳高导层成因研究综述. 地质科技情报 , 1995, 14(3): 15–22. Xu Y X. Review of the mechanism of the high conduction layers in middle to lower continental crust. Geological Science and Technology Information (in Chinese) , 1995, 14(3): 15-22. |
[31] | Jones A G, Ferguson I J. The electric Moho. Nature , 2001, 409: 331-333. DOI:10.1038/35053053 |
[32] | 毛景文, 谢桂青, 李晓峰, 等. 华南地区中生代大规模成矿作用与岩石圈多阶段伸展. 地学前缘 , 2004, 11(1): 45–55. Mao J W, Xie G Q, Li X F, et al. Mesozoic large scale mineralization and multiple lithospheric extension inSouth China. Earth Science Frontiers (in Chinese) , 2004, 11(1): 45-55. |
[33] | 舒良树. 华南构造演化的基本特征. 地质通报 , 2012, 31(7): 1035–1053. Shu L S. An analysis of principal features of tectonic evolution in South China Block. Geological Bulletin of China (in Chinese) , 2012, 31(7): 1035-1053. |
[34] | 陈沪生, 张永鸿, 徐师文, 等. 下扬子及邻区岩石圈结构构造特征与油气资源评价. 北京: 地质出版社, 1999 : 1 -287. Chen H S, Zhang Y H, Xu S W, et al. The lithospheric textural and structure features as well as oil and gas evaluation in the lower Yangtze area and its adjacent region, China (in Chinese). Beijing: Geological Publishing House, 1999 : 1 -287. |
[35] | 汪洋. 华北和华南岩石圈热状态、流变学特征与地壳成分. 北京: 地质出版社, 2006 : 1 -91. Wang Y. Lithospheric thermal state, rheology and crustal composition of north and south China (in Chinese). Beijing: Geological Publishing House, 2006 : 1 -91. |
[36] | 朱英. 中国及邻区大地构造和深部构造纲要:全国1:100万航磁异常图的初步解释. 北京: 地质出版社, 2004 : 1 -90. Zhu Y. A framework: Geotectonic and deep crustal architectural features of China and some peripheral areas: Interpretation of aeromagnetic map series 1:1000000 (in Chinese). Beijing: Geological Publishing House, 2004 : 1 -90. |