贵州省六盘水机场东北端高填方体地基稳定性多维度剖析与加固策略研究

贵州省六盘水机场东北端高填方体地基稳定性多维度剖析与加固策略研究一、绪论1.1研究背景与意义随着我国经济的飞速发展和综合国力的不断提升，航空运输业作为国家重要的基础设施和现代综合交通运输体系的重要组成部分，取得了令人瞩目的成就。机场作为航空运输的关键节点，其建设规模和数量持续增长。在西部山区，受复杂地形地貌的限制，“开山填谷”成为机场建设的常见方式，这不可避免地导致了大量高填方工程的出现，贵州省六盘水机场便是其中之一。六盘水机场位于贵州省六盘水市钟山区境内，是该地区重要的交通枢纽。机场东北侧存在高填方体作为路基，其填土路基的建设对机场地基稳定性有着至关重要的影响。高填方体地基稳定性直接关乎机场的安全运营，一旦出现问题，可能引发跑道变形、塌陷等严重后果，对飞机起降安全构成巨大威胁，进而影响整个机场的正常运行秩序，造成难以估量的经济损失和社会影响。从工程实践角度来看，深入研究六盘水机场东北端高填方体地基稳定性，能够为该机场的建设、维护和运营提供科学依据和技术支持。通过准确评估地基稳定性，合理设计地基处理方案和填方工程，可有效预防潜在的地基病害，确保机场设施的安全可靠，降低工程建设和运营成本。同时，研究成果对我国其他类似地形条件下的机场建设以及各类高填方工程具有重要的参考和借鉴价值，有助于推动我国基础设施建设领域的技术进步和工程质量提升。在学术理论层面，高填方体地基稳定性研究涉及土力学、岩石力学、工程地质学等多个学科领域，是一个复杂且具有挑战性的课题。对六盘水机场东北端高填方体地基稳定性的研究，能够丰富和完善相关学科的理论体系，为解决复杂地质条件下的地基稳定性问题提供新的思路和方法，促进学科交叉融合和创新发展。此外，研究过程中所采用的先进技术和方法，如数值模拟、现场监测等，也将为相关领域的科学研究提供有益的实践经验。1.2国内外研究现状1.2.1高填方稳定性研究现状高填方稳定性研究一直是岩土工程领域的重要课题，国内外学者在这方面开展了大量研究，取得了丰富的成果。在理论研究方面，极限平衡理论是高填方边坡稳定性分析中应用最早且最为广泛的理论之一。1916年瑞典人Petterson提出的瑞典圆弧滑动法，是条分法分析土质边坡稳定性中最古老、最简单的方法，该方法假设土质边坡的稳定性问题是平面应变问题，滑裂面为圆柱面，计算时不考虑各个土条之间的相互作用力，通过滑裂面上全部抗滑力与滑动力之比来定义土质边坡的安全系数。此后，Fellenius对条分法进一步改进，形成了极限平衡条分法的理论体系。在此基础上，国内外诸多学者通过提出各种不同假设，逐渐形成了适用范围不同的多种极限平衡条分法，如Bishop法、Janbu法、Lowe-Karafiath法、Morgenstern-Price法、Spencer法、美国陆军工团法、Sama法、不平衡推力传递法等。这些方法在不同程度上考虑了土条间的相互作用力，使计算结果更加符合实际情况。随着计算机技术的飞速发展，数值分析方法在高填方稳定性研究中得到了广泛应用。有限元法（FEM）是目前应用最为广泛的数值分析方法之一，它将连续体离散化为有限个单元，通过单元节点位移表示单元内部位移，建立整体刚度方程求解节点位移和应力，能够较好地模拟高填方的应力分布和变形情况，分析不同工况下的边坡稳定性。此外，有限差分法（FDM）、无界元法（IDEM）、数值流形元法（NMM）等数值方法也在高填方稳定性分析中发挥着重要作用。这些数值方法可以考虑土体的非线性、非均质、各向异性等复杂特性，以及边界条件和施工过程等因素对高填方稳定性的影响，为高填方工程的设计和分析提供了更加准确和全面的手段。在试验研究方面，室内试验和现场试验是获取高填方土体物理力学参数和验证理论分析与数值模拟结果的重要手段。室内试验可以对土体的基本物理力学性质进行测试，如密度、含水量、抗剪强度、压缩性等，为理论分析和数值模拟提供基础数据。现场试验则能够在实际工程条件下对高填方的变形和稳定性进行监测和研究，获取真实的工程数据，验证理论和数值模型的准确性和可靠性。例如，通过在高填方现场埋设各种监测仪器，如位移计、压力计、测斜仪等，实时监测高填方在施工过程和运营期间的变形、应力变化等情况，分析高填方的稳定性变化规律。尽管国内外在高填方稳定性研究方面取得了显著成果，但仍存在一些不足之处。例如，在理论分析方面，现有的各种分析方法都存在一定的假设条件和局限性，难以完全准确地描述高填方复杂的受力状态和破坏机制；在数值模拟方面，土体本构模型的选择和参数确定仍然是一个难题，不同的本构模型和参数对模拟结果的影响较大，且模拟过程中往往难以考虑所有的影响因素；在试验研究方面，室内试验和现场试验都存在一定的局限性，室内试验条件与实际工程情况存在差异，现场试验则受到各种条件的限制，难以获取全面和准确的数据。未来，高填方稳定性研究的发展趋势将是进一步完善理论分析方法，开发更加准确和实用的数值模型，加强试验研究与理论分析、数值模拟的结合，深入研究高填方在复杂地质条件和荷载作用下的稳定性问题，为高填方工程的设计、施工和运营提供更加可靠的技术支持。1.2.2高填方体变形研究现状高填方体变形研究对于保障工程的安全和正常使用具有重要意义，国内外学者在这一领域进行了深入的探索，取得了一系列成果。在变形机理研究方面，高填方体的变形主要包括瞬时变形、固结变形和次固结变形。瞬时变形是指在荷载施加瞬间，由于土体颗粒的重新排列而产生的变形，其变形速度快，主要受土体的弹性性质和初始应力状态影响。固结变形是由于土体孔隙中的水排出，土体发生压缩而产生的变形，这一过程与土体的渗透性、压缩性以及荷载大小和作用时间密切相关。次固结变形则是在固结变形基本完成后，由于土骨架的蠕变而产生的缓慢变形，通常在长期荷载作用下较为明显。影响高填方体变形的因素众多，主要包括土体性质、填方高度、地基条件、施工工艺和排水条件等。不同的土体性质，如颗粒级配、含水量、黏聚力和内摩擦角等，对高填方体的变形特性有着显著影响。一般来说，细粒土的压缩性较大，在相同荷载作用下产生的变形也较大；而粗粒土的透水性较好，固结速度相对较快，变形相对较小。填方高度越大，高填方体的自重应力越大，由此产生的变形也越大。地基条件是影响高填方体变形的重要因素之一，软弱地基会导致高填方体产生较大的沉降和不均匀变形。施工工艺的合理性直接关系到高填方体的压实质量，进而影响其变形特性。例如，分层填筑厚度过大、压实度不足等都会导致高填方体在后期使用过程中产生较大的变形。良好的排水条件能够加速土体孔隙水的排出，减少固结变形和次固结变形的发生。在变形预测方法方面，目前常用的方法主要有经验公式法、数值模拟法和现场监测法。经验公式法是根据大量的工程实践数据，建立起填方高度、土体性质等因素与变形之间的经验关系，从而对高填方体的变形进行预测。这种方法简单易行，但由于其建立在特定的工程条件和经验基础上，通用性和准确性相对有限。数值模拟法如有限元法、有限差分法等，通过建立高填方体的数值模型，考虑土体的力学性质、边界条件和施工过程等因素，对高填方体的变形进行模拟预测。数值模拟法能够较为全面地考虑各种因素对变形的影响，但模型的建立和参数的选取需要一定的经验和专业知识，且计算结果的准确性依赖于模型的合理性和参数的准确性。现场监测法是在高填方体施工和运营过程中，通过埋设各种监测仪器，如水准仪、全站仪、沉降计等，实时监测高填方体的变形情况，根据监测数据对变形进行分析和预测。现场监测法能够直接获取高填方体的实际变形数据，准确性高，但监测范围和时间有限，且监测成本较高。尽管国内外在高填方体变形研究方面取得了一定的进展，但仍存在一些问题有待解决。例如，目前对于高填方体在复杂应力状态和长期荷载作用下的变形机理研究还不够深入，变形预测方法的准确性和可靠性仍需进一步提高。此外，不同预测方法之间的比较和验证研究相对较少，缺乏统一的评价标准。未来，高填方体变形研究需要进一步加强理论研究，深入揭示变形机理，结合多种预测方法，提高变形预测的准确性和可靠性，同时加强现场监测数据的分析和应用，为高填方体工程的设计和施工提供更加科学的依据。1.2.3高填方软基处理技术高填方软基处理技术是确保高填方工程稳定性和正常使用的关键环节，针对不同的工程地质条件和工程要求，目前已发展出多种软基处理技术。常见的高填方软基处理技术包括换填法、强夯法、碎石桩法、CFG桩法、排水固结法等。换填法是将基础底面以下一定范围内的软弱土层挖去，然后回填强度较高、压缩性较低、并且无侵蚀性的材料，如灰土、砂、碎石等，以提高地基的承载力和减小地基变形。换填法施工简单、成本较低，但处理深度有限，一般适用于浅层软土地基处理。强夯法是利用重锤从高处自由落下，对地基土施加强大的冲击能，使地基土在冲击作用下得以压实、加固，提高地基的承载力和稳定性。强夯法具有加固效果显著、适用范围广、施工速度快等优点，可用于处理砂土、粉土、粘性土、杂填土等多种地基土。但强夯法施工过程中会产生较大的振动和噪声，对周围环境有一定影响，且不适用于饱和软粘土等特殊地基。碎石桩法是通过振动、冲击或水冲等方式在软土地基中成孔，然后将碎石填入孔中，形成碎石桩。碎石桩与桩间土共同作用，形成复合地基，提高地基的承载力和抗变形能力。碎石桩法具有施工工艺简单、造价低、加固效果好等优点，在高填方软基处理中应用广泛。例如，在贵州六盘水机场高填方软基处理工程中，采用振动沉管碎石桩，通过挤密、排水、直接参与承载等作用，有效提高了软弱地基的承载力，较好地解决了边坡稳定和不均匀沉降问题。CFG桩法（水泥粉煤灰碎石桩）是由水泥、粉煤灰、碎石、石屑或砂等混合料加水拌和形成的高粘结强度桩，与桩间土和褥垫层一起形成复合地基。CFG桩法具有承载力提高幅度大、变形小、适用范围广等优点，适用于处理粘性土、粉土、砂土和已自重固结的素填土等地基。排水固结法是通过在地基中设置竖向排水体（如砂井、塑料排水板等）和施加预压荷载，使地基土中的孔隙水排出，土体逐渐固结，强度提高，从而减小地基的沉降和不均匀变形。排水固结法适用于处理饱和软粘土、淤泥质土等透水性较差的地基土。在六盘水机场类似工程中，这些软基处理技术都有不同程度的应用。不同的处理技术在实际应用中表现出不同的效果。例如，换填法对于浅层软土的处理效果明显，能够快速改善地基的承载性能，但对于较深的软土层则难以达到预期效果；强夯法虽然能够有效加固地基，但在周边环境较为敏感的区域，其振动和噪声问题可能会受到限制；碎石桩法在提高地基承载力和增强地基稳定性方面效果显著，但对于桩体的施工质量要求较高，若施工不当可能影响加固效果；CFG桩法在提高地基承载力和控制沉降方面具有优势，但工程造价相对较高；排水固结法对于解决软土地基的沉降问题效果较好，但处理周期较长，需要合理安排施工进度。在实际工程中，需要根据具体的工程地质条件、工程要求、施工条件和经济因素等，综合考虑选择合适的软基处理技术，或者采用多种技术相结合的方式，以达到最佳的处理效果。1.3研究内容与方法1.3.1研究内容本研究聚焦于贵州省六盘水机场东北端高填方体地基稳定性，涵盖多方面关键内容。在高填方体工程性质研究方面，将深入分析填方体的物质组成，精确测定其颗粒级配，明确不同粒径颗粒的分布情况，从而了解其基本构成特性。通过标准试验方法，精准测定填方体的密实度，掌握其紧密程度；细致研究其含水量，明确水分对土体性质的影响；深入探究其孔隙比，了解孔隙结构特征。全面测定抗剪强度参数，包括黏聚力和内摩擦角，为后续的稳定性分析提供关键数据支持。针对高填方体下伏地基的地质条件，本研究将开展详细的地质勘察工作。运用先进的地质勘探技术，全面查明下伏地基的地层结构，确定不同地层的分布、厚度和岩性特征。深入分析地质构造情况，包括断层、褶皱等的分布和特征，评估其对地基稳定性的潜在影响。对地基岩土的物理力学性质进行系统测试，获取密度、压缩性、抗剪强度等参数，为后续的分析提供准确的数据基础。在地基承载力与变形分析方面，本研究将采用多种方法进行综合评估。通过理论计算，运用经典的地基承载力计算公式，结合土体参数，初步估算地基的承载力。开展现场载荷试验，在实际场地中对地基进行加载测试，直接获取地基在不同荷载作用下的变形和承载能力数据。运用数值模拟方法，借助专业的岩土工程分析软件，建立地基的数值模型，模拟不同工况下地基的应力分布和变形情况，预测地基在长期荷载作用下的沉降和变形趋势，评估其是否满足机场工程的要求。本研究还将全面分析影响高填方体地基稳定性的因素。深入探讨土体性质的影响，包括填方体和下伏地基土体的物理力学性质对稳定性的作用机制；研究填方高度的影响，分析不同填方高度下地基的受力状态和稳定性变化规律；分析地基条件的影响，如地基的岩土类型、地质构造等对稳定性的影响；考虑施工工艺的影响，研究不同施工方法和施工顺序对地基稳定性的影响；探讨排水条件的影响，分析地下水的渗流和排泄对地基稳定性的作用。通过敏感性分析，确定各因素对地基稳定性影响的敏感程度，找出影响稳定性的关键因素。基于以上研究内容，本研究将对高填方体地基的稳定性进行全面评价。运用极限平衡法，通过计算滑动力和抗滑力的平衡关系，评估地基在不同工况下的稳定性安全系数；采用数值模拟法，分析地基在各种荷载作用下的应力应变分布情况，判断其是否出现塑性区或破坏区域；结合现场监测数据，验证理论分析和数值模拟的结果，确保评价结果的准确性。根据稳定性评价结果，提出针对性的加固与处理措施。对于稳定性不足的区域，提出合理的地基加固方案，如采用强夯法、碎石桩法、CFG桩法等进行加固处理；优化填方设计，调整填方高度、坡度等参数，提高地基的稳定性；加强排水系统设计，改善地基的排水条件，减少地下水对地基稳定性的不利影响。1.3.2研究方法本研究将综合运用实地勘测、室内试验、数值模拟等多种研究方法，以全面、深入地探究贵州省六盘水机场东北端高填方体地基稳定性。实地勘测是获取第一手资料的重要手段，利用高精度的GPS定位技术，能够精确确定高填方体及周边区域的地理位置坐标，为后续的研究提供准确的空间定位基础。借助全站仪进行地形测量，可详细绘制出场地的地形地貌图，清晰展现地形的起伏变化情况，为分析地基的初始条件提供直观依据。通过地质钻探，从地表向下钻进，获取不同深度的岩土芯样，以便对地基的地层结构、岩土类型和特性进行直观观察和分析。在钻探过程中，记录钻孔的深度、岩性变化、钻进难易程度等信息，为后续的室内试验和分析提供重要的数据支持。进行原位测试，如标准贯入试验、静力触探试验等，能够在现场直接测定岩土的物理力学性质，获取更加真实可靠的数据，为地基稳定性分析提供关键参数。室内试验是研究高填方体和地基岩土物理力学性质的重要环节。通过采取具有代表性的岩土样品，将其带回实验室进行各项测试。进行基本物理性质试验，测定岩土的密度、含水量、孔隙比等参数，了解岩土的基本物理特性。开展颗粒分析试验，运用筛分法、比重计法等手段，确定岩土颗粒的大小分布情况，分析其颗粒级配特征。进行压缩试验，通过对岩土样品施加不同的压力，测定其在压力作用下的压缩变形情况，获取压缩系数、压缩模量等参数，评估岩土的压缩性。进行抗剪强度试验，采用直剪试验、三轴剪切试验等方法，测定岩土的黏聚力和内摩擦角，为地基稳定性分析提供关键的强度参数。数值模拟是本研究的重要方法之一，它能够模拟复杂的工程实际情况，预测地基的变形和稳定性。借助专业的岩土工程分析软件，如Geostudio、Plaxis等，建立高填方体地基的三维数值模型。在模型中，准确输入岩土的物理力学参数，这些参数来源于实地勘测和室内试验的结果，确保模型的真实性和可靠性。设置合理的边界条件，模拟实际工程中的约束和荷载情况，如地基的边界约束、填方的自重荷载、飞机起降产生的动荷载等。通过数值模拟，分析不同工况下地基的应力分布、变形情况以及潜在的滑动面，预测地基在各种条件下的稳定性变化趋势，为工程设计和决策提供科学依据。本研究还将结合工程类比法，收集国内外类似机场高填方工程的成功案例和经验教训，对比六盘水机场东北端高填方体地基的特点和条件，为本次研究提供参考和借鉴。在研究过程中，将多种方法相互验证和补充，确保研究结果的准确性和可靠性。通过实地勘测获取的实际数据，验证室内试验和数值模拟的结果；利用室内试验得到的岩土参数，优化数值模型的输入条件，提高模拟的精度；借助数值模拟的结果，指导实地勘测和室内试验的重点和方向，形成一个有机的研究体系，全面深入地研究六盘水机场东北端高填方体地基稳定性。二、研究区工程地质条件2.1区域地质背景与地震六盘水机场位于滇东高原向黔中高原过渡的区域，地处云贵高原的乌蒙山区，地貌以高原丘陵夹山原峡谷为主，地势起伏较大，垂直落差明显，海拔高度在区域内变化显著。这种复杂的地貌条件使得机场建设面临着诸多挑战，如地形起伏导致的土石方工程量大、填方高度大等问题。从地层单元来看，研究区出露的地层较为复杂，主要包括古生界的寒武系、奥陶系、志留系、泥盆系、石炭系、二叠系，中生界的三叠系以及新生界的第四系。各时代地层岩性差异较大，寒武系主要为一套浅海相碎屑岩和碳酸盐岩沉积，岩石致密坚硬，但受地质构造运动影响，部分区域岩石节理裂隙较为发育；奥陶系以海相沉积的石灰岩、页岩为主，石灰岩质纯，岩溶现象较为常见；志留系多为砂页岩互层，岩石抗风化能力相对较弱；泥盆系和石炭系主要为海陆交互相沉积，岩性包括砂岩、页岩、石灰岩以及煤层等；二叠系则以峨眉山玄武岩大面积喷发和海陆交互相沉积为特征，玄武岩岩性坚硬，而海陆交互相沉积地层岩性复杂，包含多种岩性组合。三叠系主要为陆相碎屑岩沉积，第四系主要为松散堆积物，包括残坡积物、冲洪积物等，其分布厚度和工程性质在不同地段存在较大差异。区域构造上，研究区处于扬子准地台黔北台隆六盘水断陷的中部，经历了多期次的构造运动，地质构造较为复杂。区内主要构造形迹为褶皱和断裂，褶皱轴向主要为北北东-南南西向，褶皱形态多样，有紧闭褶皱、开阔褶皱等，褶皱的存在使得地层产状发生变化，对地基稳定性产生影响。断裂构造主要有北北东向和北西向两组，北北东向断裂规模较大，切割较深，对区域地层的分布和地质构造格局起着重要的控制作用，北西向断裂相对规模较小，但也对局部地质条件产生一定影响。这些断裂构造不仅破坏了地层的完整性，还可能导致岩体破碎、地下水渗流条件改变等问题，增加了地基稳定性的不确定性。在区域地震安全性评价方面，根据国家地震局颁布的《中国地震动参数区划图》（GB18306-2015），研究区的地震动峰值加速度为0.05g，地震基本烈度为Ⅵ度。虽然研究区历史上未发生过强烈地震，但由于地处板块构造活动带的边缘，受到区域构造应力场的影响，仍存在一定的地震活动风险。地震可能引发地基土的液化、震陷等问题，对机场高填方体地基稳定性构成威胁。在机场建设和运营过程中，需要充分考虑地震作用对地基稳定性的影响，采取相应的抗震措施，如加强地基处理、提高建筑物的抗震性能等，以确保机场设施在地震发生时的安全稳定。2.2地形地貌六盘水机场东北端地形呈现出典型的高原山地特征，地势总体西北高、东南低。该区域内山峦起伏，沟谷纵横，地形破碎，地势高差较大。场地原始地面坡度变化较大，一般在15°-45°之间，局部地段坡度可达60°以上。这种较大的坡度使得高填方体在填筑过程中面临着较大的下滑力，对地基稳定性构成潜在威胁。在坡度较陡的区域，填方体与地基之间的摩擦力相对较小，当填方高度增加时，填方体在自重作用下容易产生沿地基坡面的滑动，从而影响地基的整体稳定性。高填方体区域的最大高差可达数十米，这导致了填方厚度的不均匀性。在地势较低的区域，填方厚度较大，而在地势较高的区域，填方厚度相对较小。填方厚度的不均匀会引起地基应力分布的不均匀，进而导致地基产生不均匀沉降。不均匀沉降可能会使机场跑道出现裂缝、变形等问题，严重影响飞机的起降安全。此外，地形的起伏还会导致填方区的排水条件差异较大，在低洼地段容易积水，增加了地基土的含水量，降低了土体的抗剪强度，进一步削弱了地基的稳定性。该区域的地貌类型主要为构造剥蚀低中山地貌和溶蚀侵蚀地貌。构造剥蚀低中山地貌主要由坚硬的岩石组成，岩石风化程度相对较低，山体较为陡峭，地形起伏较大。这种地貌条件下，地基的承载能力相对较高，但由于地形复杂，填方工程难度较大，且填方与地基的结合部位容易出现应力集中现象，对地基稳定性产生不利影响。溶蚀侵蚀地貌则是由于长期的岩溶作用形成，地表和地下存在着大量的岩溶洞穴、溶沟、溶槽等，这些岩溶现象会导致地基的局部强度降低，增加了地基塌陷的风险。在高填方体建设过程中，如果没有对岩溶洞穴等进行有效的处理，填方荷载可能会导致洞穴顶板坍塌，引发地基的突然破坏，严重威胁机场的安全运营。2.3地层岩性研究区内出露的地层主要为第四系全新统人工填土（Q4ml）、残坡积层（Q4el+dl）以及下伏基岩二叠系峨眉山玄武岩（P2β），各层岩土体在岩性、分布及工程性质上存在显著差异。第四系全新统人工填土（Q4ml）主要分布于高填方体区域，是机场建设过程中人工填筑形成的。其物质组成较为复杂，主要由开挖山体的岩石碎块、砂土、粘性土等混合而成。岩石碎块的粒径大小不一，从几厘米到几十厘米不等，含量约占30%-50%；砂土和粘性土则填充于碎块之间，起到胶结和填充的作用。人工填土的结构较为松散，密实度较低，在未经过有效压实处理的区域，其孔隙率较大，一般在30%-40%之间。由于人工填土的物质组成和结构的不均匀性，导致其工程性质差异较大。其压缩性较高，在荷载作用下容易产生较大的变形；抗剪强度较低，尤其是在含水量较高的情况下，土体的抗剪强度会显著降低，对地基稳定性产生不利影响。残坡积层（Q4el+dl）广泛分布于研究区的地表，厚度一般在0.5-3.0m之间，局部地段可达5.0m以上。该层土体主要由岩石风化残积物和坡面堆积物组成，其物质成分与下伏基岩密切相关。在玄武岩分布区域，残坡积层主要为褐红色、灰黑色的粉质粘土，含有少量的玄武岩碎块和角砾，粉质粘土的粘性较强，塑性指数一般在15-20之间。残坡积层的含水量较高，一般在20%-30%之间，天然孔隙比在0.8-1.2之间，属于中压缩性土。其抗剪强度参数相对较低，黏聚力一般在15-30kPa之间，内摩擦角在10°-20°之间。由于残坡积层的厚度和工程性质在不同地段存在较大差异，在地基稳定性分析中，需要充分考虑其不均匀性对地基变形和稳定性的影响。下伏基岩为二叠系峨眉山玄武岩（P2β），岩性致密坚硬，呈灰黑色、黑色。岩石主要由斜长石、辉石等矿物组成，具斑状结构，基质为隐晶质或微晶质。柱状节理发育，节理间距一般在0.5-2.0m之间，节理面多呈闭合状态，少数节理面有少量的次生矿物充填。玄武岩的完整性较好，岩石的抗压强度较高，饱和单轴抗压强度一般在50-100MPa之间，弹性模量在10-30GPa之间。在研究区内，玄武岩的分布较为稳定，埋深在不同地段有所差异，一般在残坡积层以下数米至数十米不等。虽然玄武岩本身的工程性质较好，但在长期的地质作用下，其表层部分岩石可能会发生风化，风化层的厚度一般在1-5m之间，风化后的岩石强度降低，工程性质变差，在地基稳定性分析中需要特别关注。2.4地质构造研究区位于扬子准地台黔北台隆六盘水断陷的中部，经历了多期复杂的构造运动，地质构造特征显著，对高填方体地基稳定性产生着重要影响。区内褶皱构造较为发育，轴向主要呈北北东-南南西向。这些褶皱的形态多样，紧闭褶皱处岩层紧密挤压，岩石破碎程度较高，导致岩体的完整性遭到破坏，力学性能降低。而开阔褶皱虽相对完整，但褶皱转折端等部位仍存在应力集中现象。在高填方体荷载作用下，褶皱构造可能导致地基的不均匀沉降，进而影响高填方体的稳定性。例如，当填方体位于褶皱的翼部时，由于岩层倾斜，填方体与地基之间的接触应力分布不均匀，容易引发填方体沿岩层层面的滑动。区内断裂构造主要有北北东向和北西向两组。北北东向断裂规模较大，切割较深，对区域地层的分布和地质构造格局起着重要的控制作用。断裂带附近岩石破碎，节理裂隙发育，岩体的强度和完整性受到严重破坏。地下水往往沿着断裂带富集和运移，这不仅会降低岩体的抗剪强度，还可能引发地基的渗透变形问题。在高填方体地基中，如果断裂带穿过填方区域，填方荷载可能会促使断裂带进一步活动，导致地基失稳。北西向断裂相对规模较小，但同样会对局部地质条件产生影响，增加地基稳定性的不确定性。断层对高填方体地基稳定性的影响较为复杂。断层的存在破坏了地层的连续性和完整性，使地基岩体的力学性质发生显著变化。断层破碎带内的岩石往往较为破碎，呈碎块状或碎屑状，其抗剪强度低，压缩性高。在高填方体的荷载作用下，断层破碎带容易产生较大的变形，甚至可能引发地基的滑动破坏。此外，断层还可能导致地下水的渗漏和富集，进一步恶化地基的工程地质条件。当高填方体下伏基岩存在断层时，需要特别关注断层的产状、规模以及破碎带的性质，采取有效的工程措施进行处理，以确保地基的稳定性。褶皱对高填方体地基稳定性的影响也不容忽视。褶皱构造使地层产状发生变化，在褶皱的不同部位，地基岩体的受力状态和变形特性存在差异。在褶皱的核部，岩层受挤压强烈，岩石破碎，容易形成软弱带，降低地基的承载能力。而在褶皱的翼部，由于岩层倾斜，填方体与地基之间的摩擦力减小，增加了填方体滑动的风险。例如，当填方高度较大且位于褶皱翼部时，填方体在自重和上部荷载作用下，可能会沿着岩层层面产生滑动，危及地基的稳定性。因此，在分析高填方体地基稳定性时，需要充分考虑褶皱构造的影响，对褶皱的形态、规模和产状进行详细的研究。综上所述，研究区复杂的地质构造对高填方体地基稳定性有着多方面的影响。在机场建设和运营过程中，需要充分认识这些影响，通过详细的地质勘察和分析，准确评估地质构造对地基稳定性的影响程度，采取针对性的工程措施，如地基加固、断层处理、合理设计填方方案等，以确保高填方体地基的稳定，保障机场的安全运营。2.5岩土体物理地质特征2.5.1坡型特征六盘水机场东北端高填方体区域的坡型复杂多样，主要包括直线型、折线型和阶梯型三种类型。直线型边坡多分布于填方区的边缘，其坡度较为均匀，一般在1:1.5-1:2.5之间。这种坡型的边坡在填方高度较小时，稳定性相对较好，但随着填方高度的增加，其下滑力也会相应增大，稳定性逐渐降低。例如，在填方高度为20m时，直线型边坡的安全系数为1.3，基本满足稳定性要求；当填方高度增加到40m时，安全系数降至1.1，接近临界稳定状态。折线型边坡则是由多个不同坡度的直线段组成，通常出现在地形起伏较大的区域。由于不同坡度段的受力情况不同，在折线转折点处容易产生应力集中现象，导致边坡的稳定性降低。例如，某折线型边坡由两段坡度分别为1:1.5和1:2.0的直线段组成，在转折点处，土体的剪应力明显增大，容易引发局部滑坡。此外，折线型边坡的整体稳定性还受到各坡度段长度和高差的影响，当高差较大且坡度较陡时，边坡的稳定性问题更为突出。阶梯型边坡是将边坡分成若干个台阶，每个台阶之间设置一定宽度的平台。这种坡型可以有效地减小边坡的整体下滑力，提高边坡的稳定性。平台的宽度一般在2-5m之间，平台的存在不仅可以增加边坡的抗滑力，还可以起到排水和截流的作用，减少雨水对边坡的冲刷。例如，在某填方高度为60m的阶梯型边坡中，设置了宽度为3m的平台，通过稳定性分析计算，其安全系数达到了1.4，明显高于相同条件下的直线型和折线型边坡。不同坡型对地基稳定性的影响主要体现在以下几个方面：直线型边坡在高填方情况下，下滑力较大，容易导致地基产生滑动破坏；折线型边坡的应力集中问题会使地基局部受力过大，引发地基的不均匀沉降和局部破坏；阶梯型边坡虽然可以提高边坡的稳定性，但在施工过程中，平台的设置会增加施工难度和成本，且如果平台的施工质量控制不当，也会影响地基的稳定性。在进行边坡稳定性分析时，需要根据不同坡型的特点，选择合适的分析方法和参数，准确评估地基的稳定性，为工程设计和施工提供科学依据。2.5.2覆盖层分布特征研究区内覆盖层主要为第四系全新统人工填土（Q4ml）和残坡积层（Q4el+dl），其分布特征对地基稳定性有着重要影响。人工填土主要分布于高填方体区域，是机场建设过程中人工填筑形成的。其厚度在不同位置差异较大，在填方区中心部位，填土厚度可达数十米，而在填方区边缘，填土厚度相对较薄，一般在5-10m之间。人工填土的物质组成复杂，主要由开挖山体的岩石碎块、砂土、粘性土等混合而成，结构松散，密实度较低，孔隙率较大，一般在30%-40%之间。由于其物质组成和结构的不均匀性，导致其工程性质差异较大，压缩性较高，在荷载作用下容易产生较大的变形；抗剪强度较低，尤其是在含水量较高的情况下，土体的抗剪强度会显著降低，对地基稳定性产生不利影响。残坡积层广泛分布于研究区的地表，厚度一般在0.5-3.0m之间，局部地段可达5.0m以上。该层土体主要由岩石风化残积物和坡面堆积物组成，其物质成分与下伏基岩密切相关。在玄武岩分布区域，残坡积层主要为褐红色、灰黑色的粉质粘土，含有少量的玄武岩碎块和角砾。残坡积层的含水量较高，一般在20%-30%之间，天然孔隙比在0.8-1.2之间，属于中压缩性土。其抗剪强度参数相对较低，黏聚力一般在15-30kPa之间，内摩擦角在10°-20°之间。由于残坡积层的厚度和工程性质在不同地段存在较大差异，在地基稳定性分析中，需要充分考虑其不均匀性对地基变形和稳定性的影响。例如，在残坡积层较厚且土体强度较低的区域，地基容易产生较大的沉降和不均匀变形，影响高填方体的稳定性。覆盖层的分布特征对地基稳定性的影响主要体现在承载能力和变形特性两个方面。人工填土由于其结构松散、强度低，其承载能力有限，在高填方荷载作用下，容易产生较大的沉降和变形，导致地基失稳。残坡积层虽然厚度相对较薄，但因其含水量高、压缩性大、抗剪强度低，也会对地基的承载能力和稳定性产生不利影响。此外，覆盖层的不均匀分布会导致地基应力分布不均匀，进而引发地基的不均匀沉降，进一步降低地基的稳定性。在工程建设中，需要对覆盖层进行详细的勘察和分析，采取有效的地基处理措施，如强夯、换填等，提高覆盖层的承载能力和稳定性，确保高填方体地基的安全。2.5.3基岩结构、分布特征研究区下伏基岩为二叠系峨眉山玄武岩（P2β），其结构和分布特征对高填方体地基稳定性起着关键作用。玄武岩岩性致密坚硬，呈灰黑色、黑色，主要由斜长石、辉石等矿物组成，具斑状结构，基质为隐晶质或微晶质。柱状节理发育，节理间距一般在0.5-2.0m之间，节理面多呈闭合状态，少数节理面有少量的次生矿物充填。岩石的完整性较好，饱和单轴抗压强度一般在50-100MPa之间，弹性模量在10-30GPa之间。在研究区内，玄武岩的分布较为稳定，埋深在不同地段有所差异，一般在残坡积层以下数米至数十米不等。然而，在长期的地质作用下，其表层部分岩石可能会发生风化，风化层的厚度一般在1-5m之间。风化后的岩石强度降低，工程性质变差，其饱和单轴抗压强度可降至20-50MPa，弹性模量也会相应减小。风化层的存在对高填方体地基稳定性产生不利影响，在填方荷载作用下，风化层容易产生变形和破坏，进而影响下部完整基岩的承载能力和稳定性。基岩的结构和分布特征对高填方体地基稳定性的影响主要体现在岩石强度和完整性方面。致密坚硬的玄武岩具有较高的强度和较好的承载能力，能够为高填方体提供稳定的基础。然而，柱状节理的发育会在一定程度上削弱岩石的完整性，增加岩石的渗透性，使得地下水更容易沿着节理面渗透，降低岩石的抗剪强度。当填方荷载较大时，节理面可能会成为潜在的滑动面，导致地基失稳。此外，风化层的存在进一步降低了基岩的强度和稳定性，在地基设计和施工中，需要对风化层进行妥善处理，如清除或加固，以确保地基的安全。例如，在某工程中，通过对风化层进行强夯加固处理，提高了风化层的强度和密实度，有效增强了地基的稳定性。2.6水文地质条件研究区的水文地质条件对六盘水机场东北端高填方体地基稳定性有着不容忽视的影响。区内地下水类型主要包括孔隙水、裂隙水和岩溶水，它们在赋存状态、补给排泄条件以及对地基稳定性的作用机制上各有特点。孔隙水主要赋存于第四系全新统人工填土（Q4ml）和残坡积层（Q4el+dl）中。人工填土结构松散，孔隙率较大，一般在30%-40%之间，具有一定的透水性，能够储存一定量的孔隙水。残坡积层的孔隙水则主要存在于土体的孔隙和裂隙中，其含水量较高，一般在20%-30%之间，天然孔隙比在0.8-1.2之间。孔隙水的补给来源主要为大气降水，降水通过地表入渗进入土体孔隙中。在地形低洼处，孔隙水容易汇聚，形成局部的滞水区域。其排泄方式主要为蒸发和向深部地层的渗流。孔隙水的存在会增加土体的重量，降低土体的抗剪强度，尤其是在高填方体的作用下，孔隙水压力升高，进一步削弱了土体的稳定性。例如，当人工填土中的孔隙水含量增加时，填方体的自重增大，同时土体的抗剪强度降低，容易导致填方体产生滑动破坏。裂隙水主要分布于下伏基岩二叠系峨眉山玄武岩（P2β）的节理裂隙中。玄武岩柱状节理发育，节理间距一般在0.5-2.0m之间，这些节理裂隙为裂隙水的赋存和运移提供了通道。裂隙水的补给来源主要为大气降水和孔隙水的下渗。在地形较高处，裂隙水接受降水补给后，沿节理裂隙向低处渗流。其排泄方式主要为在地形低洼处出露地表，形成泉水或汇入地表水体。裂隙水的存在会对基岩的力学性质产生影响，降低岩石的抗剪强度，尤其是在节理面被水浸泡后，节理面的摩擦力减小，容易引发基岩的滑动破坏。此外，裂隙水的渗流还可能导致地基的渗透变形，如管涌、流土等现象的发生。岩溶水主要存在于研究区内的岩溶洞穴、溶沟、溶槽等岩溶通道中。由于研究区存在一定的岩溶发育，岩溶水的分布较为复杂。岩溶水的补给主要来自大气降水和地表水的入渗，其排泄方式主要为通过岩溶管道排入附近的河流或地下暗河。岩溶水的水位和水量变化较大，受降水和季节影响明显。在雨季，岩溶水水位上升，水量增大；在旱季，水位下降，水量减少。岩溶水对地基稳定性的影响主要体现在可能导致地基塌陷和岩溶渗漏问题。当高填方体荷载作用于岩溶地区时，如果岩溶洞穴顶板的承载能力不足，可能会发生塌陷，导致地基失稳。此外，岩溶水的渗漏还可能引起地基土的流失，降低地基的承载能力。地下水水位在研究区内呈现出一定的变化规律。在地势较高的区域，地下水水位相对较深，一般在地面以下5-10m；在地势低洼处，地下水水位相对较浅，可能接近地表。地下水水位的变化受降水、蒸发、地表水体补给等多种因素的影响。在雨季，降水充沛，地下水水位上升；在旱季，蒸发量大，地下水水位下降。此外，地表水体如河流、湖泊等对地下水水位也有一定的影响，当河流丰水期时，河水补给地下水，水位上升；当河流枯水期时，地下水补给河水，水位下降。研究区地下水的补给主要来自大气降水，降水通过地表入渗进入地下含水层。此外，地表水体如河流、池塘等也会对地下水进行补给。排泄方式主要包括蒸发、向地表水体排泄以及人工开采等。在天然状态下，地下水的补给和排泄基本处于平衡状态，但在机场建设和运营过程中，工程活动可能会改变地下水的补给和排泄条件。例如，填方工程可能会堵塞部分地下水的渗流通道，导致地下水水位升高；排水工程的建设则可能加速地下水的排泄，降低地下水水位。水文地质条件对地基稳定性的影响是多方面的。地下水的存在会改变土体和岩石的物理力学性质，降低其抗剪强度和承载能力。孔隙水和裂隙水的压力变化会影响地基的应力状态，增加地基失稳的风险。岩溶水导致的地基塌陷和渗漏问题更是直接威胁到地基的稳定性。在高填方体地基稳定性分析和工程设计中，必须充分考虑水文地质条件的影响，采取有效的排水措施，降低地下水水位，减少地下水对地基的不利影响。例如，设置排水盲沟、排水管网等，及时排除地基中的孔隙水和裂隙水，防止地下水积聚对地基稳定性造成破坏。2.7主要工程地质问题六盘水机场东北端高填方体地基存在多种主要工程地质问题，这些问题对机场的安全运营构成潜在威胁，需予以高度重视。地基沉降是一个关键问题，研究区内第四系全新统人工填土结构松散、密实度低、孔隙率大，一般在30%-40%之间，在高填方荷载作用下，土体颗粒重新排列，孔隙被压缩，容易产生较大的沉降。下伏基岩若存在风化层，其强度降低，工程性质变差，也会在填方荷载下产生变形，进一步加剧地基沉降。地基沉降可能导致机场跑道出现裂缝、变形，影响飞机起降安全，增加飞行事故的风险。不均匀沉降还可能使机场建筑物出现倾斜、开裂等问题，危及建筑物的结构安全。边坡失稳也是一个突出问题。高填方体区域地形起伏大，边坡坡度较陡，部分地段坡度可达60°以上，在填方自重和外部荷载作用下，边坡土体容易沿潜在滑动面产生滑动破坏。区内的地质构造如褶皱和断裂，破坏了土体和岩体的完整性，降低了其抗剪强度，增加了边坡失稳的可能性。降雨入渗会使边坡土体含水量增加，重度增大，抗剪强度降低，孔隙水压力升高，进一步削弱边坡的稳定性，引发滑坡、坍塌等地质灾害。岩溶问题在研究区也较为突出。溶蚀侵蚀地貌导致区内存在大量的岩溶洞穴、溶沟、溶槽等，岩溶洞穴顶板在填方荷载作用下可能发生坍塌，导致地基突然下沉，影响机场设施的正常使用。岩溶水的渗漏会带走地基土颗粒，导致地基土流失，强度降低，增加地基变形和失稳的风险。地下水问题同样不容忽视。孔隙水和裂隙水的存在会增加土体和岩体的重量，降低其抗剪强度。在高填方体地基中，地下水水位的变化会引起孔隙水压力的改变，当孔隙水压力升高时，有效应力减小，土体抗剪强度降低，容易引发地基失稳。岩溶水的水位和水量变化较大，对地基稳定性产生不利影响，可能导致地基塌陷和岩溶渗漏等问题。此外，地基承载力不足也是一个潜在问题。人工填土和残坡积层的工程性质较差，承载能力有限，难以满足高填方体和机场设施的荷载要求。若下伏基岩存在软弱夹层或风化层，也会降低地基的整体承载能力，导致地基在荷载作用下产生过大的变形和破坏。综上所述，六盘水机场东北端高填方体地基存在的主要工程地质问题相互关联、相互影响，对机场的安全运营构成严重威胁。在机场建设和运营过程中，必须充分认识这些问题，通过详细的地质勘察、科学的分析计算和合理的工程措施，加以解决和防范，确保机场地基的稳定和安全。三、高填方体地基稳定性研究3.1高填方体地基稳定性的地质评价根据详细的地质勘察资料，六盘水机场东北端高填方体地基稳定性受多种地质因素综合影响，需从多个方面进行定性评价。从地形地貌角度来看，研究区地势起伏大，西北高、东南低，原始地面坡度多在15°-45°，局部达60°以上，高填方体最大高差可达数十米。较大的坡度使填方体在自重作用下产生较大下滑力，填方厚度的不均匀导致地基应力分布不均，容易引发不均匀沉降，对地基稳定性极为不利。如在坡度较陡的区域，填方体与地基间摩擦力小，填方高度增加时，填方体易沿地基坡面滑动；填方厚度差异大的区域，地基沉降差异明显，可能导致跑道等设施变形。地层岩性方面，第四系全新统人工填土结构松散、密实度低、孔隙率大，一般在30%-40%之间，压缩性高、抗剪强度低，在高填方荷载下易产生较大变形和滑动。残坡积层含水量高，一般在20%-30%之间，天然孔隙比在0.8-1.2之间，属于中压缩性土，抗剪强度参数低，黏聚力一般在15-30kPa之间，内摩擦角在10°-20°之间，其不均匀分布会影响地基稳定性。下伏基岩峨眉山玄武岩虽岩性致密坚硬，但表层存在1-5m厚的风化层，风化后岩石强度降低，饱和单轴抗压强度可降至20-50MPa，弹性模量减小，降低了地基的承载能力和稳定性。地质构造对地基稳定性影响显著。区内褶皱轴向主要为北北东-南南西向，褶皱使地层产状改变，在褶皱核部和翼部，地基岩体受力状态和变形特性不同，易引发地基不均匀沉降和滑动。如在褶皱核部，岩层受挤压破碎，地基承载能力降低；在翼部，填方体与地基摩擦力减小，增加滑动风险。断裂构造主要有北北东向和北西向两组，断裂带附近岩石破碎、节理裂隙发育，岩体强度和完整性破坏，地下水沿断裂带富集运移，降低岩体抗剪强度，引发地基渗透变形，填方荷载可能促使断裂带活动，导致地基失稳。岩土体物理地质特征也不容忽视。坡型方面，直线型边坡在高填方时下滑力大，稳定性随填方高度增加而降低；折线型边坡在转折点易产生应力集中，稳定性受各坡度段长度和高差影响；阶梯型边坡设置平台可提高稳定性，但施工质量控制不当会影响地基稳定性。覆盖层中人工填土厚度和工程性质差异大，残坡积层厚度和性质不均匀，都对地基稳定性产生不利影响。基岩的柱状节理发育虽增强了岩石渗透性，但也降低了其完整性，风化层进一步削弱基岩强度和稳定性。水文地质条件同样影响地基稳定性。孔隙水、裂隙水和岩溶水在研究区广泛分布，孔隙水增加土体重量、降低抗剪强度；裂隙水降低岩石抗剪强度，引发渗透变形；岩溶水可能导致地基塌陷和渗漏。地下水水位受降水、蒸发、地表水体补给等因素影响而变化，工程活动也会改变其补给和排泄条件，进而影响地基稳定性。综上所述，六盘水机场东北端高填方体地基地质条件复杂，多种地质因素相互作用，对地基稳定性产生不利影响。在机场建设和运营中，需充分认识这些因素，采取有效措施保障地基稳定，确保机场安全运营。3.2极限平衡法稳定性计算3.2.1计算方法极限平衡法是边坡稳定性分析中应用广泛的经典方法，其基本原理是基于刚体极限平衡理论，将滑动土体视为由多个条块组成的刚体系统。假定滑动面已经确定，通过分析条块上的作用力，建立力和力矩的平衡方程，从而求解边坡的稳定安全系数。安全系数被定义为滑面上的抗滑力与滑动力之比，当安全系数大于1时，表明边坡处于稳定状态；等于1时，边坡处于极限平衡状态；小于1时，则边坡不稳定。在众多极限平衡法中，瑞典圆弧法是最早提出的条分法，由Petterson于1916年提出。该方法假设边坡的滑裂面为圆柱面，将滑动土体沿滑裂面分成若干垂直土条，在分析时不考虑土条之间的相互作用力。计算时，对每个土条进行受力分析，将土条的重力分解为切向力和法向力，切向力为滑动力，法向力通过滑面的抗剪强度提供抗滑力。瑞典圆弧法计算简单，但其未考虑土条间相互作用，导致计算结果偏于保守，一般适用于均质土坡且对计算精度要求不高的情况。Bishop法是对瑞典圆弧法的改进，由Bishop于1955年提出。该方法考虑了土条间的相互作用力，假设土条两侧的条间力合力作用点位置，并假定条间力的合力方向为水平。在计算时，通过对土条进行力的平衡分析，建立平衡方程求解安全系数。Bishop法在一定程度上克服了瑞典圆弧法的不足，计算结果相对更接近实际情况，是工程中常用的边坡稳定性分析方法之一，适用于各种土质边坡，尤其对于非均质土坡的分析更为准确。在实际应用中，应根据六盘水机场东北端高填方体地基的具体情况选择合适的计算方法。考虑到该区域地基土体性质的复杂性和不均匀性，Bishop法由于其对土条间相互作用的合理考虑，能够更准确地反映地基的实际受力状态，因此在本研究中具有较高的适用性。同时，为了确保计算结果的可靠性，可结合多种方法进行对比分析，综合评估地基的稳定性。3.2.2计算剖面、工况与参数取值为准确分析六盘水机场东北端高填方体地基稳定性，需合理确定计算剖面、工况及参数取值。根据场地地形地貌、地层岩性、地质构造等条件，选取具有代表性的A-A’、B-B’、C-C’三个剖面进行计算。A-A’剖面位于填方区中部，填方高度较大，且下伏基岩起伏较大，能反映填方区中部的典型地质特征；B-B’剖面靠近填方区边缘，地形坡度较陡，可考察边坡边缘的稳定性情况；C-C’剖面穿越了一条小型断裂构造，能研究地质构造对地基稳定性的影响。针对不同工况，考虑以下荷载组合：工况一为自重工况，仅考虑高填方体和地基土体的自重作用；工况二为自重+地下水工况，在自重基础上，考虑地下水对土体的浮力和渗透压力作用；工况三为自重+地震工况，考虑高填方体和地基土体自重以及地震作用产生的惯性力。在自重工况下，土体的重度根据室内试验结果取值，人工填土重度取18kN/m³，残坡积层重度取19kN/m³，下伏基岩重度取26kN/m³。在自重+地下水工况中，根据水文地质勘察资料，确定地下水位线位置，计算地下水浮力和渗透压力。浮力根据阿基米德原理计算，渗透压力通过渗流计算确定，采用达西定律计算渗流速度，进而得到渗透压力。在自重+地震工况下，地震作用采用拟静力法考虑，根据区域地震安全性评价结果，研究区地震基本烈度为Ⅵ度，地震动峰值加速度为0.05g，地震作用系数取0.05，通过在条块上施加水平和竖向地震惯性力来考虑地震作用。计算参数的取值直接影响计算结果的准确性，因此需合理选取。抗剪强度参数根据室内试验结果和经验取值，人工填土黏聚力取15kPa，内摩擦角取20°；残坡积层黏聚力取20kPa，内摩擦角取15°；下伏基岩黏聚力取500kPa，内摩擦角取35°。这些参数取值综合考虑了土体的物质组成、结构特征以及现场实际情况，并参考了类似工程的经验数据。通过合理确定计算剖面、工况和参数取值，能够更全面、准确地分析高填方体地基在不同条件下的稳定性，为工程设计和决策提供可靠依据。3.2.3安全系数根据极限平衡法的计算原理，对于采用Bishop法计算的高填方体地基，其安全系数计算公式为：F_s=\frac{\sum_{i=1}^{n}\frac{1}{m_{\thetai}}\left(c_il_i+\left[W_i+\left(P_{i+1}-P_i\right)\sin\alpha_i\right]\cos\alpha_i\tan\varphi_i\right)}{\sum_{i=1}^{n}\left(W_i+\left(P_{i+1}-P_i\right)\sin\alpha_i\right)\sin\alpha_i}其中，F_s为安全系数；n为土条数量；c_i为第i个土条滑面上的黏聚力；l_i为第i个土条滑面的长度；W_i为第i个土条的重力；P_i、P_{i+1}分别为第i个土条两侧条间力的合力；\alpha_i为第i个土条滑面与水平面的夹角；\varphi_i为第i个土条滑面上的内摩擦角；m_{\thetai}=\cos\alpha_i+\frac{\sin\alpha_i\tan\varphi_i}{F_s}。通过对各计算剖面在不同工况下进行计算，得到的安全系数结果如下表所示：计算剖面工况一（自重）工况二（自重+地下水）工况三（自重+地震）A-A’1.351.201.05B-B’1.281.150.98C-C’1.151.020.85根据《建筑边坡工程技术规范》（GB50330-2013），对于机场高填方体地基，其稳定安全系数的最小值应满足：一级边坡安全系数不小于1.35，二级边坡安全系数不小于1.30，三级边坡安全系数不小于1.25。对比规范标准，A-A’剖面在工况一和工况二下安全系数满足要求，但在工况三（自重+地震）下安全系数为1.05，小于规范要求的最小值，表明在地震作用下该剖面存在一定的不稳定风险；B-B’剖面在工况一下安全系数接近规范要求，工况二和工况三下安全系数均小于规范要求，说明该剖面在正常工况下稳定性尚可，但在地下水和地震作用下稳定性较差；C-C’剖面在三种工况下安全系数均小于规范要求，尤其是在工况三（自重+地震）下安全系数仅为0.85，稳定性严重不足，需进行加固处理。通过对安全系数的计算和对比分析，能够直观地评估高填方体地基在不同工况下的稳定性状况，为后续的稳定性分析和工程措施制定提供重要依据。3.2.4计算结果与稳定性分析从极限平衡法的计算结果来看，各计算剖面在不同工况下的稳定性存在明显差异。A-A’剖面在自重工况下安全系数为1.35，处于相对稳定状态，这主要得益于该剖面位于填方区中部，填方高度虽较大，但地基土体分布相对均匀，且下伏基岩的承载能力较强，能够较好地支撑填方体的重量。然而，在自重+地下水工况下，安全系数降至1.20，这是因为地下水的存在增加了土体的重量，降低了土体的抗剪强度，同时地下水的渗透压力也对地基稳定性产生了不利影响。在自重+地震工况下，安全系数进一步降低至1.05，地震作用产生的惯性力使得地基土体所受的外力增大，超过了土体的抗滑能力，导致稳定性下降。B-B’剖面在自重工况下安全系数为1.28，接近规范要求，这是由于该剖面靠近填方区边缘，地形坡度较陡，填方高度相对较小，地基土体的受力相对较小。但在自重+地下水工况下，安全系数降至1.15，在自重+地震工况下，安全系数降至0.98，表明该剖面在地下水和地震作用下稳定性较差。这是因为地形坡度较陡使得土体在自重作用下就有较大的下滑趋势，地下水和地震作用进一步加剧了这种趋势，导致土体的稳定性降低。C-C’剖面在三种工况下安全系数均小于规范要求，这是因为该剖面穿越了一条小型断裂构造，断裂构造破坏了地基土体的完整性，降低了土体的强度和抗滑能力。在自重工况下，安全系数为1.15，这是由于断裂构造附近土体破碎，承载能力下降，无法有效支撑填方体的重量。在自重+地下水工况下，安全系数降至1.02，地下水沿着断裂构造渗透，进一步软化了土体，降低了土体的抗剪强度。在自重+地震工况下，安全系数降至0.85，地震作用使得断裂构造附近的土体更容易发生滑动和变形，导致稳定性严重不足。综合分析可知，影响高填方体地基稳定性的因素主要包括土体强度、坡角、荷载等。土体强度是影响地基稳定性的关键因素，土体的黏聚力和内摩擦角越大，其抗剪强度越高，地基的稳定性越好。如人工填土黏聚力和内摩擦角相对较低，其稳定性较差；而下伏基岩黏聚力和内摩擦角较高，稳定性较好。坡角越大，土体在自重作用下的下滑力越大，地基的稳定性越差。如B-B’剖面地形坡度较陡，其稳定性明显低于A-A’剖面。荷载的大小和类型也对地基稳定性有重要影响，除了自重荷载外，地下水和地震等荷载会增加土体的重量或产生附加力，降低土体的抗滑能力，从而影响地基的稳定性。在工程设计和施工中，应充分考虑这些因素，采取相应的措施提高地基的稳定性，如对地基进行加固处理、优化填方设计、加强排水系统建设等，以确保机场高填方体地基的安全稳定。3.3二维有限元计算3.3.1有限单元法计算的安全系数有限单元法（FEM）作为一种强大的数值分析方法，在高填方体地基稳定性分析中发挥着关键作用。其基本原理是将连续的求解域离散为有限个单元的组合体，通过对每个单元进行力学分析，建立单元的平衡方程，再将这些单元方程组装成整个求解域的总体平衡方程，从而求解出节点的位移和应力。在高填方体地基稳定性分析中，有限单元法具有显著的优势。它能够精确地考虑土体的非线性特性，土体在实际受力过程中，其应力-应变关系往往呈现出非线性特征，有限单元法通过采用合适的本构模型，如摩尔-库仑本构模型、邓肯-张本构模型等，能够准确地模拟土体的这种非线性行为。有限单元法还可以全面地考虑边界条件的影响。在高填方体地基中，边界条件复杂多样，包括地基与基础的接触条件、地下水的渗流边界条件等。有限单元法能够根据实际情况，准确地设定这些边界条件，使计算结果更加符合实际工程情况。此外，有限单元法可以直观地给出地基的应力应变分布情况，通过计算结果的可视化处理，能够清晰地展示地基在不同荷载作用下的应力集中区域、变形较大区域等，为工程人员分析地基稳定性提供直观的依据。在有限单元法中，安全系数的计算通常采用强度折减法。强度折减法的基本思想是通过逐步折减土体的抗剪强度参数（黏聚力c和内摩擦角\varphi），直至地基达到极限平衡状态，此时的折减系数即为安全系数。具体计算过程如下：首先，建立高填方体地基的有限元模型，赋予模型中土体初始的抗剪强度参数c_0和\varphi_0。然后，按照一定的步长逐步折减抗剪强度参数，折减后的参数为c=c_0/F，\varphi=\varphi_0/F，其中F为折减系数。在每一次折减后，进行有限元计算，求解地基的应力应变分布。当计算结果表明地基达到极限平衡状态，如出现塑性区贯通、位移急剧增大等现象时，此时的折减系数F即为地基的安全系数。强度折减法通过这种迭代计算的方式，能够准确地确定地基的安全系数，为高填方体地基稳定性评价提供了一种可靠的方法。3.3.2有限元稳定安全系数计算与结果分析利用专业的有限元软件，如Geostudio中的SLOPE/W模块，对六盘水机场东北端高填方体地基进行二维有限元计算。在建立有限元模型时，充分考虑研究区的地形地貌、地层岩性、地质构造以及水文地质条件等因素。根据实际测量数据，准确绘制高填方体的几何形状和边界条件，确保模型的真实性。将地层划分为人工填土、残坡积层和下伏基岩等不同单元，为每个单元赋予相应的物理力学参数，这些参数来源于室内试验和现场测试结果。在计算过程中，采用摩尔-库仑本构模型来描述土体的力学行为。该模型基于摩尔-库仑强度理论，能够较好地反映土体在剪切破坏时的特性。设置合适的边界条件，底部边界固定，限制水平和竖向位移；侧面边界根据实际情况，限制水平位移。考虑不同工况，包括自重工况、自重+地下水工况以及自重+地震工况，模拟地基在不同条件下的受力状态。通过有限元计算，得到不同工况下地基的应力应变分布和塑性区开展情况。在自重工况下，地基的竖向应力随着填方高度的增加而增大，在填方体底部和下伏基岩顶部出现应力集中现象。水平应力则在填方区边缘较大，向内部逐渐减小。从应变分布来看，填方体和残坡积层的变形较大，下伏基岩的变形相对较小。塑性区主要分布在填方体与残坡积层的界面处以及填方区边缘，这表明这些区域是地基的薄弱部位，容易发生破坏。在自重+地下水工况下，由于地下水的存在，地基土体的有效应力减小，抗剪强度降低。此时，地基的竖向应力分布与自重工况相似，但水平应力有所增加，这是由于地下水的渗透压力作用。应变分布也发生了变化，填方体和残坡积层的变形进一步增大，尤其是在地下水位附近，变形更为明显。塑性区的范围扩大，不仅在填方体与残坡积层的界面处和填方区边缘，还向残坡积层内部扩展，说明地下水对地基稳定性的影响较为显著。在自重+地震工况下，地震作用产生的惯性力使地基的应力应变状态更加复杂。竖向应力和水平应力在地震波的作用下呈现出动态变化，在地震力的作用方向上，应力明显增大。应变分布也出现了较大的波动，填方体和基岩的变形都有所增加，尤其是在填方区的上部和下部，变形更为突出。塑性区进一步扩展，可能出现贯通的滑动面，导致地基失稳。根据强度折减法，计算得到不同工况下地基的稳定安全系数。在自重工况下，安全系数为1.38，表明地基在自重作用下处于稳定状态，但安全储备相对较小。在自重+地下水工况下，安全系数降至1.22，说明地下水的存在对地基稳定性产生了较大影响，降低了地基的安全储备。在自重+地震工况下，安全系数为1.08，接近临界稳定状态，地震作用显著降低了地基的稳定性，存在较大的安全隐患。综合分析有限元计算结果，可知地下水和地震是影响六盘水机场东北端高填方体地基稳定性的重要因素。地下水的存在降低了土体的有效应力和抗剪强度，增大了地基的变形和塑性区范围，从而降低了地基的稳定性。地震作用产生的惯性力使地基的应力应变状态恶化，增加了地基失稳的风险。在工程设计和施工中，应充分考虑这些因素，采取有效的措施，如加强排水系统建设、提高地基的抗震性能等，以提高地基的稳定性，确保机场的安全运营。3.4稳定性综合评价综合极限平衡法和有限元法的计算结果，并结合地质评价，对六盘水机场东北端高填方体地基稳定性可进行全面、深入的评价。从地质评价来看，研究区复杂的地质条件对地基稳定性产生了诸多不利影响。地形地貌上，地势起伏大，坡度陡峭，填方厚度不均，易引发地基的不均匀沉降和滑动。地层岩性方面，人工填土结构松散、强度低，残坡积层含水量高、压缩性大，下伏基岩虽整体强度较高，但表层风化层强度降低，这些都削弱了地基的稳定性。地质构造上，褶皱和断裂破坏了土体和岩体的完整性，降低了其抗剪强度，增加了地基失稳的风险。岩土体物理地质特征中，坡型的复杂性、覆盖层的不均匀分布以及基岩节理裂隙的发育和风化层的存在，都对地基稳定性产生了负面影响。水文地质条件方面，孔隙水、裂隙水和岩溶水的存在增加了土体和岩体的重量，降低了抗剪强度，地下水水位的变化也对地基稳定性产生影响。极限平衡法计算结果表明，不同计算剖面在不同工况下的稳定性存在差异。A-A’剖面在自重工况下安全系数为1.35，相对稳定，但在自重+地下水工况下安全系数降至1.20，在自重+地震工况下安全系数仅为1.05，低于规范要求，存在不稳定风险。B-B’剖面在自重工况下安全系数为1.28，接近规范要求，但在自重+地下水和自重+地震工况下安全系数均小于规范要求，稳定性较差。C-C’剖面由于穿越断裂构造，在三种工况下安全系数均小于规范要求，尤其是在自重+地震工况下安全系数仅为0.85，稳定性严重不足。有限元法计算结果显示，在自重工况下，地基的安全系数为1.38，处于稳定状态，但安全储备较小。在自重+地下水工况下，安全系数降至1.22，说明地下水对地基稳定性影响较大。在自重+地震工况下，安全系数为1.08，接近临界稳定状态，地震作用显著降低了地基的稳定性。对比极限平衡法和有限元法的计算结果，两者在趋势上具有一致性，都表明地下水和地震是影响地基稳定性的重要因素。有限元法能够更全面地考虑土体的非线性特性和边界条件，其计算结果相对更接近实际情况，但计算过程较为复杂，需要较高的计算资源。极限平衡法计算简单，能够快速给出安全系数，但在考虑土体的复杂特性和边界条件方面存在一定局限性。综合来看，六盘水机场东北端高填方体地基在当前条件下存在一定的稳定性问题，尤其是在地下水和地震作用下，稳定性明显降低。为确保机场的安全运营，需要采取有效的加固与处理措施。对于稳定性不足的区域，可采用强夯法对地基进行加固，通过强大的冲击能使地基土密实，提高地基的承载能力和稳定性。采用碎石桩法或CFG桩法，形成复合地基，增强地基的强度和抗变形能力。优化填方设计，调整填方高度和坡度，减少填方荷载对地基的不利影响。加强排水系统设计，降低地下水水位，减少地下水对地基稳定性的影响。在机场建设和运营过程中，应加强对地基稳定性的监测，及时发现问题并采取相应的措施，以保障机场的安全稳定运行。四、变形破坏的二维有限元分析4.1有限单元法基本理论4.1.1有限单元法概述有限单元法（FiniteElementMethod，FEM）是一种用于求解数学物理问题的数值计算方法，其基本概念是将连续的求解域离散为有限个单元的组合体，通过对每个单元进行力学分析，建立单元的平衡方程，再将这些单元方程组装成整个求解域的总体平衡方程，从而求解出节点的位移和应力。有限单元法的发展历程充满了创新与突破。其起源可以追溯到20世纪40年代，Courant在1943年取定义在三角形域上分片连续函数，利用最小势能原理研究St.Venant的扭转问题，这一开创性的工作为有限单元法的发展奠定了基础。在20世纪50年代，有限单元法在固体力学领域得到了重要发展，1956年Turner、Clough把刚架位移法（直接刚度法）应用到弹性力学平面应力问题中，他们将结构划分成一个个三角形和矩形的“单元”，通过单元刚度矩阵来表示每一单元的特性，这是有限单元法的重要雏形。1960年，Clough在论文《平面应力分析的有限单元法》中首次提出了“有限单元法”这个名称，标志着有限单元法正式诞生。此后，众多数学家和工程师的加盟使得有限单元法进入了黄金发展阶段，其理论不断完善，应用范围也不断扩大。有限单元法的应用范围极为广泛，涵盖了多个学科领域。在航空航天领域，有限单元法被广泛应用于飞机结构的静、动态特性分析，通过对飞机结构进行有限元建模和分析，可以准确预测飞机在各种工况下的应力、应变和位移，为飞机的设计和优化提供重要依据。在机械工程领域，有限单元法可用于机械零部件的强度和刚度分析，例如对发动机零部件、汽车底盘等进行分析，确保其在工作过程中的可靠性和安全性。在土木工程领域，有限单元法在岩土工程、结构工程等方面发挥着重要作用。在岩土工程中，可用于分析地基沉降、边坡稳定性、地下洞室围岩稳定性等问题；在结构工程中，可用于分析建筑物、桥梁、大坝等结构的受力性能和抗震性能。在生物医学领域，有限单元法可用于模拟人体骨骼、器官等的力学行为，为医学研究和临床治疗提供支持。在热传导、电磁场、流体力学等领域，有限单元法也有着广泛的应用，能够有效解决这些领域中的复杂问题。在岩土工程中，有限单元法具有不可替代的重要作用。它能够考虑土体的非线性特性，土体在实际受力过程中，其应力-应变关系往往呈现出非线性特征，有限单元法通过采用合适的本构模型，如摩尔-库仑本构模型、邓肯-张本构模型等，能够准确地模拟土体的这种非线性行为。有限单元法可以全面考虑边界条件的影响，在岩土工程中，边界条件复杂多样，包括地基与基础的接触条件、地下水的渗流边界条件等，有限单元法能够根据实际情况，准确地设定这些边界条件，使计算结果更加符合实际工程情况。此外，有限单元法还可以直观地给出地基的应力应变分布情况，通过计算结果的可视化处理，能够清晰地展示地基在不同荷载作用下的应力集中区域、变形较大区域等，为工程人员分析地基稳定性提供直观的依据。总之，有限单元法为岩土工程的分析和设计提供了一种强大的工具，能够有效解决岩土工程中的复杂问题，提高工程的安全性和可靠性。4.1.2有限元基本方程有限元法的基本方程推导基于弹性力学的基本原理，包括平衡方程、几何方程和物理方程。以平面问题为例，假设物体在笛卡尔坐标系下，其平衡方程表示为：\begin{cases}\frac{\partial\sigma_{x}}{\partialx}+\frac{\partial\tau_{xy}}{\partialy}+f_{x}=0\\\frac{\partial\tau_{xy}}{\partialx}+\frac{\partial\sigma_{y}}{\partialy}+f_{y}=0\end{cases}其中，\sigma_{x}、\sigma_{y}分别为x、y方向的正应力，\tau_{xy}为剪应力，f_{x}、f_{y}分别为x、y方向的体力。几何方程描述了物体的应变与位移之间的关系，对于平面问题，其几何方程为：\begin{cases}\varepsilon_{x}=\frac{\partialu}{\partialx}\\\varepsilon_{y}=\frac{\partialv}{\partialy}\\\gamma_{xy}=\frac{\partialu}{\partialy}+\frac{\partialv}{\partialx}\end{cases}其中，\varepsilon_{x}、\varepsilon_{y}分别为x、y方向的正应变，\gamma_{xy}为剪应变，u、v分别为x、y方向的位移。物理方程则建立了应力与应变之间的关系，对于各向同性弹性体，其物理方程为：\begin{cases}\sigma_{x}=\frac{E}{1-\mu^{2}}(\varepsilon_{x}+\mu\varepsilon_{y})\\\sigma_{y}=\frac{E}{1-\mu^{2}}(\varepsilon_{y}+\mu\varepsilon_{x})\\\tau_{xy}=\frac{E}{2(1+\mu)}\gamma_{xy}\end{cases}其中，E为弹性模量，\mu为泊松比。在有限元分析中，将求解域离散为有限个单元，每个单元内的位移可以通过节点位移进行插值表示。