基于GIS的晋江市工程地质环境功能区划：模型构建与实践应用

基于GIS的晋江市工程地质环境功能区划：模型构建与实践应用一、引言1.1研究背景与意义晋江市作为福建省东南沿海的重要城市，地理位置优越，经济发展迅速，在区域发展中占据重要地位。然而，随着城市化进程的不断加快，城市建设规模日益扩大，对土地资源的需求持续增加，工程建设活动愈发频繁，这使得晋江市的工程地质环境面临着诸多挑战。工程地质环境是城市建设和发展的基础条件，其质量的优劣直接影响着城市工程建设的安全性、稳定性和可持续性。晋江市的地形地貌复杂多样，北部沿海平原广阔，南部山区连绵起伏，地层岩性差异较大，地质构造活动频繁，这些自然因素使得工程地质条件较为复杂。同时，人类工程活动如大规模的土地开发、高层建筑建设、地下空间利用等，对工程地质环境产生了深刻的影响，进一步加剧了工程地质环境的复杂性和不确定性。在这种背景下，开展晋江市工程地质环境功能区划研究具有重要的现实意义。首先，通过对工程地质环境的系统评价和功能区划，可以全面了解晋江市不同区域的工程地质条件，明确各区域的优势和限制因素，为城市规划和建设提供科学依据，避免因不合理的工程建设导致地质灾害的发生，保障城市建设的安全。其次，合理的功能区划有助于优化土地资源配置，提高土地利用效率，促进城市的可持续发展。根据不同区域的工程地质条件，合理布局城市功能区，如居住区、商业区、工业区等，可以充分发挥土地的潜力，减少资源浪费，实现城市的高效发展。此外，工程地质环境功能区划还可以为环境保护和地质灾害防治提供指导，通过对地质灾害易发区的识别和划分，制定相应的防治措施，减少地质灾害对人民生命财产的威胁，保护生态环境的平衡。地理信息系统（GIS）作为一种强大的空间分析工具，具有数据管理、空间分析、可视化表达等多种功能，能够有效地处理和分析大量的地理空间数据。在工程地质环境功能区划研究中，应用GIS技术可以将各种地质数据、地形数据、工程建设数据等进行整合和分析，直观地展示工程地质环境的空间分布特征和变化规律，提高研究的准确性和科学性。因此，基于GIS技术开展晋江市工程地质环境功能区划研究，具有重要的理论和实践价值。1.2国内外研究现状1.2.1工程地质环境研究现状工程地质环境研究是地质学与工程学交叉的重要领域，旨在探究地质环境与人类工程活动之间的相互作用和影响。国外在这方面起步较早，早在20世纪中叶，随着大规模工程建设的开展，工程地质环境问题逐渐受到关注。早期研究主要集中在对地质灾害的认识和简单的工程地质条件评估，如滑坡、泥石流等地质灾害对工程建设的影响。随着研究的深入，逐渐发展到对工程地质环境系统的全面分析，包括地质构造、岩土体性质、地下水条件等多方面因素对工程的综合影响。近年来，国外工程地质环境研究在理论和技术上取得了显著进展。在理论方面，强调工程地质环境的动态演化和系统分析，将工程地质环境视为一个复杂的开放系统，研究其在自然因素和人类活动双重作用下的演变规律。在技术手段上，不断引入先进的探测技术和分析方法，如地球物理勘探技术、数值模拟技术等，以提高对工程地质环境的认识和评价精度。例如，利用高精度的地球物理探测技术，可以更准确地查明地下地质结构和岩土体性质；通过数值模拟方法，可以预测工程建设对地质环境的影响，为工程设计和决策提供科学依据。国内工程地质环境研究虽然起步相对较晚，但发展迅速。20世纪70年代以来，随着我国经济建设的快速发展，大量基础设施工程的兴建促使工程地质环境研究不断深入。早期主要是学习和借鉴国外的理论和方法，结合国内工程实践开展研究。近年来，国内在工程地质环境研究方面取得了一系列具有自主知识产权的成果，在理论、方法和技术上都有了很大的创新。例如，在地质灾害防治方面，提出了一系列适合我国国情的防治理论和技术体系，有效减少了地质灾害造成的损失；在工程地质勘察方面，发展了多种先进的勘察技术和方法，提高了勘察的效率和质量。1.2.2工程地质环境功能区划研究现状工程地质环境功能区划是在工程地质环境研究的基础上，根据区域工程地质条件和工程建设需求，将研究区域划分为不同功能区的过程。国外在工程地质环境功能区划方面的研究较为系统，注重多学科交叉和综合评价。在评价指标体系方面，涵盖了地质、地形、水文、生态等多个方面的因素，采用层次分析法、模糊综合评价法等多种数学方法进行综合评价。在功能区划分上，通常结合城市规划、土地利用规划等进行，以实现工程建设与地质环境保护的协调发展。例如，美国在城市规划中，充分考虑工程地质环境因素，通过功能区划合理布局城市建设项目，避免在地质灾害易发区进行大规模建设。国内工程地质环境功能区划研究在近几十年也取得了丰硕成果。学者们针对不同地区的特点，开展了大量的实证研究。在评价指标选取上，结合我国的实际情况，注重考虑地质灾害、岩土体工程性质、地下水等因素对工程建设的影响。在评价方法上，不断创新和改进，如引入神经网络、灰色关联分析等方法，提高评价的准确性和科学性。在功能区划分结果应用方面，为城市规划、土地利用规划、工程建设选址等提供了重要依据。例如，在一些大城市的规划中，根据工程地质环境功能区划结果，合理确定不同区域的开发强度和建设类型，有效降低了工程建设风险，保护了地质环境。1.2.3GIS技术在相关领域的应用现状GIS技术作为一种强大的空间分析工具，在工程地质环境研究和功能区划中得到了广泛应用。国外在GIS技术应用方面处于领先地位，早在20世纪70年代就开始将GIS技术应用于地质领域。目前，GIS技术已广泛应用于地质数据管理、地质灾害分析、工程地质评价等多个方面。例如，利用GIS的空间分析功能，可以对地质灾害的分布规律进行分析，预测灾害的发展趋势；通过建立地质模型，结合GIS技术可以对工程地质条件进行三维可视化展示，为工程设计提供直观的依据。在国内，GIS技术在工程地质环境领域的应用也日益普及。随着GIS技术的不断发展和完善，国内学者在应用方面进行了大量的探索和实践。在工程地质环境功能区划中，利用GIS技术可以快速、准确地处理和分析海量的地质数据，实现评价指标的量化和空间分析，提高功能区划的效率和精度。例如，通过将地质、地形、水文等多源数据导入GIS系统，进行空间叠加分析和统计分析，可以直观地展示工程地质环境的空间分布特征，为功能区划分提供科学依据。同时，基于GIS技术开发的工程地质环境信息管理系统，实现了地质数据的信息化管理和共享，方便了数据的查询和更新，为工程地质环境研究和管理提供了有力支持。1.2.4存在的问题尽管国内外在工程地质环境研究、功能区划以及GIS技术应用方面取得了诸多成果，但仍存在一些问题有待解决。在评价指标体系方面，虽然目前已经考虑了多个方面的因素，但不同地区的地质条件和工程建设需求差异较大，现有的评价指标体系在通用性和针对性上还存在不足，需要进一步完善和优化。在评价方法上，虽然多种数学方法被广泛应用，但不同方法的适用条件和局限性不同，在实际应用中如何选择合适的评价方法，以及如何提高评价结果的可靠性和准确性，仍是需要深入研究的问题。在GIS技术应用方面，虽然GIS技术为工程地质环境研究和功能区划提供了强大的支持，但目前在数据质量、数据更新、模型构建等方面还存在一些问题。例如，地质数据的获取和整理工作较为繁琐，数据质量参差不齐，影响了GIS分析结果的准确性；数据更新不及时，难以反映工程地质环境的动态变化；在构建地质模型和分析模型时，缺乏统一的标准和规范，导致模型的可靠性和可重复性较差。此外，在工程地质环境功能区划结果的应用和实施方面，与城市规划、土地利用规划等的衔接还不够紧密，如何将功能区划结果更好地融入到实际的工程建设和规划管理中，还需要进一步加强研究和实践。1.3研究目标与内容1.3.1研究目标本研究旨在基于GIS技术，全面、系统地开展晋江市工程地质环境功能区划研究。通过对晋江市工程地质环境条件的深入调查和分析，建立科学合理的评价指标体系和评价模型，准确划分工程地质环境功能区，为晋江市的城市规划、土地利用、工程建设等提供科学依据，促进城市的可持续发展。具体目标如下：全面分析工程地质环境条件：详细收集和整理晋江市的地形地貌、地层岩性、地质构造、岩土体特征、地质灾害、地下水等工程地质环境资料，利用GIS强大的数据管理和分析功能，对这些资料进行综合分析，深入了解晋江市工程地质环境的现状和特征。建立科学的评价指标体系和模型：根据晋江市工程地质环境的特点和城市发展需求，遵循科学性、系统性、实用性等原则，选取合适的评价指标，运用层次分析法、熵值法等方法确定各指标的权重，构建综合评价模型，实现对工程地质环境适宜性的定量评价。准确划分工程地质环境功能区：依据综合评价结果，结合晋江市的城市规划和土地利用现状，将晋江市划分为不同的工程地质环境功能区，明确各功能区的范围、功能定位和发展方向，为城市建设和管理提供指导。为城市规划和建设提供科学依据：将工程地质环境功能区划结果与城市规划、土地利用规划等相结合，为城市的合理布局、基础设施建设、地质灾害防治等提供科学建议，促进城市建设与工程地质环境的协调发展，保障城市的安全和可持续发展。1.3.2研究内容研究区地质环境概况分析：全面收集晋江市的自然地理条件数据，包括地形地貌、河流水系等信息。详细梳理区域地质条件，涵盖地层概况、地质构造以及侵入岩体等方面。深入研究岩土体特征，分析岩体和土体的特性。调查地质灾害特征，掌握其类型、分布和发生规律。探究地下水特征，包括松散岩类孔隙潜水、风化带孔隙裂隙水和基岩裂隙水等，为后续的功能区划研究奠定基础。工程地质功能区划评价模型构建：按照科学性、全面性、独立性等原则，选取地形坡度、岩土体类型、地质灾害危险性、地下水水位等多个评价指标，构建晋江市工程地质环境功能区划评价指标体系。运用层次分析法和熵值法相结合的方式，确定各评价指标的权重，以体现不同指标对工程地质环境功能区划的影响程度。构建综合评价模型，对晋江市工程地质环境适宜性进行定量评价，为功能区划分提供数据支持。工程地质功能区划实现：采用格网剖分的方法，将晋江市研究区域划分为多个评价单元，确保评价的全面性和准确性。计算每个评价单元的各项评价指标值，并根据构建的综合评价模型，计算各评价单元的综合得分，以此作为功能区划分的依据。根据综合得分结果，运用聚类分析等方法，将研究区域划分为不同的工程地质功能区，并对各功能区的特征和功能进行详细阐述和分析。结果验证与分析：运用实地调查、历史数据对比等方法，对工程地质功能区划结果进行验证，确保划分结果的准确性和可靠性。深入分析各功能区的工程地质条件和发展潜力，结合晋江市城市规划和土地利用规划，提出各功能区的合理开发利用建议和地质灾害防治措施，为城市建设和管理提供科学指导。1.4技术路线本研究以晋江市为研究区域，围绕工程地质环境功能区划展开，技术路线如图1-1所示。具体步骤如下：资料收集与整理：全面收集晋江市地形地貌、地层岩性、地质构造、岩土体特征、地质灾害、地下水等相关资料。对收集到的资料进行整理和分析，确保数据的准确性和完整性。通过实地调查和测量，对部分关键数据进行验证和补充，提高数据的可靠性。数据处理与分析：运用ArcGIS软件强大的数据处理功能，将收集到的各类数据进行数字化处理，转化为计算机可识别和处理的格式。利用ArcGIS的空间分析功能，对地形坡度、岩土体类型、地质灾害危险性、地下水水位等数据进行分析和计算，提取关键信息。结合数理统计方法，对数据进行统计分析，了解数据的分布特征和规律，为后续的评价和区划提供数据支持。评价指标体系构建：依据科学性、全面性、独立性等原则，从地形地貌、岩土体特性、地质灾害、地下水等方面选取评价指标，构建晋江市工程地质环境功能区划评价指标体系。采用层次分析法和熵值法相结合的方式，确定各评价指标的权重。层次分析法通过构建判断矩阵，确定各指标的相对重要性；熵值法根据指标数据的变异程度确定权重，两者结合使权重确定更加科学合理。综合评价模型构建：根据确定的评价指标体系和权重，构建晋江市工程地质环境适宜性综合评价模型。运用模糊综合评价法、灰色关联分析法等方法对评价单元进行综合评价，计算各评价单元的综合得分，以量化工程地质环境的适宜性。工程地质功能区划：采用格网剖分的方法，将晋江市研究区域划分为多个大小相等的评价单元，每个评价单元具有相同的面积和形状。根据综合评价模型计算得到的各评价单元综合得分，运用聚类分析等方法，将得分相近的评价单元划分为同一功能区，从而将研究区域划分为不同的工程地质功能区。对划分出的各功能区进行特征分析和功能定位，明确各功能区的优势和限制因素，提出合理的开发利用建议。结果验证与分析：运用实地调查、历史数据对比等方法，对工程地质功能区划结果进行验证。实地调查各功能区的工程地质条件，与区划结果进行对比分析，检查区划结果的准确性和可靠性。根据验证结果，对功能区划结果进行优化和调整，确保区划结果更加符合晋江市的实际情况。深入分析各功能区的工程地质条件和发展潜力，结合晋江市城市规划和土地利用规划，提出各功能区的合理开发利用建议和地质灾害防治措施，为城市建设和管理提供科学指导。图1-1技术路线图二、晋江市地质环境概况2.1自然地理条件2.1.1地形地貌晋江市地处福建省东南沿海，位于闽东南沿海大陆边缘坳陷变质带中部，其地势呈现出西北高、东南低的态势，地形地貌丰富多样，涵盖了山地、丘陵、台地和平原等多种类型。山地主要集中在晋江市的西北部和中部地区，是戴云山向东南沿海延伸的余脉。其中，西北部与鲤城区交界处的紫帽山海拔达到517.8米，为境内最高峰，山体巍峨，植被茂密，不仅是重要的自然景观，还对区域生态环境有着重要的调节作用。中部的灵源山海拔305.0米，山上古迹众多，文化底蕴深厚，其独特的地形地貌也为当地的旅游业发展提供了资源基础。这些山地的岩石主要由二长花岗岩、花岗闪长岩和金黑云花岗岩等组成，经过长期的地质作用和风化侵蚀，形成了陡峭的山峰、深邃的山谷等复杂地形。丘陵在晋江市分布较为广泛，其海拔相对较低，一般在几十米到一两百米之间。这些丘陵起伏和缓，形态各异，与山地相互交错。丘陵地区的土壤类型多样，主要包括红壤、黄壤等，适合发展经济林和果树种植，如龙眼、荔枝等果树在丘陵地带广泛种植，成为当地农业经济的重要组成部分。台地是晋江市面积最大的地貌单元，主要分布在中部和南部地区。台地地势相对平坦，坡度较小，有利于大规模的土地开发和建设。台地的岩性主要为第四纪松散沉积物，土层较厚，土壤肥沃，是重要的农业生产区域，种植着水稻、蔬菜等农作物。同时，由于其地势平坦开阔，也为工业发展和城市建设提供了良好的基础条件，许多工业园区和城镇建设在台地之上。平原主要分布在东北部靠近泉州湾和西南部围头湾一带。晋东平原是晋江重要的平原区域，由九十九溪、晋江及海浪冲积而成，属于泉州平原的构成部分。平原地区地势低平，水源充足，土壤肥沃，是晋江市人口密集、经济发达的区域。这里交通便利，基础设施完善，集中了大量的工业企业、商业中心和居民区，是晋江市的经济核心地带。同时，平原地区也是重要的农业产区，主要种植水稻、小麦等粮食作物，为保障当地的粮食供应发挥着重要作用。晋江市的海岸线长121千米，海岸地貌丰富多样。其中，海滩是海岸地貌的重要组成部分，主要由砂质沉积物组成，沙滩细腻柔软，是休闲旅游的好去处，如深沪湾的海滩吸引了众多游客前来观光度假。海湾则为海洋生物提供了良好的栖息环境，也是渔业和港口发展的重要区域，围头湾是重要的对台贸易港口，在区域经济交流中发挥着重要作用。此外，晋江还有一些海岛，这些海岛生态环境独特，拥有丰富的自然资源和生物多样性，具有重要的生态保护和旅游开发价值。2.1.2河流水系晋江市的河流水系较为发达，主要河流包括晋江、九十九溪等，这些河流对晋江市的经济发展、生态环境和居民生活都有着至关重要的影响。晋江是泉州地区的第一大河，也是福建省第三大河，被誉为泉州人的母亲河。它发源于戴云山脉东南坡，在中上游分为东西两溪，东溪发源于永春县呈祥乡云路村，流域面积1917平方千米，河长120千米；西溪发源于安溪县感德乡挑舟村，流域面积3101平方千米，河长14千米。两支流于南安市丰州镇井兜村双溪口汇合始称晋江，并于晋江市淄滨纳入九十九溪至丰泽区浔浦入海，全长约162千米，流域面积5629平方千米。晋江水量变化较大，含沙量较福建省其他河流为多，多年平均径流量为53.49亿立方米。其流域为亚热带季风气候，温暖湿润，年均气温20.8℃，年降雨量自西向东南递减，平均在1600毫米左右，集中于4-9月，占全年的76%-80%。晋江年平均径流深为950毫米，径流量约47.8亿立方米，汛期径流量占全年的75%，其中5-7月占60%。由于晋江下游地势低洼，面临东海，台风频繁，涨潮顶托，以及人多地少、水土流失严重，每年逢6-9月台风季节，经常发生洪潮灾害，水旱交错。九十九溪是晋江市境内最长的河流，发源于南安市大旗尾山，由官桥镇流入晋江市注入泉州湾。在晋江市境内长20千米，流域面积137平方千米。九十九溪的河水主要来源于降水和上游的山涧溪流，其水量受季节影响较大，雨季时水量充沛，旱季时水量相对减少。河流两岸地势平坦，土壤肥沃，是重要的农业灌溉水源，灌溉着周边大量的农田，为当地的农业生产提供了有力保障。同时，九十九溪也是晋江市生态环境的重要组成部分，河流两岸植被丰富，为众多生物提供了栖息和繁衍的场所。除了晋江和九十九溪，晋江市还有一些独流入海的时令溪流，约有19条。这些溪流大多短小，受降水影响明显，在雨季时水流较大，而在旱季时可能会出现干涸的情况。它们虽然规模较小，但在调节区域地表径流、补充地下水等方面也发挥着一定的作用。晋江市内还有多处湖泊和水库，如东湖、西湖、科任水库等。这些湖泊和水库在调节洪水、提供水源、改善生态环境等方面发挥着重要作用。例如，在雨季时，湖泊和水库可以储存多余的水量，减轻洪水对下游地区的压力；在旱季时，它们又可以为周边地区提供生产和生活用水，保障居民的正常生活和工农业生产的顺利进行。同时，湖泊和水库周边的湿地生态系统也是许多珍稀鸟类和动植物的栖息地，对于维护区域生物多样性具有重要意义。2.2区域地质条件2.2.1地层概况晋江市的地层发育较为复杂，从老到新主要包括中生代侏罗系、白垩系，新生代古近系、新近系以及第四系。中生代侏罗系主要出露于晋江市的西北部山区，岩性主要为一套陆相碎屑岩沉积，包括砂岩、页岩、粉砂岩等。这些岩石在长期的地质作用下，形成了较为致密的结构，抗压强度较高，但由于经历了多次构造运动，岩石中节理、裂隙较为发育，对工程建设存在一定影响。例如，在进行山区道路建设时，需要充分考虑岩石的节理裂隙分布，避免因岩体破碎导致边坡失稳。白垩系地层在晋江市也有一定分布，主要岩性为火山碎屑岩和熔岩。火山碎屑岩如凝灰岩，具有孔隙率较大、强度较低的特点；而熔岩如流纹岩，结构致密，强度相对较高。在工程建设中，对于白垩系地层的处理需要根据具体岩性进行分析。如在进行基础工程建设时，对于凝灰岩地层，可能需要采取加固措施来提高地基的承载能力。新生代古近系和新近系地层主要分布在晋江市的中部和南部地区，多为一套河湖相沉积，岩性主要包括泥岩、砂岩、砾岩等。泥岩具有遇水易软化、强度降低的特性，在工程建设中，如果基础位于泥岩地层上，需要注意防水和地基处理，防止因泥岩软化导致基础沉降。砂岩和砾岩的强度相对较高，但砾岩的颗粒较大，在进行地下工程施工时，需要考虑其透水性和稳定性。第四系地层在晋江市广泛分布，主要包括冲积层、洪积层、残积层和风积层等。冲积层主要分布在河流两岸和河口地区，由河流携带的泥沙堆积而成，岩性主要为砂、砾石和黏土，其颗粒大小和分选性因河流的流速和流量不同而有所差异。在这些地区进行工程建设时，需要考虑地基的承载力和不均匀沉降问题，例如在建设桥梁基础时，需要对冲积层进行详细的勘察和分析，确保基础的稳定性。洪积层分布在山前地带，由洪水携带的碎屑物质堆积而成，其颗粒较粗，透水性较强。残积层是岩石风化后残留在原地的产物，岩性与下伏基岩密切相关，其工程性质受风化程度影响较大。风积层主要分布在沿海地区，由风力搬运的砂粒堆积而成，如沿海的沙丘，其结构松散，在工程建设中需要进行加固处理。2.2.2地质构造晋江市位于闽东南沿海大陆边缘坳陷变质带中部，地质构造较为复杂，主要受东北新华系结构控制。由于地处长乐-南澳大断裂中段，境内发育有青阳-安海、西坑-古厝、祥芝-围头三条断裂带。这些断裂带对晋江市的工程地质条件产生了重要影响。断裂带附近的岩石受到强烈的构造应力作用，岩体破碎，节理裂隙发育，导致岩石的完整性和强度降低。这使得在断裂带附近进行工程建设时，地基的稳定性较差，容易发生地基沉降、塌陷等问题。例如，在进行高层建筑建设时，如果基础位于断裂带附近，需要采取特殊的地基处理措施，如采用桩基础等，以提高地基的承载能力和稳定性。同时，断裂带还可能影响地下水的分布和径流，导致地下水水位变化较大，给工程建设中的地下水处理带来困难。褶皱构造在晋江市也有一定表现，主要为一些规模较小的褶皱。褶皱构造会使地层发生弯曲变形，导致地层的产状发生变化。在进行工程建设时，需要准确掌握地层的产状，以合理设计工程结构。例如，在进行隧道工程建设时，如果隧道穿越褶皱构造区域，需要根据地层的弯曲情况，合理确定隧道的走向和支护方式，避免因地层压力不均匀导致隧道坍塌。此外，褶皱构造还可能影响岩石的力学性质，使得褶皱核部和翼部的岩石强度和稳定性存在差异，在工程建设中需要分别进行评估和处理。2.2.3侵入岩体晋江市的侵入岩体主要为花岗岩类，包括二长花岗岩、花岗闪长岩和金黑云花岗岩等。这些侵入岩体主要分布在晋江市的西北部和中部山区。花岗岩类岩石具有结构致密、强度高、抗风化能力强等优点，是良好的建筑材料和地基基础。在工程建设中，对于大型建筑物的基础，如果位于花岗岩岩体上，一般可以直接利用其作为天然地基，能够承受较大的上部荷载。例如，在建设大型工业厂房时，花岗岩地基可以为厂房提供稳定的支撑。然而，花岗岩岩体中也可能存在一些节理、裂隙等结构面，这些结构面会在一定程度上影响岩体的完整性和强度。在进行工程建设前，需要对花岗岩岩体进行详细的勘察，评估结构面的分布和性质，采取相应的工程措施，如对节理裂隙进行灌浆处理等，以确保岩体的稳定性。同时，花岗岩的风化程度也会对其工程性质产生影响，风化程度较深的花岗岩，其强度会降低，在工程建设中需要进行特殊处理。2.3岩土体特征2.3.1岩体特征晋江市的岩体主要为花岗岩类，包括二长花岗岩、花岗闪长岩和金黑云花岗岩等，这些岩体广泛分布于西北部和中部山区。花岗岩类岩石具有较高的抗压强度和抗剪强度，一般抗压强度可达100-200MPa，抗剪强度在10-20MPa之间。其结构致密，矿物结晶良好，颗粒间的胶结力较强，使得岩体具有较好的整体性和稳定性。这种良好的力学性质使得花岗岩成为大型工程建设的理想基础材料，如在建设大型桥梁基础、高层建筑基础时，花岗岩地基能够提供稳定的承载能力。岩体中的结构面特征对其工程性质有着重要影响。晋江市的花岗岩岩体中存在着节理、裂隙等结构面，这些结构面的分布具有一定的规律性。节理的走向主要有北北东向、北东向和北西向等，其中北北东向节理较为发育，其延伸长度一般在数米至数十米之间。裂隙的密度和张开度因区域而异，在构造活动较为强烈的地区，裂隙密度较大，张开度也相对较大。结构面的存在降低了岩体的完整性和强度，使得岩体在受到外力作用时，容易沿着结构面发生破坏。例如，在进行露天采矿时，节理和裂隙会影响岩体的爆破效果和边坡稳定性；在隧道工程中，结构面可能导致岩体的坍塌，因此在工程建设中需要对结构面进行详细的勘察和分析，并采取相应的加固措施。此外，岩体的风化程度也是影响其工程性质的重要因素。晋江市的花岗岩岩体风化程度在不同区域有所差异，一般来说，山区顶部和山坡地带的岩体风化程度较深，而山谷和低洼地带的岩体风化程度相对较浅。风化作用使得岩体的矿物成分发生改变，结构变得疏松，强度降低。全风化花岗岩的抗压强度一般在1-5MPa之间，中风化花岗岩的抗压强度在10-50MPa之间。在工程建设中，对于风化程度较深的岩体，需要进行特殊的处理，如采用换填、灌浆等方法来提高其承载能力和稳定性。2.3.2土体特征晋江市的土体类型丰富多样，主要包括残积土、冲积土、洪积土和风积土等。残积土是岩石经过长期风化作用后残留在原地的产物，其岩性与下伏基岩密切相关。在花岗岩分布区，残积土主要为花岗岩残积土，其颗粒组成以黏性土和砂质土为主，含有少量的砾石。花岗岩残积土的工程性质受风化程度影响较大，一般来说，风化程度越深，其压缩性越高，强度越低。残积土的含水量较高，孔隙比大，透水性较差，在工程建设中需要注意其地基承载力和稳定性问题。例如，在进行建筑物基础设计时，需要对残积土地基进行详细的勘察和分析，根据其工程性质采取相应的基础形式和处理措施。冲积土主要分布在河流两岸和河口地区，是由河流携带的泥沙堆积而成。其颗粒大小和分选性因河流的流速和流量不同而有所差异，一般靠近河流上游，颗粒较粗，以砾石和粗砂为主；靠近河流下游，颗粒较细，以细砂和粉土为主。冲积土的工程性质较好，具有较高的承载力和较低的压缩性。在进行工程建设时，冲积土可以作为一般建筑物的天然地基，但需要注意其不均匀沉降问题，特别是在颗粒粗细变化较大的区域。例如，在建设道路工程时，需要对冲积土地基进行压实处理，以提高其承载能力和稳定性，防止道路出现不均匀沉降和开裂。洪积土分布在山前地带，是由洪水携带的碎屑物质堆积而成。其颗粒较粗，多为砾石、粗砂和碎石等，分选性较差，具有较大的孔隙率和透水性。洪积土的承载力较高，但由于其颗粒间的胶结力较弱，在受到振动或水流作用时，容易发生变形和失稳。在工程建设中，对于洪积土地基，需要进行适当的加固处理，如采用强夯法、砂石桩法等，以提高其稳定性。例如，在建设大型工业厂房时，如果地基位于洪积土上，需要对地基进行加固处理，确保厂房的安全稳定。风积土主要分布在沿海地区，是由风力搬运的砂粒堆积而成，如沿海的沙丘。风积土的颗粒均匀，以细砂为主，结构松散，孔隙率大，透水性强。风积土的承载力较低，压缩性较高，在工程建设中需要进行加固处理。常用的加固方法有土工织物加筋、压实法、化学加固法等。例如，在沿海地区进行建筑工程时，对于风积土地基，可采用土工织物加筋的方法，提高地基的承载能力和稳定性；也可以通过压实法，增加土体的密实度，降低其压缩性。2.4地质灾害特征晋江市的地质灾害类型多样，主要包括滑坡、泥石流、地面沉降等，这些地质灾害的发生给当地的人民生命财产安全和经济发展带来了严重威胁。滑坡是晋江市较为常见的地质灾害之一，主要分布在山区和丘陵地带，如西北部的紫帽山、中部的灵源山等周边区域。这些地区地形起伏较大，坡度较陡，岩土体在重力作用下，容易沿着软弱结构面发生滑动。滑坡的成因主要与地形地貌、岩土体性质、降水以及人类工程活动等因素密切相关。从地形地貌来看，山区和丘陵的陡峭地形为滑坡的发生提供了势能条件，坡度越大，岩土体下滑的可能性就越大。岩土体性质方面，若岩土体的抗剪强度较低，如一些风化严重的岩石或土体，在受到外力作用时，就容易发生变形和滑动。降水是滑坡发生的重要诱发因素，晋江市属于亚热带季风气候，降水集中且强度大，尤其是在雨季和台风季节，大量的雨水渗入地下，使岩土体的含水量增加，重度增大，同时降低了岩土体的抗剪强度，从而触发滑坡。人类工程活动如切坡建房、修路等，破坏了山体的原有稳定性，也增加了滑坡发生的风险。例如，在山区进行道路建设时，不合理的切坡会使山体的边坡变陡，改变了岩土体的应力分布，在降雨等因素的作用下，容易引发滑坡灾害。泥石流也是晋江市需要重点防范的地质灾害，多发生在沟谷地区，这些区域地形狭窄，沟床比降大，有利于泥石流的形成和汇聚。泥石流的形成通常需要丰富的松散固体物质、充足的水源和较大的地形坡度等条件。晋江市的山区和丘陵地区岩石风化强烈，加之人类工程活动如采矿、修路等产生了大量的废渣和弃土，为泥石流的发生提供了丰富的固体物质来源。在降水集中的季节，大量的雨水迅速汇聚成地表径流，与松散固体物质混合，形成强大的泥石流。当泥石流发生时，具有强大的冲击力和破坏力，会冲毁房屋、道路、桥梁等基础设施，掩埋农田，造成人员伤亡和财产损失。地面沉降在晋江市的部分地区也有出现，主要分布在城市中心区和沿海平原等人口密集、经济活动频繁的区域。地面沉降的成因较为复杂，主要与地下水开采、工程建设以及地质构造等因素有关。随着城市的发展，对水资源的需求不断增加，大量抽取地下水导致地下水位下降，土体有效应力增加，引起土层压缩，从而导致地面沉降。工程建设方面，大规模的高层建筑建设和地下空间开发，改变了地基土的应力状态，也可能引发地面沉降。此外，地质构造活动导致的地层变形，也会对地面沉降产生一定影响。地面沉降会导致地面建筑物开裂、倾斜，道路、桥梁等基础设施损坏，还会使沿海地区的防潮能力降低，加剧洪涝灾害的威胁。晋江市的地质灾害分布具有明显的地域性特征，山区和丘陵地带以滑坡、泥石流等地质灾害为主，而城市中心区和沿海平原则主要面临地面沉降等问题。这些地质灾害的发生不仅对当地的生态环境造成破坏，还严重影响了城市的可持续发展和居民的生活质量，因此，加强地质灾害的监测、预警和防治工作至关重要。2.5地下水特征晋江市的地下水类型主要包括松散岩类孔隙潜水、风化带孔隙裂隙水和基岩裂隙水。松散岩类孔隙潜水主要分布在沿海平原和河谷地区，含水层岩性主要为砂、砾石等，其透水性较好，富水性较强。该类潜水的水位变化受降水和潮汐影响较大，雨季时，降水充沛，地下水水位上升；旱季时，降水减少，水位下降。在靠近海岸的区域，潮汐作用也会使地下水水位产生周期性变化。由于含水层颗粒较粗，孔隙大，其水质一般较好，矿化度较低，适合作为生活用水和农业灌溉用水。但在一些人口密集、工业发达的地区，由于人类活动的影响，如工业废水排放、生活污水下渗等，可能导致地下水受到污染，水质变差。风化带孔隙裂隙水赋存于岩石的风化带中，在晋江市的丘陵和山区分布较为广泛。风化带岩石破碎，孔隙和裂隙发育，为地下水的储存和运移提供了空间。这类地下水的水位和水量受岩石风化程度、地形地貌和降水等因素的影响。一般来说，风化程度越深，孔隙裂隙越发育，地下水的储存量和补给量就越大。在地形低洼处，地下水容易汇聚，水位相对较高。其水质与岩石成分和风化程度有关，部分地区的风化带孔隙裂隙水可能含有较多的矿物质，矿化度相对较高。基岩裂隙水主要存在于基岩的裂隙中，分布在晋江市的山区和丘陵地带。基岩裂隙的发育程度和连通性决定了地下水的赋存和运移条件。在构造运动强烈的区域，基岩裂隙较为发育，地下水相对丰富；而在构造稳定的区域，裂隙较少，地下水含量也较少。基岩裂隙水的水位变化相对较为稳定，但在降水集中的季节，也会有所上升。其水质一般较好，受外界污染的影响较小，但由于其埋藏较深，开采难度相对较大。地下水对晋江市的工程建设有着重要影响。在进行基础工程建设时，地下水水位过高可能导致地基土的强度降低，增加基础的沉降量，甚至引发地基失稳。例如，在沿海平原地区进行高层建筑建设时，如果地下水水位较高，地基土处于饱和状态，其抗剪强度会明显下降，容易使基础产生不均匀沉降，导致建筑物倾斜、开裂。此外，地下水的腐蚀性也会对工程设施造成损害，如一些含有酸性物质或盐分的地下水，会腐蚀地下管道、基础钢筋等，缩短工程设施的使用寿命。在地下工程施工中，如隧道、地下商场等建设，地下水的涌入会增加施工难度和安全风险，可能引发坍塌等事故。2.6本章小结本章全面且系统地分析了晋江市的地质环境概况，涵盖自然地理条件、区域地质条件、岩土体特征、地质灾害特征以及地下水特征等多个方面。在自然地理条件方面，晋江市地势西北高、东南低，地形地貌类型丰富，包括山地、丘陵、台地和平原等，不同地形地貌在区域内的分布具有一定规律，且对工程建设和生态环境有着不同程度的影响。同时，晋江市河流水系发达，主要河流晋江和九十九溪等在区域水资源调配、生态调节和经济发展中发挥着重要作用。区域地质条件复杂，地层发育多样，从老到新涵盖中生代侏罗系、白垩系，新生代古近系、新近系以及第四系。地质构造受东北新华系结构控制，境内存在多条断裂带和褶皱构造，这些构造对岩体的完整性和稳定性产生影响，进而影响工程建设的选址和设计。侵入岩体主要为花岗岩类，其力学性质和风化程度在不同区域存在差异，对工程建设的影响也各不相同。岩土体特征方面，岩体主要为花岗岩类，具有较高的抗压强度和抗剪强度，但结构面和风化程度会降低其完整性和强度。土体类型多样，包括残积土、冲积土、洪积土和风积土等，不同土体的工程性质各异，在工程建设中需要根据其特性采取相应的处理措施。地质灾害类型多样，滑坡、泥石流和地面沉降等灾害在不同区域有不同的分布特征和成因。滑坡主要分布在山区和丘陵地带，与地形、岩土体性质、降水和人类工程活动等因素密切相关；泥石流多发生在沟谷地区，受松散固体物质、水源和地形坡度等条件影响；地面沉降主要出现在城市中心区和沿海平原，与地下水开采、工程建设和地质构造等因素有关。地下水类型包括松散岩类孔隙潜水、风化带孔隙裂隙水和基岩裂隙水，其分布、水位变化和水质受地形地貌、地层岩性和降水等因素影响，对工程建设中的基础稳定性、地下工程施工和工程设施耐久性等方面具有重要影响。本章对晋江市地质环境概况的详细分析，为后续工程地质功能区划评价模型的构建以及功能区的划分提供了坚实的基础，有助于深入了解晋江市工程地质环境的特点和规律，为城市规划和建设提供科学依据。三、工程地质功能区划评价模型构建3.1评价指标体系确立3.1.1指标体系确立原则科学性原则：评价指标体系应基于科学的理论和方法构建，准确反映工程地质环境的本质特征和内在规律。指标的选取要具有明确的科学依据，能够客观地衡量工程地质环境的各个方面，避免主观随意性。例如，在选取岩土体性质指标时，应依据岩土力学理论，选择能够反映岩土体强度、变形特性等关键力学性质的参数，如岩土体的抗压强度、抗剪强度、压缩模量等，确保指标能够准确地反映岩土体在工程建设中的稳定性和承载能力。系统性原则：工程地质环境是一个复杂的系统，包含多个相互关联、相互影响的因素。评价指标体系应从系统的角度出发，全面考虑各种因素对工程地质环境的影响，形成一个有机的整体。不仅要涵盖地形地貌、岩土体性质、地质灾害、地下水等主要方面，还要考虑各方面因素之间的相互作用和协同效应。例如，地形地貌会影响地下水的分布和径流，进而影响岩土体的稳定性；地质灾害的发生与岩土体性质、地形地貌以及人类工程活动等多种因素密切相关。因此，在构建指标体系时，要充分考虑这些因素之间的内在联系，使指标体系能够全面、系统地反映工程地质环境的整体状况。可操作性原则：评价指标应具有可获取性和可量化性，便于在实际工作中进行数据采集和分析。选取的指标应能够通过现有的技术手段和方法进行测量或估算，数据来源可靠，且能够用具体的数值进行表达。例如，地形坡度可以通过地形数据利用GIS软件进行计算得到；岩土体类型可以通过地质勘察资料直接确定；地质灾害危险性可以通过对历史灾害数据的统计分析以及地质条件的评估进行量化。同时，指标的计算方法应简单明了，易于操作，避免过于复杂的计算过程和难以理解的概念，以保证评价工作的高效性和准确性。独立性原则：各评价指标之间应具有相对独立性，避免指标之间存在过多的重叠或相关性。如果指标之间相关性过高，会导致信息重复，影响评价结果的准确性和可靠性。在选取指标时，要对指标进行严格的筛选和分析，确保每个指标都能独立地反映工程地质环境的某一方面特征。例如，在选择地质灾害相关指标时，滑坡、泥石流等不同类型的地质灾害指标应分别选取，避免将滑坡和泥石流的发生频率合并为一个指标，因为它们的形成机制和影响因素有所不同，分别选取可以更全面地反映地质灾害的情况。动态性原则：工程地质环境是一个动态变化的系统，受到自然因素和人类工程活动的双重影响。评价指标体系应具有一定的动态性，能够反映工程地质环境的变化趋势。随着时间的推移和工程建设的进行，地形地貌可能会因人类活动而改变，岩土体性质可能会因地下水的变化或工程荷载的作用而发生变化，地质灾害的发生频率和强度也可能会有所改变。因此，在构建指标体系时，要考虑到这些动态变化因素，预留一定的指标调整空间，以便能够及时更新和完善指标体系，使其更好地适应工程地质环境的变化。3.1.2评价指标体系确立基于上述原则，结合晋江市的工程地质环境特点，确定以下评价指标：地形地貌指标：地形坡度是影响工程建设的重要地形地貌因素之一。较大的地形坡度会增加工程建设的难度和成本，如在进行道路建设时，坡度较大的区域需要进行大量的土石方开挖和填方工程，同时也会增加边坡失稳的风险。通过对晋江市地形数据的分析，利用GIS软件的空间分析功能，计算出不同区域的地形坡度，并将其划分为不同的等级，如小于5°、5°-15°、15°-25°、大于25°等，以反映地形坡度对工程建设的影响程度。地形起伏度反映了地形的高低变化程度，起伏度较大的区域，工程建设的场地平整难度较大，基础设施建设成本也会相应增加。利用数字高程模型（DEM）数据，通过计算一定范围内的地形高差，得到地形起伏度指标，同样将其划分为不同等级，用于评价地形地貌条件对工程建设的适宜性。岩土体性质指标：岩土体类型直接影响着地基的承载能力和稳定性。晋江市的岩土体类型多样，包括花岗岩、砂岩、泥岩、砂土、黏土等。不同类型的岩土体具有不同的物理力学性质，如花岗岩强度高、稳定性好，适合作为大型建筑物的基础；而泥岩遇水易软化，强度降低，在工程建设中需要采取特殊的处理措施。根据岩土体的工程性质，将其分为优良、良好、中等、较差四个等级，为工程建设提供基础依据。岩土体的抗压强度和抗剪强度是衡量其力学性能的重要指标，直接关系到工程地基的承载能力和抗变形能力。通过室内试验和现场测试等方法，获取不同岩土体的抗压强度和抗剪强度数据，并将其纳入评价指标体系，以量化岩土体性质对工程建设的影响。地质灾害指标：地质灾害危险性是工程地质环境功能区划中需要重点考虑的因素之一。晋江市主要的地质灾害类型有滑坡、泥石流、地面沉降等。通过对地质灾害的历史数据进行统计分析，结合地质条件、地形地貌和人类工程活动等因素，利用层次分析法、模糊综合评价法等方法，对不同区域的地质灾害危险性进行评价，划分为低危险区、较低危险区、中等危险区、较高危险区和高危险区五个等级，为工程建设选址和地质灾害防治提供参考。地质灾害的发生频率和规模也是衡量其危险性的重要指标。统计不同类型地质灾害在不同区域的发生次数，以及灾害发生时造成的破坏范围和损失程度，将这些数据作为地质灾害指标的补充，更全面地反映地质灾害对工程地质环境的影响。地下水指标：地下水水位的高低对工程建设有着重要影响。水位过高会导致地基土的含水量增加，强度降低，增加基础的沉降量，甚至引发地基失稳；水位过低则可能影响工程的供水条件。通过对晋江市地下水水位监测数据的分析，结合工程建设的要求，确定不同的水位区间，如水位埋深小于1m、1m-3m、3m-5m、大于5m等，并评估不同水位区间对工程建设的影响程度。地下水的腐蚀性也是需要考虑的重要因素。地下水含有多种化学成分，如硫酸根离子、氯离子、氢离子等，当这些成分含量较高时，会对工程设施中的金属构件、混凝土等产生腐蚀作用，缩短工程设施的使用寿命。通过对地下水水质的检测分析，确定其腐蚀性等级，如弱腐蚀性、中等腐蚀性、强腐蚀性等，将其作为评价指标之一，为工程建设中的防腐设计提供依据。通过以上评价指标体系的构建，能够全面、系统地反映晋江市工程地质环境的特征，为后续的工程地质功能区划评价模型构建和功能区划分提供科学依据。3.2评价方法确定3.2.1层次分析法简介层次分析法（AnalyticHierarchyProcess，简称AHP）是美国运筹学家T.L.Saaty在20世纪70年代提出的一种定性与定量相结合的多目标决策方法。该方法将复杂的决策问题分解为多个层次，通过对各层次元素的两两比较，确定其相对重要性，进而计算出各元素的权重。其基本原理和步骤如下：建立层次结构模型：将决策问题分解为目标层、准则层和方案层等不同层次。在晋江市工程地质环境功能区划研究中，目标层为工程地质环境功能区划，准则层包括地形地貌、岩土体性质、地质灾害、地下水等方面的评价指标，方案层则为不同的工程地质功能区划分方案。以地形地貌指标为例，在层次结构模型中，它作为准则层的一个重要因素，影响着目标层工程地质环境功能区划的结果，同时又受到其下具体的地形坡度、地形起伏度等子因素的影响。构造判断矩阵：针对同一层次的元素，根据其对上层元素的相对重要性，采用1-9标度法进行两两比较，构建判断矩阵。例如，对于准则层中的地形地貌和岩土体性质两个元素，通过专家打分等方式，判断它们对于工程地质环境功能区划目标层的相对重要程度，若认为地形地貌比岩土体性质稍微重要，则在判断矩阵中相应位置赋值为3，反之则赋值为1/3。判断矩阵中的元素a_{ij}表示第i个元素相对于第j个元素的重要性程度，且满足a_{ij}\gt0，a_{ij}=1/a_{ji}，a_{ii}=1。层次单排序及其一致性检验：通过计算判断矩阵的最大特征根\lambda_{max}和对应的特征向量W，对特征向量进行归一化处理，得到该层次元素对于上一层次某元素的相对重要性排序权值，即层次单排序。同时，为了确保判断矩阵的一致性，需要进行一致性检验。一致性指标CI=(\lambda_{max}-n)/(n-1)，其中n为判断矩阵的阶数。引入随机一致性指标RI，根据判断矩阵的阶数从标准值表中查得。计算一致性比例CR=CI/RI，当CR\lt0.1时，认为判断矩阵具有满意的一致性，否则需要重新调整判断矩阵。例如，对于一个5阶的判断矩阵，计算得到\lambda_{max}=5.2，则CI=(5.2-5)/(5-1)=0.05，查得RI=1.12，CR=0.05/1.12\approx0.045\lt0.1，说明该判断矩阵通过一致性检验。层次总排序及其一致性检验：计算某一层次所有因素对于最高层（总目标）相对重要性的权值，即层次总排序。这一过程从最高层次到最低层次依次进行。同样需要进行一致性检验，以确保层次总排序结果的可靠性。例如，在确定了准则层各指标相对于目标层的权重，以及方案层各方案相对于准则层各指标的权重后，通过加权计算得到方案层各方案相对于目标层的总权重，同时对总排序进行一致性检验，若检验通过，则得到的总权重可用于后续的功能区划分等决策分析。层次分析法能够将复杂的决策问题层次化，充分考虑各因素之间的相互关系，使决策过程更加系统、科学。然而，该方法也存在一定的主观性，判断矩阵的构建依赖于专家的经验和判断，可能会受到人为因素的影响。3.2.2熵值法简介熵值法是一种基于信息熵理论的客观赋权方法。在信息论中，熵是对不确定性的一种度量，它可以反映信息的无序程度或者信息的效用价值。在权重确定中，熵值法通过计算各个评价指标的信息熵，来度量各个指标值的离散程度，从而确定各个指标的权重。其原理如下：数据标准化处理：由于不同评价指标的量纲和数量级可能不同，为了消除量纲的影响，需要对原始数据进行标准化处理。对于正向指标，标准化公式为x_{ij}'=\frac{x_{ij}-\min(x_j)}{\max(x_j)-\min(x_j)}；对于逆向指标，标准化公式为x_{ij}'=\frac{\max(x_j)-x_{ij}}{\max(x_j)-\min(x_j)}，其中x_{ij}为第i个评价单元的第j个指标的原始值，\max(x_j)和\min(x_j)分别为第j个指标的最大值和最小值。以地形坡度指标为例，若其为正向指标，某评价单元的地形坡度原始值为10°，该指标的最小值为0°，最大值为30°，则标准化后的值为x_{ij}'=\frac{10-0}{30-0}=\frac{1}{3}。计算指标熵值：根据标准化后的数据，计算每个指标的熵值e_j，公式为e_j=-k\sum_{i=1}^{n}p_{ij}\ln(p_{ij})，其中k=\frac{1}{\ln(n)}，p_{ij}=\frac{x_{ij}'}{\sum_{i=1}^{n}x_{ij}'}，n为评价单元的数量。熵值e_j反映了该指标值的离散程度，熵值越大，指标的离散程度越大，该指标对综合评价的影响越小。计算指标差异系数：用1减去熵值，得到指标的差异系数g_j=1-e_j。差异系数越大，说明该指标在评价中提供的有效信息量越大，对综合评价的影响越大。确定指标权重：根据差异系数的大小，确定各指标的权重w_j=\frac{g_j}{\sum_{j=1}^{m}g_j}，其中m为评价指标的数量。例如，对于地形地貌、岩土体性质、地质灾害、地下水等多个评价指标，通过上述步骤计算出各自的差异系数后，再计算出它们的权重，以反映各指标在工程地质环境功能区划评价中的相对重要性。熵值法的优点在于其客观性强，不需要事先设定权重，而是根据数据的实际情况来确定权重，能够有效避免人为因素的干扰。然而，熵值法也存在一定的局限性，它忽略了指标之间的相关性，并且对于数据的要求较高，需要数据量足够大且分布均匀。在实际应用中，通常将熵值法与层次分析法相结合，以充分发挥两种方法的优势，提高权重确定的准确性和科学性。3.3运用层次分析法对指标赋权3.3.1层次分析法步骤建立层次结构模型：将晋江市工程地质环境功能区划问题分解为目标层、准则层和指标层。目标层为晋江市工程地质环境功能区划；准则层包括地形地貌、岩土体性质、地质灾害、地下水等方面，这些准则是影响工程地质环境功能区划的主要因素。以地形地貌准则为例，其下又细分地形坡度、地形起伏度等指标层因素，地形坡度和地形起伏度直接影响着工程建设的难度和适宜性，进而影响功能区的划分。各层次之间具有明确的隶属关系，下层因素受上层因素的支配和影响，同时又对上层因素的实现起到支撑作用。构造判断矩阵：针对准则层的每个准则，邀请地质、工程等领域的专家对其下的指标层因素进行两两比较，采用1-9标度法构建判断矩阵。例如，对于地形地貌准则下的地形坡度和地形起伏度两个因素，若专家认为地形坡度对工程地质环境功能区划的影响比地形起伏度稍重要，则在判断矩阵中对应的元素赋值为3，反之则赋值为1/3。判断矩阵A=(a_{ij})_{n\timesn}中的元素a_{ij}表示第i个因素相对于第j个因素的重要性程度，满足a_{ij}\gt0，a_{ij}=1/a_{ji}，a_{ii}=1。通过这种方式，能够将专家的主观判断转化为具体的数值，为后续的权重计算提供依据。计算权重：**计算判断矩阵的最大特征根\lambda_{max}和对应的特征向量W**：利用方根法或和积法等方法计算判断矩阵的最大特征根和特征向量。以方根法为例，首先计算判断矩阵每一行元素的乘积M_i=\prod_{j=1}^{n}a_{ij}，然后计算M_i的n次方根\overline{W}_i=\sqrt[n]{M_i}，对\overline{W}_i进行归一化处理，得到特征向量W_i=\frac{\overline{W}_i}{\sum_{i=1}^{n}\overline{W}_i}，此时得到的W即为对应于最大特征根\lambda_{max}的特征向量，其元素W_i表示各指标对于上一层次某准则的相对重要性排序权值。一致性检验：计算一致性指标CI=(\lambda_{max}-n)/(n-1)，其中n为判断矩阵的阶数。引入随机一致性指标RI，根据判断矩阵的阶数从标准值表中查得。计算一致性比例CR=CI/RI，当CR\lt0.1时，认为判断矩阵具有满意的一致性，即专家的判断具有合理性和逻辑性，此时得到的权重是可靠的；若CR\geq0.1，则需要重新调整判断矩阵，直至通过一致性检验。例如，对于一个4阶的判断矩阵，计算得到\lambda_{max}=4.1，则CI=(4.1-4)/(4-1)\approx0.033，查得RI=0.90，CR=0.033/0.90\approx0.037\lt0.1，说明该判断矩阵通过一致性检验，计算得到的权重有效。层次总排序：计算指标层各指标对于目标层（工程地质环境功能区划）的总权重。将准则层各准则相对于目标层的权重与指标层各指标相对于相应准则的权重进行加权求和，得到指标层各指标对于目标层的总权重。例如，地形地貌准则相对于目标层的权重为w_1，地形坡度指标相对于地形地貌准则的权重为w_{11}，地形起伏度指标相对于地形地貌准则的权重为w_{12}，则地形坡度指标对于目标层的总权重为w_1\timesw_{11}，地形起伏度指标对于目标层的总权重为w_1\timesw_{12}。通过层次总排序，能够确定各指标在整个工程地质环境功能区划中的相对重要程度，为综合评价提供关键的权重数据。3.3.2层次分析法赋权以晋江市工程地质环境功能区划评价指标体系为例，运用层次分析法进行赋权。建立层次结构模型：构建的层次结构模型如图3-1所示，目标层为晋江市工程地质环境功能区划；准则层包括地形地貌（C_1）、岩土体性质（C_2）、地质灾害（C_3）、地下水（C_4）；指标层包括地形坡度（D_1）、地形起伏度（D_2）、岩土体类型（D_3）、岩土体抗压强度（D_4）、岩土体抗剪强度（D_5）、地质灾害危险性（D_6）、地质灾害发生频率（D_7）、地质灾害规模（D_8）、地下水水位（D_9）、地下水腐蚀性（D_{10}）。图3-1晋江市工程地质环境功能区划层次结构模型构造判断矩阵：邀请5位地质、工程领域的专家，对准则层和指标层的因素进行两两比较，采用1-9标度法构造判断矩阵。以准则层判断矩阵为例，如下表所示：||C_1|C_2|C_3|C_4||---|---|---|---|---||C_1|1|1/3|1/5|1/7||C_2|3|1|1/3|1/5||C_3|5|3|1|1/3||C_4|7|5|3|1|计算权重及一致性检验：利用方根法计算判断矩阵的最大特征根和特征向量，并进行一致性检验。对于上述准则层判断矩阵，计算过程如下：计算每一行元素的乘积M_i：M_1=1\times\frac{1}{3}\times\frac{1}{5}\times\frac{1}{7}\approx0.0095M_2=3\times1\times\frac{1}{3}\times\frac{1}{5}=0.2M_3=5\times3\times1\times\frac{1}{3}=5M_4=7\times5\times3\times1=105计算M_i的n次方根\overline{W}_i：\overline{W}_1=\sqrt[4]{0.0095}\approx0.313\overline{W}_2=\sqrt[4]{0.2}\approx0.669\overline{W}_3=\sqrt[4]{5}\approx1.495\overline{W}_4=\sqrt[4]{105}\approx3.202对\overline{W}_i进行归一化处理，得到特征向量W_i：W_1=\frac{0.313}{0.313+0.669+1.495+3.202}\approx0.057W_2=\frac{0.669}{0.313+0.669+1.495+3.202}\approx0.122W_3=\frac{1.495}{0.313+0.669+1.495+3.202}\approx0.273W_4=\frac{3.202}{0.313+0.669+1.495+3.202}\approx0.548计算最大特征根\lambda_{max}：A\timesW=\begin{pmatrix}1&\frac{1}{3}&\frac{1}{5}&\frac{1}{7}\\3&1&\frac{1}{3}&\frac{1}{5}\\5&3&1&\frac{1}{3}\\7&5&3&1\end{pmatrix}\times\begin{pmatrix}0.057\\0.122\\0.273\\0.548\end{pmatrix}=\begin{pmatrix}0.233\\0.504\\1.123\\2.274\end{pmatrix}\lambda_{max}=\frac{1}{4}\times(\frac{0.233}{0.057}+\frac{0.504}{0.122}+\frac{1.123}{0.273}+\frac{2.274}{0.548})\approx4.102计算一致性指标CI：CI=(\lambda_{max}-n)/(n-1)=(4.102-4)/(4-1)\approx0.034查随机一致性指标RI表，4阶矩阵的RI=0.90。计算一致性比例CR：CR=CI/RI=0.034/0.90\approx0.038\lt0.1，通过一致性检验。同理，对指标层各判断矩阵进行计算和一致性检验，得到各指标相对于准则层的权重。例如，地形地貌准则下地形坡度和地形起伏度的判断矩阵及权重计算结果如下：D_1D_2D_113D_21/31计算得到地形坡度（D_1）的权重为0.75，地形起伏度（D_2）的权重为0.25，且通过一致性检验。层次总排序：计算指标层各指标对于目标层的总权重，结果如下表所示：|准则层|权重|指标层|相对于准则层权重|总权重||---|---|---|---|---||C_1|0.057|D_1|0.75|0.057\times0.75=0.043||||D_2|0.25|0.057\times0.25=0.014||C_2|0.122|D_3|0.4|0.122\times0.4=0.049||||D_4|0.3|0.122\times0.3=0.037||||D_5|0.3|0.122\times0.3=0.037||C_3|0.273|D_6|0.5|0.273\times0.5=0.137||||D_7|0.3|0.273\times0.3=0.082||||D_8|0.2|0.273\times0.2=0.055||C_4|0.548|D_9|0.6|0.548\times0.6=0.329||||D_{10}|0.4|0.548\times0.4=0.219|通过上述层次分析法赋权过程，确定了晋江市工程地质环境功能区划评价指标体系中各指标的权重，其中地下水水位（D_9）的总权重最大，为0.329，说明在晋江市工程地质环境功能区划中，地下水水位对功能区划分的影响最为重要；其次是地质灾害危险性（D_6），总权重为0.137；而地形起伏度（D_2）的总权重相对较小，为0.014。这些权重结果将为后续的工程地质环境适宜性综合评价提供重要依据，使评价结果更加科学、合理。3.4利用熵值法对指标赋权3.4.1熵值法步骤数据标准化处理：由于晋江市工程地质环境功能区划评价指标体系中的各指标量纲和数量级不同，为了消除量纲的影响，需要对原始数据进行标准化处理。对于正向指标，即指标值越大对工程地质环境功能区划越有利的指标，如岩土体抗压强度，标准化公式为x_{ij}'=\frac{x_{ij}-\min(x_j)}{\max(x_j)-\min(x_j)}；对于逆向指标，即指标值越小对工程地质环境功能区划越有利的指标，如地质灾害危险性，标准化公式为x_{ij}'=\frac{\max(x_j)-x_{ij}}{\max(x_j)-\min(x_j)}。其中x_{ij}为第i个评价单元的第j个指标的原始值，\max(x_j)和\min(x_j)分别为第j个指标的最大值和最小值。以地形坡度指标为例，若某评价单元的地形坡度原始值为15°，该指标在所有评价单元中的最小值为0°，最大值为30°，由于地形坡度为逆向指标，根据标准化公式可得x_{ij}'=\frac{30-15}{30-0}=0.5。计算指标熵值：根据标准化后的数据，计算每个指标的熵值e_j，公式为e_j=-k\sum_{i=1}^{n}p_{ij}\ln(p_{ij})，其中k=\frac{1}{\ln(n)}，p_{ij}=\frac{x_{ij}'}{\sum_{i=1}^{n}x_{ij}'}，n为评价单元的数量。熵值e_j反映了该指标值的离散程度，熵值越大，指标的离散程度越大，该指标对综合评价的影响越小。例如，对于岩土体类型指标，假设有n=100个评价单元，经过标准化处理后，计算出p_{ij}的值，再代入熵值计算公式，可得到该指标的熵值。计算指标差异系数：用1减去熵值，得到指标的差异系数g_j=1-e_j。差异系数越大，说明该指标在评价中提供的有效信息量越大，对综合评价的影响越大。例如，若某指标的熵值为0.8，则其差异系数为g_j=1-0.8=0.2。确定指标权重：根据差异系数的大小，确定各指标的权重w_j=\frac{g_j}{\sum_{j=1}^{m}g_j}，其中m为评价指标的数量。通过该公式，可计算出每个指标在工程地质环境功能区划评价中的相对重要性权重。例如，对于包含地形地貌、岩土体性质、地质灾害、地下水等多个评价指标的体系，计算出各指标的差异系数后，再计算出它们的权重，以确定各指标在功能区划中的权重分配。3.4.2熵值法赋权收集晋江市工程地质环境相关数据，包括地形坡度、地形起伏度、岩土体类型、岩土体抗压强度、岩土体抗剪强度、地质灾害危险性、地质灾害发生频率、地质灾害规模、地下水水位、地下水腐蚀性等指标的原始数据，共涉及n=500个评价单元。首先进行数据标准化处理，以岩土体抗压强度（正向指标）为例，该指标原始数据中最小值为50MPa，最大值为200MPa，某评价单元的原始值为100MPa，则标准化后的值为x_{ij}'=\frac{100-50}{200-50}\approx0.33。对于地质灾害危险性（逆向指标），原始数据中最小值为1（低危险等级），最大值为5（高危险等级），某评价单元的原始值为3，标准化后的值为x_{ij}'=\frac{5-3}{5-1}=0.5。接着计算指标熵值，对于地形坡度指标，计算p_{ij}后，k=\frac{1}{\ln(500)}，经计算得到其熵值e_{坡度}\approx0.92。同理计算其他指标熵值，如岩土体抗压强度熵值e_{抗压}\approx0.85。然后计算指标差异系数，地形坡度的差异系数g_{坡度}=1-0.92=0.08，岩土体抗压强度的差异系数g_{抗压}=1-0.85=0.15。最后确定指标权重，假设计算得到所有m=10个指标的差异系数总和\sum_{j=1}^{10}g_j=1.2，则地形坡度的权重w_{坡度}=\frac{0.08}{1.2}\approx0.067，岩土体抗压强度的权重w_{抗压}=\frac{0.15}{1.2}=0.125。以此类推，计算出其他指标的权重，结果如下表所示：评价指标熵值法权重地形坡度0.067地形起伏度0.055岩土体类型0.112岩土体抗压强度0.125岩土体抗剪强度0.120地质灾害危险性0.135地质灾害发生频率0.095地质灾害规模0.088地下水水位0.140地下水腐蚀性0.163通过熵值法赋权，确定了各评价指标在晋江市工程地质环境功能区划中的相对重要性权重。其中，地下水腐蚀性的权重最大，为0.163，说明其在功能区划中对评价结果的影响相对较大；地形起伏度的权重相对较小，为0.055。这些权重结果将为后续的工程地质环境适宜性综合评价提供客观的数据支持，使得评价结果更能反映晋江市工程地质环境的实际情况。3.5组合权重3.5.1组合权重原理层次分析法（AHP）是一种定性与定量相结合的多目标决策方法，它通过构建层次结构模型，将复杂的决策问题分解为多个层次，通过专家的主观判断对各层次元素进行两两比较，从而确定各指标的相对重要性权重。这种方法能够充分考虑决策者的经验和偏好，体现各指标的主观重要性。然而，由于其依赖专家判断，主观性较强，不同专家的判断可能存在差异，导致权重结果不够客观。熵值法是一种基于信息熵理论的客观赋权方法。它通过计算各评价指标的信息熵来度量指标值的离散程度，从而确定指标权重。信息熵越大，说明该指标的数据离散程度越小，提供的有效信息量越少，其权重也就越小；反之，信息熵越小，指标权重越大。熵值法完全依据数据本身的特征来确定权重，不受主观因素的影响，具有较强的客观性。但熵值法也存在一定局限性，它只考虑了数据的离散程度，忽略了指标之间的相关性，且对数据量和数据分布有一定要求。为了充分发挥两种方法的优势，弥补各自的不足，本研究采用组合权重法，将层次分析法得到的主观权重与熵值法得到的客观权重进行有机结合。组合权重法的基本原理是通过一定的数学方法，将两种权重进行加权求和，得到综合考虑主观和客观因素的组合权重。其计算公式为：W_{组合}=αW_{AHP}+(1-α)W_{熵值法}其中，W_{组合}为组合权重，W_{AHP}为层次分析法确定的权重，W_{熵值法}为熵值法确定的权重，α为权重系数，取值范围为[0,1]，它反映了主观权重和客观权重在组合权重中的相对重要程度。α的确定通常需要综合考虑研究目的、数据特点以及决策者对主观和客观因素的重视程度等因素。一般可以通过专家打分、经验判断或采用优化算法（如遗传算法、粒子群优化算法等）来确定α的值，以使得组合权重能够更准确地反映各指标在工程地质环境功能区划中的实际重要性。3.5.2组合权重计算在确定α值时，本研究邀请了5位在工程地质领域具有丰富经验的专家进行打分。专家们根据晋江市工程地质环境的特点、数据的可靠性以及评价的目的，对主观权重和客观权重的相对重要程度进行评价，打分范围为0-1。经过统计分析，专家们给出的α值平均值为0.4，即认为主观权重在组合权重中占40\%的比重，客观权重占60\%的比重。根据组合权重公式W_{组合}=αW_{AHP}+(1-α)W_{熵值法}，将层次分析法得到的权重（W_{AHP}）和熵值法得到