胶莱盆地古地热场：基于数值模拟的深部热历史解析

胶莱盆地古地热场：基于数值模拟的深部热历史解析一、引言1.1研究背景与意义胶莱盆地作为中国东部重要的沉积盆地，其复杂的地质演化历程蕴含着丰富的地球科学信息，对研究区域地质构造和资源分布具有不可替代的作用。胶莱盆地位于山东半岛中部，北临渤海湾，西靠泰山山脉，东部和南部与山东半岛边缘山脉接壤。它经历了多期构造运动和沉积作用，记录了从远古时期到现今的漫长地质历史变迁，这些地质过程对盆地的形成、发展和改造产生了深远影响。其独特的地理位置和复杂的地质背景，使其成为揭示中国东部板块构造演化、沉积环境变迁以及深部地质作用的关键区域。在地质演化研究领域，胶莱盆地提供了一个天然的实验室，帮助我们深入了解地球内部动力学过程和构造变形机制。通过研究盆地的地层序列、沉积相和构造特征，可以重建古地理环境和古气候条件，揭示地球历史上重大地质事件的发生过程和机制，为全球地质演化研究提供重要的区域实例。例如，胶莱盆地在中生代经历了强烈的构造活动和岩浆侵入，这些事件对盆地的地层沉积和构造变形产生了深刻影响，研究这些过程有助于我们理解板块运动和深部地质作用的相互关系。从能源勘探角度来看，胶莱盆地具有巨大的资源潜力，特别是在油气和页岩气等能源资源方面。近年来，随着全球能源需求的不断增长，非常规油气资源的勘探开发受到了广泛关注。胶莱盆地白垩系地层中发育的富有机质页岩，具备形成大型页岩气田的地质条件，有望成为我国未来能源供应的重要接替区。此外，盆地内的沉积岩系中还可能蕴藏着丰富的石油、天然气和煤炭等常规能源资源，对这些资源的勘探开发将有助于缓解我国能源短缺的压力，保障国家能源安全。因此，深入研究胶莱盆地的地质构造和沉积演化，对于准确评估其能源资源潜力、指导能源勘探开发具有重要的现实意义。古地热场作为控制盆地内物质能量交换和地质作用的关键因素，对理解盆地的演化历史和资源形成过程起着至关重要的作用。地热能是地球内部的一种重要能源形式，它通过热传导、对流和辐射等方式在地球内部和地表之间传递，对地球的物理和化学过程产生深远影响。在盆地演化过程中，古地热场的变化控制了沉积物的压实、成岩作用、有机质的热演化以及油气的生成、运移和聚集等关键过程。例如，较高的古地温可以加速有机质的热降解，促进油气的生成；而古地热场的不均匀分布则可能导致油气在盆地内的差异聚集。因此，重建胶莱盆地的古地热场，能够为我们深入理解盆地的构造演化、沉积作用和能源资源形成机制提供关键线索。数值模拟技术作为一种强大的研究手段，为我们定量研究古地热场提供了可能。通过建立数学模型，模拟古地热场的演化过程，可以克服传统地质研究方法在时间和空间尺度上的局限性，更加准确地揭示古地热场的分布特征和演化规律。数值模拟能够综合考虑多种地质因素的相互作用，如岩石的热物理性质、构造运动、沉积速率和热流边界条件等，从而对古地热场进行更加全面和深入的研究。此外，数值模拟还可以预测不同地质条件下古地热场的变化趋势，为地质勘探和资源开发提供科学依据。综上所述，开展胶莱盆地古地热场数值模拟研究具有重要的科学意义和现实价值。通过这项研究，我们可以深化对胶莱盆地地质演化历史的认识，揭示古地热场与构造运动、沉积作用之间的内在联系，为区域地质研究提供新的视角和方法。同时，准确重建古地热场有助于我们更准确地评估胶莱盆地的能源资源潜力，指导油气和页岩气等能源资源的勘探开发，为我国能源战略的实施提供有力支持。1.2国内外研究现状在胶莱盆地地质研究领域，国内外学者已取得了一定成果。国外方面，部分学者通过地质调查和地球物理勘探技术，对胶莱盆地的地层结构和构造特征进行了初步分析。他们利用地震反射剖面，识别出盆地内不同地层的界面和厚度变化，初步揭示了盆地的沉积序列和构造格局，如[国外学者姓名1]等在研究中发现，胶莱盆地存在多个沉积旋回，反映了其复杂的地质演化历史。在盆地构造演化方面，[国外学者姓名2]通过对盆地内构造变形的研究，提出胶莱盆地经历了多期构造运动，包括伸展、挤压等不同构造应力场的作用，这些研究为理解胶莱盆地的地质演化提供了重要的理论基础。国内对胶莱盆地的研究更为深入和全面。在沉积特征研究中，学者们通过对盆地内大量钻井资料和露头的分析，详细阐述了不同地层的沉积相类型和沉积环境变迁。例如，[国内学者姓名1]研究表明，胶莱盆地白垩系地层经历了从湖泊相到河流相再到三角洲相的沉积演化过程，不同沉积相的分布受构造运动和古气候的双重控制。在构造演化方面，众多学者利用多种研究方法，深入探讨了胶莱盆地的形成机制和演化历史。[国内学者姓名2]通过构造应力场模拟和平衡剖面分析，认为胶莱盆地在中生代主要受太平洋板块俯冲的影响，经历了强烈的伸展和断陷作用，形成了现今的构造格局。此外，国内学者还对胶莱盆地的油气地质特征进行了大量研究，为盆地的能源勘探提供了重要依据。然而，现有研究仍存在一定不足。在古地热场研究方面，虽然国内外学者已认识到古地热场对盆地演化和能源资源形成的重要性，但相关研究还相对较少。目前对胶莱盆地古地热场的认识主要基于有限的地热梯度测量和岩石热物性参数分析，缺乏系统的、定量的研究。在数值模拟方面，虽然已有一些学者尝试对盆地的构造演化和沉积过程进行数值模拟，但针对古地热场的数值模拟研究还十分薄弱。已有的数值模拟研究往往未能充分考虑多种地质因素对古地热场的综合影响，如构造运动、岩浆活动、沉积速率等，导致模拟结果的准确性和可靠性有待提高。此外，在古地热场与构造演化、沉积作用之间的耦合关系研究方面，也存在明显的不足，尚未建立起完善的理论模型和研究体系。综上所述，尽管目前对胶莱盆地的地质研究已取得了一定进展，但在古地热场研究领域仍存在诸多空白和不足。本研究旨在通过系统的数值模拟方法，综合考虑多种地质因素，深入研究胶莱盆地的古地热场演化特征，填补该领域在定量研究方面的空白，为全面理解胶莱盆地的地质演化历史和能源资源形成机制提供新的视角和方法。1.3研究内容与方法1.3.1研究内容本研究围绕胶莱盆地古地热场展开多方面研究。首先，对胶莱盆地的地质特征进行全面分析，通过收集整理大量地质资料，包括地层、构造、岩浆活动等信息，明确盆地的地质演化历程。深入研究地层的沉积序列和岩性特征，确定不同地层的沉积环境和时代，为后续古地热场模拟提供基础地质框架。详细分析构造运动对盆地的影响，识别主要断裂构造和褶皱形态，了解构造应力场的演化过程，因为构造活动与古地热场的变化密切相关，断裂和褶皱的形成会改变岩石的热传导路径和热边界条件。同时，关注岩浆活动的时空分布，岩浆侵入会带来大量热能，对古地热场产生显著影响，确定岩浆活动的时期、规模和侵入位置，有助于准确模拟古地热场的演化。在模拟参数选取方面，测定岩石的热物理性质参数，如热导率、比热容和热扩散率等。这些参数是古地热场数值模拟的关键输入，它们决定了岩石中热量的传导和储存能力。不同岩石类型的热物理性质存在差异，例如砂岩和泥岩的热导率不同，会导致热量在不同地层中的传播速度和分布情况不同。因此，通过实验分析获取准确的岩石热物理性质参数至关重要。此外，还需确定边界条件和初始条件，边界条件包括盆地边缘的热流输入和输出，以及与周边地质体的热交换情况；初始条件则是指模拟开始时盆地内的温度分布。这些条件的准确设定对于模拟结果的可靠性至关重要，需要综合考虑地质历史时期的各种因素来确定。构建胶莱盆地古地热场数值模型是研究的核心内容之一。基于有限元法或有限差分法等数值计算方法，建立二维或三维的古地热场模型。在模型中，将盆地划分为多个网格单元，每个单元赋予相应的岩石热物理性质参数和边界条件。通过求解热传导方程，模拟不同地质时期古地热场的分布特征和演化规律。考虑多种地质因素对古地热场的影响，如构造运动导致的岩石变形和地层抬升沉降，会改变热传导路径和地热梯度；岩浆活动带来的热扰动会使局部地区温度升高；沉积作用导致地层厚度增加，会影响热量的积累和散失。通过在模型中合理设置这些因素，能够更真实地模拟古地热场的演化过程。对模拟结果进行分析和验证也是研究的重要环节。将模拟结果与实际地质数据进行对比，如地层温度测量数据、有机质成熟度指标等，评估模拟结果的准确性和可靠性。通过对比分析，发现模拟结果与实际数据之间的差异，分析差异产生的原因，对模型进行优化和改进。进一步探讨古地热场对盆地内地质作用的影响，如对有机质热演化的影响，通过模拟不同古地热场条件下有机质的成熟度变化，确定油气生成的有利区域；对岩石成岩作用的影响，研究古地热场如何影响岩石的压实、胶结和矿物转化等过程，从而揭示古地热场在盆地地质演化中的重要作用。1.3.2研究方法本研究采用多种研究方法相结合的方式。地质调查是基础研究方法之一，通过野外实地考察，详细观察盆地内的地层露头、构造形态和岩浆岩出露情况。绘制地质图件，记录地层的厚度、岩性变化、接触关系以及构造的走向、倾角和规模等信息。采集岩石样品，用于后续的实验分析。在野外调查过程中，注意观察不同地层的沉积构造，如交错层理、波痕等，以推断沉积环境和古水流方向；仔细研究构造变形特征，如断裂的性质、褶皱的紧闭程度等，为分析构造演化提供依据。实验分析是获取关键数据的重要手段。对采集的岩石样品进行实验室测试，测定岩石的热物理性质参数。采用稳态法或瞬态法测量岩石的热导率，通过量热实验测定比热容，利用热扩散率测试仪测定热扩散率。分析岩石的矿物组成和化学成分，因为矿物组成和化学成分会影响岩石的热物理性质，例如石英含量高的岩石热导率相对较高。通过显微镜观察岩石的微观结构，了解岩石的孔隙度和渗透率等储层特征，这些特征也与热传导和热量储存有关。数值模拟是本研究的核心方法。运用专业的数值模拟软件，如COMSOLMultiphysics、FEFLOW等，建立胶莱盆地古地热场数值模型。根据地质调查和实验分析获取的数据，合理设置模型参数和边界条件。通过模拟不同地质时期古地热场的演化，预测古地温的分布和变化趋势。在数值模拟过程中，对模型进行敏感性分析，研究不同参数对模拟结果的影响程度，确定关键参数，提高模拟结果的准确性。利用模拟结果进行地质解释，探讨古地热场与地质作用之间的内在联系。二、胶莱盆地地质特征剖析2.1盆地构造特征2.1.1区域构造背景胶莱盆地位于中国东部，处于华北板块与扬子板块的碰撞带附近，其形成与演化深受区域构造运动的影响。盆地夹持于胶北隆起和胶南造山带之间，这种特殊的大地构造位置，使其成为研究板块相互作用和区域构造演化的关键区域。盆地的西部边界为沂沭断裂带，这是中国东部规模宏大的走滑断裂，也是鲁东与鲁西地区的重要地质分界线，自中生代以来，对山东东部的构造格局产生了深远影响。沂沭断裂带具有复杂的地质构造，它不仅是中新生代岩浆活动、成矿作用和地震活动的重要场所，还控制了胶莱盆地西部的边界形态和沉积演化。众多学者对沂沭断裂带进行了深入研究，[学者姓名1]通过地质调查和地球物理勘探发现，沂沭断裂带在中生代经历了多期构造活动，包括强烈的走滑运动和伸展作用，这些活动导致了断裂带两侧岩石的变形和位移，也为胶莱盆地的形成提供了构造背景。[学者姓名2]的研究进一步表明，沂沭断裂带的活动与太平洋板块向欧亚板块的俯冲密切相关，俯冲作用引发的深部构造应力变化，促使沂沭断裂带发生强烈的构造变形，从而控制了胶莱盆地的形成和演化。盆地的南部边界为五一荣断裂，它是著名的苏鲁大别超高压变质带东端——胶南造山带的北部边界。胶南造山带是华北板块南部活动大陆边缘与扬子板块北部被动大陆边缘碰撞、拼贴的产物，在中国地质发展史上具有极其重要的地位。五一荣断裂的活动对胶莱盆地的南部边界和沉积演化同样具有重要控制作用。[学者姓名3]通过对五一荣断裂带的岩石变形和变质作用研究发现，该断裂在碰撞造山过程中经历了强烈的挤压和逆冲作用，导致断裂带两侧岩石发生复杂的变形和变质，形成了一系列高压、超高压变质矿物和构造形迹。这些构造活动不仅改变了盆地南部的岩石性质和构造格局，还影响了盆地内的沉积环境和地层分布。在中生代，太平洋板块向欧亚板块的强烈俯冲，是影响胶莱盆地形成和演化的关键动力学因素。俯冲作用导致了区域构造应力场的强烈变化，使胶莱盆地所在区域经历了伸展、挤压等复杂的构造变形过程。在伸展作用下，地壳发生拉伸减薄，形成了一系列断陷盆地，胶莱盆地便是在这一背景下逐渐形成的。而在挤压作用阶段，盆地内的地层受到强烈挤压，发生褶皱和断裂，进一步塑造了盆地的构造格局。[学者姓名4]通过构造模拟和地质分析指出，太平洋板块俯冲引发的地幔对流和深部构造应力变化，在胶莱盆地及其周边地区产生了复杂的构造应力场，这种应力场的时空变化控制了盆地内断裂的活动、地层的沉积和变形，对盆地的演化产生了深远影响。综上所述，胶莱盆地的区域构造背景复杂，其形成和演化受到沂沭断裂带、五一荣断裂以及太平洋板块俯冲等多种因素的综合控制。这些因素的相互作用，塑造了胶莱盆地独特的构造格局和地质演化历史，为深入研究盆地的古地热场提供了重要的地质基础。2.1.2盆地内部构造格局胶莱盆地内部构造格局复杂，发育多个次级凹陷和凸起，这些次级构造单元对盆地内的沉积和地热场分布产生了重要影响。盆地主要由近南北向排列的莱阳凹陷、高密凹陷以及夹于其间的大野头凸起组成。莱阳凹陷位于盆地北部，是盆地内重要的沉积中心之一，其沉积地层厚度较大，岩性主要为白垩系的陆相碎屑岩，包括砂岩、泥岩和砾岩等。这些沉积物记录了盆地在不同演化阶段的沉积环境和构造背景信息。高密凹陷则处于盆地南部，同样堆积了较厚的白垩系地层，其沉积特征与莱阳凹陷既有相似之处，也存在一定差异，反映了盆地内部沉积环境的复杂性和多样性。大野头凸起作为盆地内的相对隆起区，其地层遭受了一定程度的剥蚀，出露的地层相对较老，主要为前寒武纪变质岩系，这些岩石的出露揭示了盆地深部的地质构造信息。盆地内部褶皱和断层构造发育，进一步复杂化了盆地的构造格局。褶皱构造主要表现为一系列紧闭褶皱和开阔褶皱，其轴向多为近东西向或北东向，与区域构造应力场方向密切相关。[学者姓名5]通过野外地质调查和构造解析发现，这些褶皱构造是在多期构造运动的作用下逐渐形成的，早期的褶皱受到后期构造运动的叠加和改造，导致褶皱形态复杂多样。断层则以正断层和逆断层为主，其中正断层多为控制盆地边界和次级构造单元的边界断层，如控制莱阳凹陷和高密凹陷边界的断层，它们在盆地演化过程中起到了重要的控盆作用，控制了沉积地层的厚度和分布范围。逆断层则多发育于盆地内部，它们的活动导致地层的错断和变形，对沉积地层的连续性和完整性产生了破坏作用。例如，[学者姓名6]通过地震剖面解释和构造分析发现，一些逆断层的活动使得白垩系地层发生强烈变形和错动，形成了复杂的构造样式，这些构造样式对盆地内的油气运移和聚集产生了重要影响。盆地内部构造格局对沉积和地热场有着显著影响。在沉积方面，次级凹陷为沉积物的堆积提供了场所，凹陷的形态和规模控制了沉积地层的厚度和分布范围。例如，莱阳凹陷和高密凹陷由于其较大的沉降幅度和相对稳定的沉积环境，堆积了较厚的白垩系地层，形成了良好的生油层和储油层。而大野头凸起由于地势较高，沉积物堆积较少，且遭受了剥蚀作用，不利于油气的生成和保存。在古地热场方面，构造活动导致岩石的变形和破裂，改变了岩石的热导率和热传递路径，从而影响了古地热场的分布。断裂的存在为深部热流的向上运移提供了通道，使得断裂附近地区的地温相对较高。例如，[学者姓名7]通过热传导数值模拟研究发现，在靠近断层的区域，由于热流的集中运移，地温梯度明显升高，这种古地热场的不均匀分布对盆地内有机质的热演化和油气生成具有重要影响。褶皱构造也会影响热传递，褶皱的紧闭程度和轴向方向会改变岩石的热传导方向和效率，进而影响古地热场的分布格局。综上所述，胶莱盆地内部构造格局复杂，次级凹陷、凸起以及褶皱、断层等构造的发育，对盆地内的沉积和古地热场产生了重要影响。深入研究这些构造特征及其对沉积和地热场的控制作用，对于理解胶莱盆地的地质演化历史和能源资源形成机制具有重要意义。2.2地层发育特征2.2.1地层序列及岩性组合胶莱盆地地层主要由白垩系组成，其沉积序列记录了盆地复杂的地质演化历史，对研究盆地的形成和发展具有重要意义。自下而上，白垩系可划分为莱阳群、青山群和王氏群，各群组具有独特的岩性组合和沉积环境特征。莱阳群作为盆地内最古老的白垩系地层，主要形成于早白垩世早期。该群岩性组合以陆相碎屑岩为主，自下而上可进一步细分为多个组。底部的逍仙庄组主要由砾岩、含砾砂岩组成，砾石成分复杂，包括石英岩、花岗岩、变质岩等，磨圆度较差，分选性也不佳，显示出近源快速堆积的特征，反映了当时盆地处于强烈的构造活动期，物源区距离较近，搬运距离短，沉积环境较为动荡。向上的止凤庄组以砂岩、粉砂岩为主，夹少量泥岩，发育交错层理、波痕等沉积构造，表明沉积环境为河流-三角洲相，水体能量适中，水流方向多变。马耳山组则以泥岩、页岩为主，夹薄层砂岩，富含植物化石和淡水双壳类化石，反映了当时为相对稳定的湖泊相沉积环境，水体较深，气候湿润，生物繁盛。水南组岩性以砂岩、粉砂岩为主，发育平行层理和小型交错层理，显示沉积环境为滨浅湖相，水体能量相对较低。曲格庄组主要由泥岩、页岩组成，夹薄层砂岩，含有机质丰富，是良好的烃源岩层，其沉积环境为半深湖-深湖相，水体安静，有利于有机质的保存和富集。最上部的龙旺庄组以砂岩、砾岩为主，夹少量泥岩，沉积构造主要为大型交错层理，表明沉积环境又转变为河流相，盆地水体变浅，构造活动增强。青山群形成于早白垩世晚期，是一套以火山岩为主的地层。其岩性组合主要包括火山熔岩和火山碎屑岩，如安山岩、英安岩、流纹岩以及凝灰岩、火山角砾岩等。这些火山岩的喷发反映了当时盆地处于强烈的火山活动期，构造运动剧烈，深部岩浆沿断裂上涌喷发至地表。安山岩和英安岩具有斑状结构，斑晶主要为斜长石、角闪石等，基质为隐晶质或玻璃质，显示出岩浆快速冷凝的特征。流纹岩则具有流纹构造，矿物结晶程度较好，反映了岩浆在相对较高的温度和压力条件下缓慢流动和冷凝。火山碎屑岩中，凝灰岩粒度较细，成分以火山灰为主，常含有晶屑、玻屑等；火山角砾岩粒度较粗，角砾成分复杂，包括火山岩、围岩等，分选性和磨圆度均较差，是火山喷发时近距离堆积的产物。青山群的火山活动对盆地的沉积环境和地层演化产生了重大影响，火山物质的大量堆积改变了盆地的地形和水体条件，同时也为后续的沉积作用提供了新的物源。王氏群形成于晚白垩世，是胶莱盆地白垩系的最上部地层。其岩性组合以河湖相红色碎屑岩为主，包括砂岩、泥岩、砾岩等。底部常发育一套底砾岩，砾石成分主要为石英岩、燧石等，磨圆度较好，分选性中等，表明物源区相对较远，经过了一定距离的搬运。向上砂岩和泥岩互层，砂岩以中粗粒为主，发育交错层理、平行层理等沉积构造，反映了河流相的沉积环境，水流能量较强。泥岩则以紫红色为主，含有石膏、岩盐等蒸发矿物，表明当时气候干旱，湖水蒸发强烈，水体盐度较高，沉积环境为干旱的湖泊相。王氏群的沉积特征表明，晚白垩世时期，胶莱盆地的构造活动相对减弱，沉积环境逐渐由湿润向干旱转变，气候条件对沉积作用的影响更为显著。综上所述，胶莱盆地白垩系地层序列完整，岩性组合复杂多样，不同群组的沉积特征反映了盆地在不同地质时期的构造运动、沉积环境和古气候条件的变化，为研究盆地的地质演化提供了丰富的信息。2.2.2地层沉积演化过程胶莱盆地地层的沉积演化过程是多种地质因素相互作用的结果，与区域构造运动密切相关，经历了多个阶段，每个阶段都有其独特的沉积特征和演化规律。早白垩世早期，太平洋板块向欧亚板块俯冲，导致胶莱盆地所在区域地壳发生伸展减薄，形成一系列断陷盆地，莱阳群开始沉积。在这一时期，盆地边缘断裂活动强烈，地形高差大，物源区岩石遭受快速剥蚀，大量粗碎屑物质被搬运至盆地内，在盆地边缘形成了以砾岩、含砾砂岩为主的逍仙庄组沉积。随着盆地的进一步沉降，水体逐渐加深，沉积环境转变为河流-三角洲相和湖泊相，止凤庄组、马耳山组、水南组、曲格庄组和龙旺庄组依次沉积。这一阶段，沉积作用受构造运动和古气候的双重控制，构造运动控制了盆地的沉降速率和地形地貌，古气候则影响了物源区的风化程度和植被覆盖情况，进而影响了沉积物的类型和性质。例如，在气候湿润时期，物源区风化作用强烈，植被茂盛，土壤侵蚀相对较弱，河流搬运的碎屑物质粒度较细，有利于形成以砂岩、泥岩为主的沉积地层；而在气候干旱时期，物源区风化作用相对较弱，植被稀疏，土壤侵蚀强烈，河流搬运的碎屑物质粒度较粗，有利于形成以砾岩、砂岩为主的沉积地层。早白垩世晚期，区域构造应力场发生转变，太平洋板块俯冲作用加剧，导致深部岩浆活动强烈，青山群火山岩开始喷发。火山活动使得盆地内的沉积环境发生了巨大变化，大量火山物质堆积在盆地内，改变了盆地的地形和水体条件。火山岩的喷发还对生物群落产生了重大影响，许多生物在火山喷发过程中灭绝，同时也为新的生物群落的形成创造了条件。在这一时期，盆地内的沉积作用以火山物质的堆积为主，火山熔岩和火山碎屑岩广泛分布，形成了青山群独特的岩性组合。火山活动的间歇期，也有少量正常沉积作用发生，形成了一些火山碎屑岩与正常沉积岩的互层。晚白垩世，太平洋板块俯冲作用相对减弱，胶莱盆地进入坳陷阶段，王氏群开始沉积。此时，盆地构造活动相对稳定，地形高差减小，物源区距离较远，沉积物搬运距离增大。在干旱的气候条件下，盆地内形成了以河湖相红色碎屑岩为主的沉积地层。河流携带的碎屑物质在盆地内逐渐堆积，形成了砂岩、泥岩互层的沉积组合。由于气候干旱，湖水蒸发强烈，水体盐度升高，在泥岩中形成了石膏、岩盐等蒸发矿物。王氏群的沉积特征表明，这一时期盆地的沉积作用主要受古气候和物源区的控制，构造运动对沉积作用的影响相对减弱。新生代以来，受喜马拉雅运动的影响，胶莱盆地整体抬升，遭受剥蚀，沉积作用基本停止。地层遭受不同程度的侵蚀，盆地边缘地层出露较多，而盆地中心部分地层保存相对较好。在抬升过程中，盆地内的地层发生了褶皱和断裂，进一步改变了地层的形态和分布格局。例如，一些原本水平的地层被褶皱成背斜和向斜，断裂则导致地层的错动和位移，使得地层的连续性受到破坏。这些构造变形对盆地内的油气运移和聚集产生了重要影响，也为研究盆地的构造演化提供了重要线索。综上所述，胶莱盆地地层的沉积演化过程与区域构造运动密切相关，经历了断陷、火山活动、坳陷和抬升剥蚀等多个阶段，每个阶段的沉积特征反映了当时的构造背景、沉积环境和古气候条件。通过对地层沉积演化过程的研究，可以深入了解盆地的地质演化历史，为研究区域地质构造和资源分布提供重要依据。2.3岩石学特征2.3.1主要岩石类型及特征胶莱盆地内岩石类型丰富多样，主要包括页岩、砂岩、砾岩和火山岩等，这些岩石类型记录了盆地复杂的地质演化历史，对研究盆地的沉积环境和构造背景具有重要意义。页岩在胶莱盆地白垩系地层中广泛分布，尤其是在莱阳群的曲格庄组和马耳山组。页岩颜色多为黑色、灰黑色，质地细腻，页理发育，具有良好的薄层状构造。其矿物组成主要以黏土矿物为主，如伊利石、蒙脱石和高岭石等，这些黏土矿物在沉积过程中，由于水体环境的相对稳定，细小颗粒得以缓慢沉降并有序排列，形成了页岩独特的页理构造。页岩中还含有一定量的石英、长石等碎屑矿物，以及少量的黄铁矿、方解石等自生矿物。石英和长石碎屑主要来源于物源区岩石的风化剥蚀，经河流等搬运作用带入盆地沉积；黄铁矿的存在表明沉积环境为还原环境，有利于有机质的保存；方解石则可能是在成岩过程中，由孔隙水中的碳酸钙沉淀形成。页岩的结构为泥质结构，颗粒细小，分选性好，磨圆度差，这是由于黏土矿物在悬浮状态下缓慢沉积，没有经过长距离的搬运和强烈的水流分选。页岩中发育水平层理和微细层理，反映了沉积时水体平静、能量较低的环境，如在曲格庄组页岩中，水平层理清晰可见，表明当时为半深湖-深湖相沉积环境，水体安静，有利于有机质的大量堆积和保存。砂岩也是盆地内常见的岩石类型，在莱阳群、青山群和王氏群中均有分布。砂岩颜色多样，有灰白色、浅黄色、紫红色等，其颜色主要取决于碎屑颗粒的成分和胶结物的性质。矿物组成以石英、长石为主，石英含量通常在50%以上，长石含量次之，此外还含有少量的云母、岩屑等。石英和长石碎屑的来源与物源区的岩石类型密切相关，如物源区为花岗岩等酸性侵入岩时，砂岩中长石含量相对较高；若物源区为变质岩或沉积岩，则石英含量可能更高。砂岩的结构为碎屑结构，根据碎屑颗粒大小可分为粗砂岩、中砂岩和细砂岩。粗砂岩颗粒直径一般大于0.5mm，中砂岩颗粒直径在0.25-0.5mm之间，细砂岩颗粒直径小于0.25mm。砂岩的分选性和磨圆度因沉积环境而异，在河流相沉积中，砂岩分选性中等-好，磨圆度次棱角状-次圆状，这是由于河流搬运过程中，颗粒之间相互碰撞、摩擦，使得颗粒逐渐磨圆，同时水流的分选作用也使得颗粒按大小进行一定程度的分选；而在三角洲相沉积中，砂岩分选性较好，磨圆度多为次圆状-圆状，因为三角洲地区水流能量相对稳定，搬运距离适中，颗粒有足够的时间被磨圆和分选。砂岩中发育多种沉积构造，如交错层理、平行层理、波痕等。交错层理在河流相和三角洲相砂岩中较为常见，反映了水流方向的变化和水动力条件的强弱交替，如在王氏群的砂岩中，大型交错层理发育，表明当时沉积环境为河流相，水流能量较强，且存在侧向迁移；平行层理则常见于滨岸相砂岩中，反映了水动力较强且稳定的沉积环境，如在青山群部分砂岩中，平行层理清晰，说明沉积时水体能量较高且较为稳定。砾岩主要分布在莱阳群的底部和王氏群的部分层位。砾岩颜色多为灰色、灰白色，由砾石和填隙物组成。砾石成分复杂，包括石英岩、花岗岩、变质岩等，这些砾石来源于盆地周边的隆起区，是物源区岩石遭受强烈风化剥蚀后，经短距离搬运快速堆积形成的。砾石的磨圆度和分选性较差，多呈棱角状-次棱角状，这是因为砾石搬运距离短，没有经过充分的磨圆和分选作用。填隙物主要为砂、粉砂和黏土，胶结类型多为基底式胶结和孔隙式胶结。基底式胶结中，填隙物含量较多，砾石被填隙物完全包裹，表明沉积时水流能量较高，砾石快速堆积，填隙物随即填充其中；孔隙式胶结中，填隙物较少，主要填充在砾石的孔隙中，说明沉积后有一定的成岩作用，孔隙中的填隙物是后期由地下水带入沉淀形成的。砾岩中常见大型交错层理和块状层理，大型交错层理反映了水流能量较强且不稳定，如在莱阳群底部的砾岩中，大型交错层理发育，表明当时为近源快速堆积的扇三角洲相沉积环境，水流湍急，携带大量砾石快速沉积；块状层理则表明沉积过程中水流能量相对稳定，沉积物快速堆积，没有明显的层理分异。火山岩主要出现在青山群，是盆地早白垩世晚期强烈火山活动的产物。火山岩类型主要有安山岩、英安岩和流纹岩等。安山岩呈灰绿色、紫红色，具斑状结构，斑晶主要为斜长石、角闪石等，基质为隐晶质或玻璃质。斜长石斑晶常呈板状，表面可见聚片双晶，反映了岩浆在上升过程中的结晶分异作用；角闪石斑晶呈长柱状，多色性明显，其存在表明岩浆源区含有一定量的幔源物质。英安岩颜色较浅，常为灰白色、浅肉红色，斑状结构发育，斑晶主要为石英、长石，基质为隐晶质或细粒结构。石英斑晶呈他形粒状，表面干净，具波状消光，反映了其在岩浆中相对较晚结晶；长石斑晶以钾长石和斜长石为主，常具环带结构，是岩浆结晶过程中物理化学条件变化的结果。流纹岩呈灰白色、浅粉色，具有典型的流纹构造，矿物结晶程度较好，主要矿物为石英、钾长石和酸性斜长石。流纹构造是岩浆在流动过程中，不同成分和颜色的物质呈条带状排列形成的，反映了岩浆在相对较高的温度和压力条件下缓慢流动和冷凝。火山岩中还常发育气孔构造和杏仁构造，气孔构造是岩浆喷发时，气体逸出留下的空洞，气孔大小不一，形状不规则；杏仁构造则是气孔被后期矿物质充填形成的，常见的充填矿物有方解石、绿泥石等，这些矿物的充填反映了火山岩形成后的热液活动。综上所述，胶莱盆地内不同岩石类型具有各自独特的矿物组成、结构构造和沉积特征，这些特征反映了盆地在不同地质时期的沉积环境、构造背景和火山活动等信息，为研究盆地的地质演化提供了重要的岩石学依据。2.3.2岩石热物理性质岩石热物理性质是古地热场数值模拟的关键参数，其对热量在岩石中的传导、储存和散失过程起着决定性作用，进而影响古地热场的分布和演化。在胶莱盆地古地热场研究中，准确获取岩石的热导率、比热容等热物理性质参数至关重要。岩石热导率是衡量岩石传导热量能力的重要参数，它决定了热量在岩石中传递的速度和效率。不同岩石类型的热导率存在显著差异，这主要取决于岩石的矿物组成、结构构造以及孔隙度等因素。一般来说，矿物的热导率与其晶体结构和化学键强度密切相关。例如，石英是常见的造岩矿物，具有较高的热导率，其晶体结构紧密，化学键强度较大，有利于热量的传导；而黏土矿物如伊利石、蒙脱石等，热导率相对较低，这是因为它们的晶体结构较为松散，层间存在大量的水分子，阻碍了热量的传递。在胶莱盆地的岩石中，砂岩由于主要矿物为石英和长石，且颗粒之间接触紧密，其热导率相对较高，一般在2-5W/(m・K)之间。其中，石英含量较高的砂岩热导率更接近上限，如在一些石英砂岩中，热导率可达4-5W/(m・K)；而长石含量较高的砂岩热导率则相对较低，约为2-3W/(m・K)。页岩的热导率较低，通常在1-2W/(m・K)之间，这主要是由于其富含黏土矿物，且页理发育，层间结构疏松，不利于热量的传导。火山岩的热导率则因其矿物组成和结构的不同而有所变化，安山岩和英安岩的热导率一般在2-4W/(m・K)之间，流纹岩的热导率相对较低，约为1.5-3W/(m・K)。火山岩中的气孔构造和杏仁构造会降低岩石的整体热导率，因为气体和充填矿物的热导率远低于岩石基质。岩石比热容是指单位质量的岩石温度升高1℃所吸收的热量，它反映了岩石储存热量的能力。比热容的大小主要取决于岩石的化学成分和矿物组成。一般来说，岩石中所含矿物的比热容越大，岩石的比热容也越大。例如，水的比热容较大，因此富含水的岩石比热容相对较高。在胶莱盆地的岩石中，页岩由于含有一定量的黏土矿物和吸附水，其比热容相对较大，一般在0.8-1.2J/(g・K)之间。黏土矿物的晶体结构中存在大量的层间水和吸附水，这些水分子在温度变化时会吸收或释放大量的热量，从而使页岩具有较高的比热容。砂岩的比热容相对较小，约为0.7-1.0J/(g・K)，这是因为砂岩主要由石英和长石等矿物组成，这些矿物的比热容相对较低。火山岩的比热容因矿物组成和结构的不同而有所差异，安山岩和英安岩的比热容一般在0.8-1.1J/(g・K)之间，流纹岩的比热容约为0.7-1.0J/(g・K)。火山岩中的玻璃质成分和气孔构造会影响其比热容，玻璃质的比热容相对较低，而气孔中的气体则会降低岩石的整体比热容。在古地热场模拟中，岩石热物理性质起着至关重要的作用。热导率决定了热量在岩石中的传导速度，影响着古地温梯度的分布。在热导率较高的岩石中，热量能够快速传导，使得地温梯度相对较小；而在热导率较低的岩石中，热量传导缓慢，容易导致地温梯度增大。例如，在胶莱盆地中，砂岩的热导率较高，热量能够较快地在砂岩中传导，使得砂岩地层中的地温梯度相对较小；而页岩的热导率较低，热量在页岩中传导较慢，容易在页岩层中形成较高的地温梯度，这对有机质的热演化和油气生成具有重要影响。比热容则影响着岩石储存热量的能力，进而影响古地热场的演化过程。在构造运动或岩浆活动等热源作用下，比热容较大的岩石能够吸收更多的热量，减缓温度的升高速度；而比热容较小的岩石则容易升温，导致局部地区温度变化较快。例如，在岩浆侵入过程中，页岩由于比热容较大，能够吸收岩浆带来的部分热量，使得岩浆侵入区域周围的温度升高相对缓慢；而砂岩由于比热容较小，在岩浆热作用下，温度升高较快，这可能导致砂岩中矿物的重结晶和岩石结构的改变。综上所述，岩石热物理性质在胶莱盆地古地热场模拟中具有重要作用，不同岩石类型的热导率和比热容差异显著，这些差异对古地热场的分布和演化产生了重要影响。准确测定和分析岩石热物理性质参数，是建立准确的古地热场数值模型、深入研究胶莱盆地古地热场演化特征的关键。三、古地热场研究方法与数据获取3.1古地热场研究方法概述3.1.1“将今论古”法原理与应用“将今论古”法，又称历史比较法，是古地热场研究中的重要方法之一。其基本原理基于地质作用的相对稳定性和连续性，认为发生在地质历史时期的地质作用及其结果，与现代正在进行的地质作用及其产物具有相似性。通过对现代地质作用过程和产物的研究，总结出相关规律，进而运用这些规律来分析保留在地层及岩石中的各种地质现象，以此推断古代地质作用的过程和古地理环境，包括古地热场的特征。这一原理由英国地质学家莱伊尔在赫顿的均变论学说基础上提出，在地质学研究中具有重要地位。在胶莱盆地古地热场研究中，“将今论古”法具有一定的可行性。从地质构造角度来看，胶莱盆地现今的构造格局是长期地质演化的结果，其内部的断裂、褶皱等构造特征在一定程度上反映了地质历史时期的构造应力场变化。例如，通过对现代构造活动的研究，发现盆地内一些主要断裂现今仍有微弱活动，这些断裂在地质历史时期可能也是热液活动和热量传输的通道。因此，可以根据现今断裂附近的热异常现象，推断在过去地质时期，这些断裂同样可能对古地热场产生重要影响，导致古地温的局部升高或热流的异常分布。在沉积作用方面，现代沉积过程中不同沉积环境下形成的沉积物具有特定的特征，这些特征与沉积时的水动力条件、物源区性质等密切相关。在胶莱盆地中，通过对现代河流、湖泊等沉积环境下沉积物的研究，总结出不同沉积相的特征和分布规律。然后，将这些规律应用于对盆地内白垩系地层沉积相的分析，推断出不同地质时期的沉积环境和古地理格局。由于沉积环境与古地热场存在相互影响关系，稳定的沉积环境有利于热量的积累和保存，而快速的沉积速率可能导致沉积物的压实和热传导变化。因此，通过对沉积相的分析，可以间接了解古地热场在不同沉积时期的大致特征和变化趋势。然而，“将今论古”法在胶莱盆地应用时也存在一定局限性。胶莱盆地经历了复杂的地质演化历史，在漫长的地质时期内，地质条件发生了显著变化。从构造运动角度来看，中生代时期太平洋板块向欧亚板块的强烈俯冲，导致胶莱盆地经历了强烈的构造变形和岩浆活动，这一时期的构造应力场和深部热动力条件与现代有很大差异。在中生代强烈的构造运动期间，大量岩浆侵入盆地，带来了巨大的热量，使古地热场发生了剧烈变化。而现代胶莱盆地的构造活动相对较弱，岩浆活动基本停止，无法直接用现代构造活动和热流条件来推断中生代时期的古地热场特征。从沉积环境角度来看，地质历史时期的古气候和海平面变化与现代也存在明显不同。在白垩纪时期，全球气候经历了多次冷暖交替，胶莱盆地的沉积环境也随之发生了显著变化。在气候温暖湿润时期，盆地内可能形成大面积的湖泊，沉积物以细粒的泥岩和页岩为主；而在气候干旱时期，沉积环境可能转变为河流相或三角洲相，沉积物以粗粒的砂岩和砾岩为主。这些不同的沉积环境对古地热场的影响机制较为复杂，且与现代沉积环境下的热传导和热量分布规律存在差异，因此仅依靠“将今论古”法难以准确重建古地热场在不同古气候条件下的变化情况。综上所述，“将今论古”法在胶莱盆地古地热场研究中具有一定的应用价值，能够为古地热场的初步分析提供重要线索。但由于盆地地质演化的复杂性和地质条件的时空变化，在应用该方法时需要充分考虑其局限性，结合其他研究方法和多源数据，综合分析和推断古地热场的特征和演化历史。3.1.2地质温度计法及常用指标地质温度计法是通过测定沉积物中某些对温度变化敏感的指标，来推断古地温的方法。这些指标如同记录古地温信息的“温度计”，能够反映地质历史时期岩石所经历的最高温度或温度变化范围。在胶莱盆地古地温研究中，镜质体反射率和流体包裹体是常用的地质温度计指标。镜质体反射率是目前应用最广泛的有机质成熟度指标之一，也是推断古地温的重要参数。镜质体是由高等植物的根、茎、叶在覆水的还原条件下，经凝胶化作用而形成的。随着地层埋深的增加和温度的升高，镜质体经历了一系列的物理化学变化，其反射率也随之增大。镜质体反射率与古地温之间存在密切的正相关关系，一般来说，古地温越高，镜质体反射率越大。这是因为在较高的温度下，镜质体中的有机质发生热降解和缩聚反应，芳香结构逐渐增多，分子排列更加有序，从而导致反射率升高。在实际应用中，通过对胶莱盆地内钻孔岩芯或露头样品中的镜质体反射率进行测定，可以利用相关的经验公式或数学模型，反演古地温。例如，常用的Lopatin法和Waples法，通过考虑镜质体反射率随时间和温度的变化关系，计算出温度-时间指数（TTI），进而估算古地温。然而，镜质体反射率的应用也存在一定局限性。在石炭系以前的地层中，由于高等植物尚未繁盛，镜质体匮乏，限制了该方法的应用。在一些特殊岩性如缺乏高等植物输入源的海相碳酸盐岩中，镜质体也十分贫乏。此外，镜质体的反射率还受到其本身的生源构成、有机质类型以及沉积环境等因素的影响。富氢的镜质体反射率偏低，在低煤化阶段伴有较弱的棕褐色荧光，这是因为其形成过程中有低等生物和其他类脂化合物的参与。再循环（或再沉积）镜质体来自剥蚀区搬运而来的较老地层有机质中，不能代表该深度处有机质的真实成熟度，会对反射率测值产生干扰。流体包裹体是矿物结晶过程中捕获的成矿流体，它们被封闭在矿物晶格中，记录了成矿时的温度、压力、成分等信息。流体包裹体的均一温度是指包裹体受热后，其中的气相和液相完全混为一相时的温度，被认为是矿物形成时的最低温度，可用于推断古地温。在胶莱盆地中，通过对砂岩、页岩等岩石中矿物（如石英、方解石等）内的流体包裹体进行研究，测定其均一温度，能够获取古地温信息。当矿物结晶时，捕获的流体包裹体处于高温高压状态，随着温度和压力的降低，流体包裹体中的气相和液相会发生分离。通过在实验室中对流体包裹体进行加热，观察其均一化过程，测定均一温度，即可推断矿物形成时的古地温。流体包裹体均一温度的测定结果受到多种因素的影响。包裹体在形成后可能会受到后期构造运动、热液活动等的改造，导致其内部压力和温度发生变化，从而使均一温度不能准确反映原始成矿温度。不同类型的流体包裹体（如原生包裹体、次生包裹体等）其均一温度所代表的地质意义不同，需要准确识别和区分。此外，在测定过程中，实验条件的控制和仪器的精度也会对均一温度的测定结果产生影响。综上所述，镜质体反射率和流体包裹体作为地质温度计法的常用指标，在胶莱盆地古地温研究中具有重要应用价值。但它们各自存在一定的局限性，在实际研究中需要综合考虑多种因素，结合其他地质资料和研究方法，以提高古地温推断的准确性和可靠性。3.2数据获取与分析3.2.1野外地质调查与样品采集野外地质调查是研究胶莱盆地古地热场的基础工作，通过系统的调查和样品采集，能够获取丰富的地质信息，为后续的室内实验分析和数值模拟提供关键数据。调查区域覆盖了胶莱盆地的主要构造单元，包括莱阳凹陷、高密凹陷以及大野头凸起等。在莱阳凹陷，重点对莱阳群和青山群地层的露头进行观测，详细记录地层的岩性、厚度、接触关系以及沉积构造等特征。莱阳群底部的砾岩和含砾砂岩，通过观察其砾石成分、磨圆度和分选性，推断物源区的位置和沉积环境的能量条件。在高密凹陷，对王氏群地层进行了深入调查，关注其红色碎屑岩的沉积特征和古生物化石分布情况，以了解晚白垩世时期的沉积环境和古气候条件。大野头凸起由于出露地层相对较老，主要对前寒武纪变质岩系进行了观察，分析其岩石变形和变质程度，为研究盆地深部地质构造提供线索。在样品采集方面，根据不同的研究目的，在盆地内多个关键位置采集了岩石样品。在莱阳凹陷的莱阳群曲格庄组，采集了大量页岩样品，这些页岩富含有机质，是研究有机质热演化和古地温的重要样品。在青山群火山岩出露区，采集了安山岩、英安岩和流纹岩等火山岩样品，用于分析火山岩的矿物组成、结构构造以及岩石热物理性质。在王氏群地层，采集了砂岩和泥岩样品，用于研究沉积岩的岩石学特征和古地磁性质。样品采集过程中，严格按照相关规范进行操作，确保样品的代表性和完整性。对每个样品进行详细编号，记录采集位置、地层信息和样品特征等，为后续实验分析提供准确的背景资料。通过野外地质调查和样品采集，获取了丰富的地质信息和岩石样品。这些数据和样品为深入研究胶莱盆地的地质构造、地层发育和岩石学特征提供了第一手资料，为室内实验分析和古地热场数值模拟奠定了坚实的基础。通过对野外地质现象的观察和分析，初步了解了盆地内不同地层的沉积环境和构造背景，为后续研究提供了重要的地质依据。采集的岩石样品将在实验室中进行多种测试分析，以获取岩石的热物理性质、有机质成熟度等关键参数，进一步揭示胶莱盆地古地热场的演化特征。3.2.2室内实验分析对采集的岩石样品进行了一系列室内实验分析，旨在获取岩石的热物理性质、有机质成熟度等关键参数，为古地热场研究提供重要数据支持。镜质体反射率测试是确定有机质成熟度和推断古地温的重要方法之一。在实验过程中，首先对样品进行精心制备，将岩石样品切割成厚度均匀的薄片，确保镜质体在显微镜下能够清晰观察。利用高倍显微镜，选择无结构镜质体中最为常见的均质镜质体和基质镜质体两种亚组分进行反射率测定。采用国际标准的反射率测定方法，确保测试结果的准确性和可靠性。对莱阳凹陷莱阳群曲格庄组页岩样品的镜质体反射率测试结果显示，其反射率值在一定范围内变化，反映了该地层在地质历史时期经历的不同热演化程度。通过与相关经验公式和模型相结合，利用镜质体反射率数据初步估算了该地区的古地温，为研究古地热场的演化提供了重要依据。然而，在测试过程中也发现了一些问题，如同一样品中镜质体反射率测值存在一定的离散程度，这可能与镜质体本身的生源构成、有机质类型以及沉积环境等因素有关。富氢的镜质体反射率偏低，在低煤化阶段伴有较弱的棕褐色荧光，这是由于其形成过程中有低等生物和其他类脂化合物的参与。再循环（或再沉积）镜质体来自剥蚀区搬运而来的较老地层有机质中，不能代表该深度处有机质的真实成熟度，会对反射率测值产生干扰。流体包裹体分析是另一种重要的古地温研究方法。对砂岩、页岩等岩石中矿物（如石英、方解石等）内的流体包裹体进行了详细研究。在显微镜下，仔细识别和挑选原生流体包裹体，因为原生包裹体能够更准确地反映矿物形成时的物理化学条件。利用冷热台对流体包裹体进行加热，观察其均一化过程，测定均一温度。通过对多个样品中流体包裹体均一温度的测定，获得了不同地层在不同地质时期的古地温信息。在青山群火山岩相关样品中，流体包裹体均一温度的测定结果显示，部分包裹体的均一温度较高，这可能与火山活动期间的高温环境有关。然而，流体包裹体分析也存在一定的局限性，包裹体在形成后可能会受到后期构造运动、热液活动等的改造，导致其内部压力和温度发生变化，从而使均一温度不能准确反映原始成矿温度。不同类型的流体包裹体（如原生包裹体、次生包裹体等）其均一温度所代表的地质意义不同，需要准确识别和区分。在测定过程中，实验条件的控制和仪器的精度也会对均一温度的测定结果产生影响。岩石热物理性质测试也是室内实验分析的重要内容。通过实验测定了岩石的热导率、比热容等热物理性质参数。对于热导率的测定，采用稳态法或瞬态法，根据岩石样品的特点选择合适的测试方法。在测试过程中，严格控制实验条件，确保测试结果的准确性。对不同岩石类型的热导率测试结果表明，砂岩的热导率相对较高，一般在2-5W/(m・K)之间，这是由于其主要矿物为石英和长石，颗粒之间接触紧密，有利于热量的传导；而页岩的热导率较低，通常在1-2W/(m・K)之间，主要是因为其富含黏土矿物，且页理发育，层间结构疏松，阻碍了热量的传递。在比热容测试中，利用量热实验，准确测量单位质量岩石温度升高1℃所吸收的热量。页岩由于含有一定量的黏土矿物和吸附水，其比热容相对较大，一般在0.8-1.2J/(g・K)之间；砂岩的比热容相对较小，约为0.7-1.0J/(g・K)。这些岩石热物理性质参数的获取，为古地热场数值模拟提供了关键输入，能够更准确地模拟热量在岩石中的传导和储存过程。综上所述，通过室内实验分析，获取了胶莱盆地岩石的镜质体反射率、流体包裹体均一温度以及热物理性质等重要数据。这些数据为研究古地热场的演化提供了有力支持，但同时也需要认识到各种实验方法的局限性，在数据分析和解释过程中综合考虑多种因素，以提高研究结果的准确性和可靠性。3.2.3地球物理数据收集与处理地球物理数据在古地热场研究中具有重要作用，通过收集和处理区域重力、磁力等地球物理数据，能够获取深部地质构造信息，为古地热场研究提供重要约束。在重力数据收集方面，广泛收集了胶莱盆地及周边区域的重力测量数据，这些数据来源于不同时期的地质调查项目，涵盖了多种比例尺的测量成果。对收集到的重力数据进行了严格的质量筛选，去除了明显异常和误差较大的数据。利用专业的数据处理软件，对重力数据进行了地形校正、布格校正等预处理工作。地形校正的目的是消除地形起伏对重力测量结果的影响，因为地形的高低变化会导致重力值的虚假变化。布格校正则是为了消除地球内部密度不均匀和地球形状等因素对重力值的影响，使重力数据更能真实反映地下地质体的密度差异。经过校正后，绘制了高精度的重力异常图。从重力异常图上可以清晰地看出，胶莱盆地内存在明显的重力异常变化。在莱阳凹陷和高密凹陷等沉积中心区域，重力值相对较低，这是由于这些区域沉积地层较厚，岩石密度相对较小；而在大野头凸起等隆起区域，重力值相对较高，反映了其深部岩石密度较大。通过对重力异常的分析，初步推断了盆地内不同构造单元的深部结构和地层分布情况。重力异常的变化还与断裂构造有关，在一些重力梯度变化较大的区域，往往存在断裂构造，这些断裂可能对深部热流的运移和古地热场的分布产生重要影响。磁力数据的收集与处理同样重要。收集了区域内的航空磁力测量数据和地面磁力测量数据，对这些数据进行了系统的整理和分析。在数据处理过程中，首先进行了噪声滤波处理，去除了由于测量仪器噪声、电磁干扰等因素产生的异常信号。对磁力数据进行了日变校正，以消除地球磁场的日变化对测量结果的影响。通过一系列的数据处理，得到了准确的磁力异常图。从磁力异常图上可以发现，胶莱盆地内的磁力异常分布与地质构造密切相关。在青山群火山岩分布区域，由于火山岩具有较高的磁性，表现出明显的正磁异常。这些磁力异常信息有助于确定火山岩的分布范围和厚度变化，进而了解火山活动对古地热场的影响。在断裂构造附近，磁力异常也会出现明显的变化，这是因为断裂带两侧岩石的磁性差异以及断裂活动导致的岩石变形和磁化率变化。通过对磁力异常的分析，可以进一步识别和追踪盆地内的断裂构造，为研究古地热场的热传递通道提供重要线索。地球物理数据对古地热场研究具有多方面的作用。重力和磁力数据所反映的深部地质构造信息，为古地热场数值模拟提供了重要的边界条件和初始条件。通过对重力和磁力异常的分析，确定了盆地内不同构造单元的边界和深部结构，这些信息在数值模拟中用于设置岩石的热物理性质和热边界条件，能够更真实地模拟古地热场的分布和演化。地球物理数据可以帮助解释古地热场的异常现象。在某些区域，古地热场的异常可能与深部地质构造有关，通过地球物理数据的分析，可以确定这些构造因素对古地热场的影响机制。在重力异常低值区，可能存在深部热流上涌的现象，导致该区域古地温升高；而在磁力异常变化较大的区域，由于岩石磁性的变化可能影响热传导，进而影响古地热场的分布。地球物理数据还可以与其他地质数据（如地质调查数据、室内实验数据等）相结合，进行综合分析，提高古地热场研究的准确性和可靠性。将重力、磁力数据与岩石热物理性质数据相结合，可以更准确地推断深部岩石的热导率和比热容等参数，从而优化古地热场数值模型。综上所述，地球物理数据的收集与处理在胶莱盆地古地热场研究中具有重要意义。通过对重力、磁力数据的系统分析，获取了深部地质构造信息，为古地热场研究提供了重要的约束和支持。地球物理数据与其他地质数据的综合应用，有助于深入理解古地热场的分布特征和演化规律，为全面研究胶莱盆地的地质演化提供了新的视角和方法。四、胶莱盆地古地热场数值模拟4.1数值模拟原理与模型选择4.1.1传热基本原理与方程在胶莱盆地古地热场数值模拟中，热传导、对流和辐射是热量传递的三种基本方式，它们在盆地的地质演化过程中起着不同程度的作用。热传导是指热量通过物质分子的热运动，从高温区域向低温区域传递的过程，是固体内部热量传递的主要方式。在胶莱盆地的岩石中，热传导是热量传递的重要机制之一。根据傅里叶定律，热传导的基本公式为：q=-k\nablaT，其中q表示热流密度，单位为W/m^2；k为岩石的热导率，单位为W/(m·K)，它反映了岩石传导热量的能力，不同岩石类型的热导率差异较大，如前文所述，砂岩的热导率一般在2-5W/(m・K)之间，页岩的热导率通常在1-2W/(m・K)之间；\nablaT是温度梯度，单位为K/m，表示温度在空间上的变化率。该公式表明，热流密度与热导率和温度梯度成正比，热导率越高，温度梯度越大，热流密度就越大。对流是指流体（液体或气体）中由于温度差异引起的宏观运动，从而导致热量的传递。在胶莱盆地中，对流主要发生在地下流体（如地下水、油气等）中。当流体受热时，其密度减小，会向上运动；而较冷的流体则会向下运动，形成对流循环。对流过程中，热量随着流体的运动而被携带和传递。对流的强度与流体的性质、温度差以及流动通道的几何形状等因素有关。在盆地的沉积地层中，如果存在连通的孔隙或裂缝，地下水等流体就可以在其中流动，形成对流，从而影响古地热场的分布。例如，在一些断裂附近，由于岩石破碎，孔隙和裂缝发育，地下水的对流作用可能较为强烈，导致该区域的热量传递加快，古地温分布也会发生相应变化。辐射是指物体通过发射电磁波的方式传递热量的过程。在胶莱盆地的地质环境中，辐射的作用相对较小，但在某些特殊情况下，如高温岩浆侵入时，辐射传热可能会对周围岩石的温度分布产生一定影响。辐射传热的强度与物体的温度、发射率以及周围环境的辐射特性等因素有关。一般来说，温度越高的物体，辐射出的热量越多。在岩浆侵入地层时，高温岩浆会向周围岩石辐射热量，使周围岩石的温度升高，这种辐射传热在岩浆侵入初期可能是热量传递的重要方式之一，但随着时间的推移，热传导和对流的作用会逐渐增强。热传导方程是描述热传导过程中温度随时间和空间变化的偏微分方程。对于三维各向同性且均匀的介质，在有内热源的情况下，热传导方程的一般形式为：\rhoc\frac{\partialT}{\partialt}=\nabla\cdot(k\nablaT)+Q，其中\rho是岩石的密度，单位为kg/m^3，不同岩石的密度也有所差异，如砂岩的密度一般在2.5-2.7kg/m^3左右，页岩的密度约为2.4-2.6kg/m^3；c为比热容，单位为J/(kg·K)，前文已提及不同岩石比热容的范围；\frac{\partialT}{\partialt}表示温度对时间的变化率，单位为K/s；\nabla\cdot(k\nablaT)是热传导项，表示热量在空间中的传导，它考虑了热导率和温度梯度的影响；Q为内热源强度，单位为W/m^3，表示单位体积内热源产生的热量，在胶莱盆地中，岩浆活动、放射性元素衰变等都可能成为内热源。如果介质中没有内热源，即Q=0，则方程简化为齐次热传导方程：\rhoc\frac{\partialT}{\partialt}=\nabla\cdot(k\nablaT)。在数值模拟中，需要确定初始条件和边界条件来求解热传导方程。初始条件是指模拟开始时的温度分布，通常根据研究区域的地质背景和已有数据来确定。在胶莱盆地古地热场模拟中，可以利用地质调查获取的地层温度信息、古地温研究得到的古地温数据等，结合盆地的地质演化历史，合理设定初始温度场。边界条件则描述了研究区域边界上的温度或热流状况。常见的边界条件有三种类型：第一类边界条件（Dirichlet边界条件），给定边界上的温度值，即T(x,y,z,t)|_{\Gamma}=T_0(x,y,z,t)，其中\Gamma表示边界，T_0(x,y,z,t)是已知的边界温度函数；第二类边界条件（Neumann边界条件），给定边界上的热流密度值，即k\frac{\partialT}{\partialn}|_{\Gamma}=q_0(x,y,z,t)，其中\frac{\partialT}{\partialn}是温度沿边界外法线方向的导数，q_0(x,y,z,t)是已知的边界热流密度函数；第三类边界条件（Robin边界条件），给定边界上的热流密度与温度的线性关系，即(k\frac{\partialT}{\partialn}+\alphaT)|_{\Gamma}=g(x,y,z,t)，其中\alpha是与边界热交换相关的系数，g(x,y,z,t)是已知函数。在胶莱盆地的实际模拟中，需要根据盆地的边界特征和与周边地质体的热交换情况，选择合适的边界条件。例如，在盆地与周边山脉或其他地质体的接触边界上，可能存在热交换，此时可以采用第三类边界条件来描述边界上的热流和温度关系。4.1.2数值模拟方法与软件选择有限差分法和有限元法是求解热传导方程常用的数值模拟方法，它们各有特点，适用于不同的地质模型和研究需求。有限差分法是将求解区域划分为规则的网格，用差商来近似代替偏导数，将热传导方程离散化为代数方程组进行求解。在一维热传导问题中，对于热传导方程\frac{\partialT}{\partialt}=\alpha\frac{\partial^2T}{\partialx^2}（其中\alpha=\frac{k}{\rhoc}为热扩散率），采用显式差分格式，可将时间和空间进行离散。设时间步长为\Deltat，空间步长为\Deltax，在第n个时间步和第i个空间节点处的温度为T_{i}^n，则根据显式差分公式T_{i}^{n+1}=T_{i}^n+\frac{\alpha\Deltat}{(\Deltax)^2}(T_{i+1}^n-2T_{i}^n+T_{i-1}^n)，可以通过已知的当前时间步的温度值计算下一个时间步的温度。有限差分法的优点是计算简单、直观，易于编程实现，对于规则形状的求解区域和简单的边界条件，能够快速得到计算结果。在模拟简单的地层结构时，有限差分法可以快速计算出温度分布。但该方法对复杂边界条件的处理较为困难，当求解区域的边界不规则或存在复杂的地质构造时，网格划分会变得复杂，计算精度也会受到影响。有限元法是将求解区域划分为有限个单元，通过变分原理或加权余量法将热传导方程转化为一组线性代数方程组进行求解。在有限元法中，首先对求解区域进行离散化，将其划分为三角形、四边形等单元，然后在每个单元内构造插值函数来逼近温度分布。通过将单元的热传导方程进行组装，得到整个求解区域的方程组。有限元法的优点是对复杂边界条件和不规则求解区域具有很强的适应性，能够精确模拟各种复杂的地质构造和热边界条件。在模拟胶莱盆地这种具有复杂构造和边界条件的区域时，有限元法能够更好地处理地层的弯曲、断裂等情况，提高模拟的准确性。它还可以方便地处理不同岩石类型的热物理性质差异，通过在不同单元中赋予相应的参数，更真实地反映地质实际情况。然而，有限元法的计算量较大，对计算机的内存和计算速度要求较高，且前处理过程（如网格划分）相对复杂，需要一定的专业知识和技能。在本次胶莱盆地古地热场数值模拟中，选择COMSOLMultiphysics软件进行模拟。COMSOLMultiphysics是一款功能强大的多物理场耦合分析软件，具有丰富的物理模型库和高效的求解器。它能够方便地处理热传导、对流等多种物理过程的耦合问题，这对于模拟胶莱盆地复杂的地质环境非常重要。在模拟盆地内地下水对流与热传导的耦合时，COMSOLMultiphysics可以准确地考虑流体的流动和热量的传递相互作用。该软件具有友好的用户界面和强大的前后处理功能。在建模过程中，用户可以通过图形化界面直观地构建地质模型，设置材料参数和边界条件。在模拟完成后，软件提供了丰富的后处理工具，能够对模拟结果进行可视化分析，如绘制温度分布云图、温度随深度或时间的变化曲线等，方便研究人员直观地了解古地热场的分布和演化特征。COMSOLMultiphysics还支持并行计算，能够利用多核处理器提高计算效率，缩短模拟时间，这对于处理大规模的数值模拟问题至关重要。由于胶莱盆地古地热场模拟涉及到大量的网格单元和复杂的计算过程，并行计算功能可以大大提高模拟的效率，使研究人员能够在更短的时间内得到模拟结果。4.2模型参数设置4.2.1岩石热物理参数输入岩石热物理参数是古地热场数值模拟的关键输入，其准确性直接影响模拟结果的可靠性。在胶莱盆地古地热场模拟中，主要输入的岩石热物理参数包括热导率、比热容和密度等。热导率是反映岩石传导热量能力的重要参数。在胶莱盆地，不同岩石类型的热导率差异显著。通过对采集的岩石样品进行实验室测试，结合前人研究成果，确定了主要岩石类型的热导率取值范围。砂岩的热导率较高，一般在2-5W/(m・K)之间。其中，石英含量较高的砂岩，由于石英具有良好的热传导性能，其热导率更接近5W/(m・K)；而长石含量相对较多的砂岩，热导率则偏向2W/(m・K)。页岩的热导率较低，通常在1-2W/(m・K)之间。这是因为页岩主要由黏土矿物组成，黏土矿物的晶体结构疏松，层间存在大量吸附水，阻碍了热量的传导。火山岩的热导率因矿物组成和结构不同而有所变化，安山岩和英安岩的热导率一般在2-4W/(m・K)之间，流纹岩的热导率相对较低，约为1.5-3W/(m・K)。火山岩中的气孔构造和杏仁构造会降低岩石的整体热导率，因为气体和充填矿物的热导率远低于岩石基质。比热容是指单位质量的岩石温度升高1℃所吸收的热量，它反映了岩石储存热量的能力。页岩由于含有一定量的黏土矿物和吸附水，其比热容相对较大，一般在0.8-1.2J/(g・K)之间。黏土矿物的层间水和吸附水在温度变化时会吸收或释放大量热量，从而使页岩具有较高的比热容。砂岩的比热容相对较小，约为0.7-1.0J/(g・K)，这是因为砂岩主要由石英和长石等矿物组成，这些矿物的比热容相对较低。火山岩的比热容因矿物组成和结构的不同而有所差异，安山岩和英安岩的比热容一般在0.8-1.1J/(g・K)之间，流纹岩的比热容约为0.7-1.0J/(g・K)。火山岩中的玻璃质成分和气孔构造会影响其比热容，玻璃质的比热容相对较低，而气孔中的气体则会降低岩石的整体比热容。岩石密度也是重要的热物理参数之一，它影响着热传导方程中的惯性项和重力项。砂岩的密度一般在2.5-2.7kg/m^3左右，页岩的密度约为2.4-2.6kg/m^3，火山岩的密度因具体岩石类型而异，安山岩和英安岩的密度通常在2.7-2.9kg/m^3之间，流纹岩的密度相对较小，约为2.5-2.7kg/m^3。这些密度数据是通过对岩石样品的测量和分析得到的，在数值模拟中用于准确计算热传导方程中的各项参数。岩石热物理参数的不确定性对模拟结果有重要影响。热导率的不确定性会直接影响热量在岩石中的传导速度和分布。如果热导率取值偏高，模拟结果中热量传导速度会加快，导致地温梯度减小，古地温分布相对均匀；反之，如果热导率取值偏低，热量传导速度会减慢，地温梯度增大，可能会出现局部高温异常。在模拟页岩地层的古地热场时，若热导率取值比实际值高，会使页岩中热量快速传导，无法准确反映页岩由于热导率低而导致的热量积聚现象，从而影响对有机质热演化和油气生成的模拟结果。比热容的不确定性会影响岩石储存热量的能力，进而影响古地热场的演化过程。比热容取值偏大，岩石在受热时温度升高缓慢，会使模拟的古地温偏低；比热容取值偏小，岩石温度升高迅速，古地温会偏高。在岩浆侵入过程中，若比热容取值不准确，会导致对岩浆热影响范围和程度的模拟偏差，无法准确反映岩浆活动对古地热场的改造作用。综上所述，在胶莱盆地古地热场数值模拟中，准确输入岩石热物理参数至关重要。需要充分考虑不同岩石类型热物理参数的差异，并认识到参数不确定性对模拟结果的影响。通过多种方法获取准确的参数值，并进行敏感性分析，以提高模拟结果的可靠性和准确性。4.2.2边界条件与初始条件设定边界条件和初始条件的设定是古地热场数值模拟的重要环节，它们对模拟结果的准确性和可靠性起着关键作用。在胶莱盆地古地热场模拟中，需要合理设定盆地边界的热流、温度等边界条件以及初始地温场。对于边界条件的设定，主要考虑盆地与周边地质体的热交换情况。在盆地的底部边界，假设热流为一恒定值。这一假设的依据是基于区域地质研究，认为盆地底部的热流主要来源于地球深部的热传导，在较长的地质时期内相对稳定。通过参考区域大地热流数据以及周边类似盆地的研究成果，确定盆地底部热流值为60-80mW/m²。在盆地的侧向边界，考虑到与周边山脉和其他地质体的热交换，采用第三类边界条件。即给定边界上的热流密度与温度的线性关系：(k\frac{\partialT}{\partialn}+\alphaT)|_{\Gamma}=g(x,y,z,t)，其中\alpha是与边界热交换相关的系数，通过对盆地周边地质体的岩石热物理性质和热交换过程的分析，结合相关文献资料，确定\alpha的取值范围为0.01-0.1W/(m²・K)；g(x,y,z,t)是已知函数，根据盆地周边地质体的温度分布和热流情况确定其具体表达式。在盆地与周边山脉的接触边界，由于山脉岩石的热导率和热容量与盆地内岩石不同，会存在一定的热交换，通过这种边界条件的设定，可以更准确地模拟热量在边界处的传递和交换过程。在盆地的顶部边界，考虑到与大气的热交换以及地表温度的变化，采用与地表温度相关的边界条件。地表温度受到太阳辐射、大气环流和地形等多种因素的影响，在地质历史时期也会发生变化。通过对区域气候数据和地质历史时期古气候研究的综合分析，确定地表温度的变化范围。在新生代以前，根据古气候重建研究，地表温度相对较高，在白垩纪时期，平均地表温度可能在20-30℃之间；新生代以来，随着全球气候的变化，地表温度逐渐降低，现代地表温度平均约为10-15℃。在数值模拟中，将地表温度作为时间的函数输入，以反映其在地质历史时期的变化。考虑到大气与地表之间的热交换，采用热交换系数来描述这种热传递过程。通过对大气热物理性质和地表热交换过程的研究，确定热交换系数的取值，使得模拟能够更真实地反映盆地顶部边界的热状况。初始条件的设定是指模拟开始时盆地内的温度分布。在胶莱盆地古地热场模拟中，初始地温场的设定主要依据地质调查和古地温研究成果。通过对盆地内钻孔岩芯的温度测量以及镜质体反射率、流体包裹体等古地温指标的分析，结合地层沉积和构造演化历史，初步确定初始地温场。在盆地的沉积中心，由于沉积地层较厚，热量积聚，初始地温相对较高；而在盆地边缘和隆起区域，地层较薄，热量散失较快，初始地温相对较低。利用这些信息，采用插值方法构建初始地温场。以盆地内多个钻孔的温度数据为控制点，通过线性插值或样条插值等方法，在整个模拟区域内生成初始温度分布。在构建初始地温场时，还考虑了地层的岩性和热物理性质对温度分布的影响。对于热导率较高的砂岩地层，温度分布相对均匀；而对于热导率较低的页岩地层，可能存在温度梯度较大的区域。通过合理考虑这些因素，可以使初始地温场更符合实际地质情况。综上所述，边界条件和初始条件的设定需要综合考虑多种地质因素和研究成果。通过合理设定边界条件，能够准确描述盆地与周边地质体的热交换过程；通过科学设定初始条件，能够为古地热场模拟提供符合实际地质情况的初始状态。这些条件的准确设定是保证古地热场数值模拟结果准确性和可靠性的关键。4.3模拟过程与结果分析4.3.1模拟过程与步骤在利用COMSOLMultiphysics软件对胶莱盆地古地热场进行数值模拟时，首先进行了详细的模型构建。基于对胶莱盆地地质构