余悠然 王冉 李振 劉世奇 王鶴

摘要：黃土高原地區(qū)地表破碎，水土流失嚴(yán)重，易發(fā)生滑坡、崩塌、地面塌陷等地質(zhì)災(zāi)害，損毀農(nóng)田，威脅農(nóng)業(yè)生產(chǎn)安全。以山西省古交市為研究區(qū)，采用信息量耦合二元邏輯（Logistic）回歸模型和地理信息系統(tǒng)（GIS）空間分析工具，在劃定地質(zhì)災(zāi)害易發(fā)性分區(qū)的基礎(chǔ)上，開(kāi)展古交市基本農(nóng)田地質(zhì)環(huán)境安全性評(píng)價(jià)研究，研究結(jié)果表明：（1）古交市地質(zhì)災(zāi)害極高、高、中、低易發(fā)區(qū)的面積占比分別為9.70%、12.27%、21.17%、56.86%，中心城區(qū)、城鎮(zhèn)村居民區(qū)和交通干線附近地質(zhì)災(zāi)害發(fā)生概率最高;（2）古交市基本農(nóng)田地質(zhì)環(huán)境安全性分為極低、低、中、高4個(gè)等級(jí)，其對(duì)應(yīng)范圍內(nèi)基本農(nóng)田面積占比分別為11.82%、17.02%、26.67%、44.49%。古交市地質(zhì)環(huán)境已對(duì)基本農(nóng)田安全造成威脅，黃土高原地區(qū)基本農(nóng)田劃定和布局優(yōu)化調(diào)整時(shí)，應(yīng)深入考慮地質(zhì)災(zāi)害、水土流失等地質(zhì)環(huán)境問(wèn)題造成的影響。

關(guān)鍵詞：黃土高原;古交市;地質(zhì)災(zāi)害;基本農(nóng)田;地質(zhì)環(huán)境安全性評(píng)價(jià)

中圖分類(lèi)號(hào)：X43 文獻(xiàn)標(biāo)志碼： A 文章編號(hào)：1002-1302（2021）16-0220-07

我國(guó)耕地資源稀缺，嚴(yán)守耕地紅線，保證糧食安全是我國(guó)的基本國(guó)策，為此我國(guó)實(shí)行基本農(nóng)田保護(hù)制度[1]，而基本農(nóng)田保護(hù)的核心內(nèi)容就是要保存基本農(nóng)田的生產(chǎn)力[2]。黃土垂直節(jié)理發(fā)育，具有較強(qiáng)的濕陷性，因此黃土高原是我國(guó)地質(zhì)災(zāi)害的頻發(fā)地區(qū)[3]。各類(lèi)地質(zhì)災(zāi)害不僅會(huì)破壞農(nóng)業(yè)基礎(chǔ)設(shè)施，還會(huì)導(dǎo)致土壤肥力下降、耕作層受到破壞等后果，進(jìn)而造成農(nóng)田生產(chǎn)力下降。目前，我國(guó)基本農(nóng)田的劃定多是從自然稟賦、區(qū)位條件[4]等方面考慮，我國(guó)學(xué)者對(duì)基本農(nóng)田質(zhì)量的評(píng)價(jià)和布局的優(yōu)化也主要從自然氣候、地力大小、基礎(chǔ)設(shè)施、經(jīng)濟(jì)效益、區(qū)位因素、環(huán)境清潔程度和可持續(xù)利用幾個(gè)方面進(jìn)行[5-8]，而對(duì)地質(zhì)災(zāi)害、水土流失等可能造成農(nóng)田損毀、產(chǎn)量減少的地質(zhì)環(huán)境影響考慮不足[9]。受滑坡、崩塌、泥石流等地質(zhì)災(zāi)害破壞的農(nóng)田，短期內(nèi)難以恢復(fù)耕種，造成耕地面積減少。山西省古交市位于黃土高原東部，地表溝壑縱橫，多黃土塬、梁、峁，耕地面積少，是我國(guó)重要的產(chǎn)煤基地和最大的煉焦基地，然而長(zhǎng)期以來(lái)各種人類(lèi)活動(dòng)對(duì)當(dāng)?shù)氐乇碓斐闪藝?yán)重破壞，地質(zhì)災(zāi)害頻繁發(fā)生，給當(dāng)?shù)鼗巨r(nóng)田保有量和糧食安全產(chǎn)生了威脅。目前評(píng)價(jià)地質(zhì)災(zāi)害易發(fā)性的方法有信息量法、邏輯回歸法、隨機(jī)森林算法等，綜合考慮研究區(qū)實(shí)際情況以及資料的豐富程度。本研究采用信息量耦合二元邏輯（Logistic）回歸方法分析山西省古交市地質(zhì)災(zāi)害易發(fā)性，并在此基礎(chǔ)上開(kāi)展基本農(nóng)田地質(zhì)環(huán)境安全性評(píng)價(jià)。目前針對(duì)研究區(qū)基本農(nóng)田展開(kāi)的研究較少，當(dāng)?shù)貙?duì)基本農(nóng)田受地質(zhì)災(zāi)害威脅的重視不足，本研究將為古交市基本農(nóng)田保護(hù)工作提供參考，也將為其他地質(zhì)災(zāi)害易發(fā)區(qū)的基本農(nóng)田布局調(diào)整優(yōu)化提供科學(xué)依據(jù)，幫助減少因地質(zhì)災(zāi)害而造成的基本農(nóng)田折損。

1 研究區(qū)域及數(shù)據(jù)

1.1 研究區(qū)概況

古交市地處黃土高原東部、山西省中部（111°43′08″～112°21′05″E，37°40′06″～38°08′09″N）（圖1），隸屬太原市下轄縣級(jí)市，總面積為1 512.345 km2，現(xiàn)轄3鎮(zhèn)7鄉(xiāng)4個(gè)街道，總?cè)丝跒?1.18萬(wàn)人，國(guó)內(nèi)生產(chǎn)總值（GDP）為36.61億元，糧食總產(chǎn)量為 10 955 t（2018年）。2019年，全市共有耕地 19 740.26 hm2，基本農(nóng)田15 792.33 hm2，主要分布于古交市北部的河口鎮(zhèn)和南部的岔口鄉(xiāng)、邢家社鄉(xiāng)，占基本農(nóng)田總數(shù)的51.71%。古交市屬于溫帶大陸性氣候，年均降水量為391.58 mm，四面環(huán)山，中部地勢(shì)低緩，平均海拔為1 604 m，形成了以中山為主的地貌形態(tài)。地表黃土厚度小于50 m，多分布于低山丘陵區(qū)和中部梁峁及其邊坡地帶，區(qū)域內(nèi)河谷縱橫，水土流失嚴(yán)重。境內(nèi)礦產(chǎn)資源豐富，多年來(lái)工礦開(kāi)采給當(dāng)?shù)卦斐芍脖黄茐?、地面塌陷等?yán)重環(huán)境問(wèn)題。據(jù)《山西省2019年度地質(zhì)災(zāi)害防治方案》，山西省全境均屬地質(zhì)災(zāi)害易發(fā)區(qū)。古交市地質(zhì)災(zāi)害頻繁發(fā)生，截至2017年底，共有287個(gè)地質(zhì)災(zāi)害（隱患）點(diǎn)，其中已發(fā)生的各類(lèi)地質(zhì)災(zāi)害點(diǎn)有117個(gè)，各類(lèi)潛在地質(zhì)災(zāi)害隱患點(diǎn)有170個(gè)。

1.2 數(shù)據(jù)來(lái)源

數(shù)字高程模型（DEM）數(shù)據(jù)來(lái)源于地理空間數(shù)據(jù)云平臺(tái)的ASTER GDEM 30 m分辨率數(shù)字高程數(shù)據(jù)，在地理信息系統(tǒng)（GIS）軟件中提取坡度、高程、坡向等信息。地質(zhì)災(zāi)害（隱患）點(diǎn)來(lái)源于《山西省古交市地質(zhì)災(zāi)害詳細(xì)調(diào)查報(bào)告》?；巨r(nóng)田、道路、河流、城鎮(zhèn)村居民區(qū)等數(shù)據(jù)來(lái)源于古交市土地利用現(xiàn)狀圖，本研究結(jié)合遙感影像對(duì)其進(jìn)行了修正。地貌、地質(zhì)數(shù)據(jù)來(lái)源于BIGMAP軟件，用于提取地貌類(lèi)型、地質(zhì)構(gòu)造、巖土體類(lèi)型信息，上述數(shù)據(jù)均通過(guò)GIS軟件進(jìn)行配準(zhǔn)、投影變換、矢量化等處理，坐標(biāo)系統(tǒng)一為2000國(guó)家大地坐標(biāo)系（China Geodetic Coordinate System 2000，CGCS2000）。

2 研究方法

2.1 地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)方法

基本農(nóng)田劃定會(huì)對(duì)耕地土壤清潔程度、田面坡度、集中連片度、水利和水土保持設(shè)施[10-12]等提出具體要求，但對(duì)地質(zhì)災(zāi)害的重視程度尚有不足。地表坡度、水土保持狀況等在一定程度上可以表征地質(zhì)環(huán)境對(duì)基本農(nóng)田的約束作用，但倘若周邊發(fā)生地質(zhì)災(zāi)害，依然會(huì)對(duì)農(nóng)田造成負(fù)面影響。本研究采用構(gòu)建信息量法和二元Logistic回歸方法的耦合模型，開(kāi)展古交市地質(zhì)災(zāi)害易發(fā)性研究，同單一模型相比，耦合模型可以在不過(guò)多增加工作量的前提下顯著提高模型預(yù)測(cè)精度[13]。在此基礎(chǔ)上，進(jìn)一步分析地質(zhì)災(zāi)害易發(fā)性對(duì)基本農(nóng)田安全的影響。本研究以287個(gè)地質(zhì)災(zāi)害（隱患）點(diǎn)和287個(gè)隨機(jī)生成的非地質(zhì)災(zāi)害點(diǎn)作為樣本，將地質(zhì)災(zāi)害發(fā)生定義為1，不發(fā)生定義為0，計(jì)算各分級(jí)影響因子的信息量值，并通過(guò)相關(guān)性分析和共線性診斷剔除相關(guān)性過(guò)強(qiáng)的多余因子，提高模型預(yù)測(cè)精度，再利用二元Logistic回歸模型得到最終的地質(zhì)災(zāi)害易發(fā)性預(yù)測(cè)結(jié)果。信息量法和二元Logistic回歸模型表達(dá)式如下：

式中：A表示地質(zhì)災(zāi)害發(fā)生;I（Xi，A）表示在A條件下致災(zāi)因子（Xi）提供的信息量值;Ni為分布在因素Xi內(nèi)變形破壞單元總數(shù);N為研究區(qū)已知變形破壞單元總數(shù);Si為研究區(qū)內(nèi)含有評(píng)價(jià)因素Xi的面積;S為研究區(qū)總面積;β1～βk為影響系數(shù);α為常量;Pi是地質(zhì)災(zāi)害發(fā)生的概率;1-Pi為災(zāi)害不發(fā)生的概率。

2.2 影響因子分析

2.2.1 地質(zhì)災(zāi)害影響因子分析 根據(jù)研究區(qū)現(xiàn)場(chǎng)調(diào)查情況和地質(zhì)特點(diǎn)，本研究將從自然因素和人為因素論述各要素對(duì)地質(zhì)災(zāi)害發(fā)育的影響。自然因素包括地貌類(lèi)型、高程、坡度、坡向、巖土體類(lèi)型、斷層、河流等，人為因素包括城鎮(zhèn)村居民區(qū)建設(shè)、道路建設(shè)和采礦活動(dòng)等。

從地貌類(lèi)型來(lái)看，古交市梁狀黃土丘陵區(qū)的地質(zhì)災(zāi)害最為集中，該區(qū)域黃土垂直節(jié)理發(fā)育、土質(zhì)松散，是地質(zhì)災(zāi)害形成的有利條件，同時(shí)又與采礦等人類(lèi)活動(dòng)頻繁區(qū)有較多重疊，更易誘發(fā)地質(zhì)災(zāi)害。單從坡度上來(lái)看，坡度越大的區(qū)域地質(zhì)災(zāi)害發(fā)生的可能性越大，但高程高、坡度大的區(qū)域多屬于建設(shè)用地不適宜區(qū)，人類(lèi)活動(dòng)極少，誘發(fā)地質(zhì)災(zāi)害的可能性較低。從坡向上來(lái)看，陽(yáng)坡受太陽(yáng)輻照時(shí)間更長(zhǎng)，巖石風(fēng)化程度更高，更有利于地質(zhì)災(zāi)害發(fā)育。從巖土體類(lèi)型來(lái)看，粉土巖組分布的區(qū)域地質(zhì)災(zāi)害數(shù)量最多，粉土垂直節(jié)理發(fā)育且具有濕陷性，當(dāng)大量雨水沖刷時(shí)，極易形成泥石流和滑坡等地質(zhì)災(zāi)害;斷層是巖層所受應(yīng)力超過(guò)自身極限而產(chǎn)生的斷裂錯(cuò)位，在距離斷層較近的地區(qū)，巖石破碎易被風(fēng)化，巖土結(jié)構(gòu)更為脆弱，也更易發(fā)生地質(zhì)災(zāi)害;河流的浸泡、侵蝕對(duì)河谷兩岸坡體的穩(wěn)定性和強(qiáng)度也會(huì)有影響。古交市全境均為地質(zhì)災(zāi)害易發(fā)區(qū)，各項(xiàng)工程建設(shè)過(guò)程中可能出現(xiàn)不合理開(kāi)挖坡腳、隨意堆放土石土方等問(wèn)題，易誘發(fā)滑坡、崩塌等地質(zhì)災(zāi)害;跟道路建設(shè)和居民區(qū)建設(shè)相比，采礦活動(dòng)對(duì)地表及周邊環(huán)境的改造最為強(qiáng)烈，生產(chǎn)過(guò)程中產(chǎn)生的采空區(qū)和廢渣的隨意堆積，為各類(lèi)地質(zhì)災(zāi)害的發(fā)生創(chuàng)造了條件。

2.2.2 影響因子分級(jí) 為便于計(jì)算影響因子的信息量值，本研究對(duì)10個(gè)影響因子進(jìn)行分級(jí)。對(duì)于離散型數(shù)據(jù)，直接依據(jù)野外實(shí)地調(diào)查時(shí)的分級(jí);對(duì)于連續(xù)型數(shù)據(jù)，依據(jù)各分級(jí)區(qū)間內(nèi)災(zāi)害點(diǎn)比例與點(diǎn)密度曲線的突變點(diǎn)分級(jí)[14]。首先將影響因子按步長(zhǎng)進(jìn)行分割，即高程、距道路距離、距斷層距離按 200 m，坡度按5°，距河流距離按100 m為步長(zhǎng)進(jìn)行災(zāi)害點(diǎn)比例與點(diǎn)密度統(tǒng)計(jì)。

如圖2所示，將高程按1 000、1 200、1 600分級(jí);坡度按10°、25°、45°分級(jí);距斷層距離按400、800、1 800 m分級(jí);距河流距離按200、500、800 m分級(jí);距道路距離按200、800、1 200 m分級(jí)。

2.2.3 信息量值計(jì)算 根據(jù)影響因子分級(jí)結(jié)果和公式（1）計(jì)算各級(jí)影響因子的信息量值，如表1所示。

2.2.4 影響因子相關(guān)性分析 為保證影響因子間的相互獨(dú)立，本研究在SPSS中采用斯皮爾曼方法進(jìn)行因子間的相關(guān)性分析，并結(jié)合研究區(qū)實(shí)際情況剔除相關(guān)性較大的因子。當(dāng)相關(guān)系數(shù)大于0.5時(shí)，認(rèn)為2個(gè)因子間存在一定的相關(guān)性[15]。如表2所示，地貌類(lèi)型、高程、距斷層距離和距道路距離與其他因子間的相關(guān)系數(shù)有的超過(guò)0.5，考慮到斷層對(duì)區(qū)域地質(zhì)穩(wěn)定性的影響，且研究區(qū)交通干線兩側(cè)存在大量無(wú)防護(hù)的高陡邊坡，這2個(gè)因子對(duì)地質(zhì)災(zāi)害的發(fā)生有較重要的影響，應(yīng)予以保留，最終選取巖土體類(lèi)型、坡度、坡向、距斷層距離、距道路距離、距水系距離、距居民區(qū)距離、距礦區(qū)距離為地質(zhì)災(zāi)害易發(fā)性模型評(píng)價(jià)因子。

2.3 地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)模型計(jì)算

在進(jìn)行二元Logistic回歸分析前，須要對(duì)影響因子進(jìn)行共線性診斷，各要素的方差膨脹因子（VIF）值均小于10，表明因子間不存在多重共線性，可以進(jìn)行回歸分析。導(dǎo)入287個(gè)地質(zhì)災(zāi)害點(diǎn)和287個(gè)非

地質(zhì)災(zāi)害點(diǎn)作為樣本，得到Logistic回歸方程：

lgP=0.468+0.559x1+0.483x2+0.761x3+0.535x4-0.118x5+0.999x6+1.070x7+0.873x8。

式中：x1、x2、…、 x8分別對(duì)應(yīng)巖土體類(lèi)型、坡度、坡向、距斷層距離、距水系距離、距道路距離、距居民區(qū)距離和距礦區(qū)距離。從回歸系數(shù)大小來(lái)看，距道路距離、距居民區(qū)距離和距礦區(qū)距離的系數(shù)分別為0.999、1.070、0.873，表明這3個(gè)因子對(duì)地質(zhì)災(zāi)害發(fā)生起主要控制作用，說(shuō)明研究區(qū)內(nèi)地質(zhì)災(zāi)害發(fā)生的可能性與人類(lèi)活動(dòng)的頻繁程度有較大相關(guān)性。而巖土體類(lèi)型、坡度、坡向等影響因子不是古交市地質(zhì)災(zāi)害發(fā)生的最主要影響因素。

為檢驗(yàn)?zāi)Ｐ途?，本研究選用接受者操作特征（ROC）曲線對(duì)模型預(yù)測(cè)結(jié)果進(jìn)行驗(yàn)證，當(dāng)曲線下面積（AUC）為0.5～1.0時(shí)，說(shuō)明結(jié)果可靠[15-16]。將模型預(yù)測(cè)結(jié)果和樣本數(shù)據(jù)檢驗(yàn)值代入得到AUC為0.897 8，遠(yuǎn)大于0.5，證明該模型計(jì)算結(jié)果可靠。

3 結(jié)果與分析

3.1 地質(zhì)災(zāi)害易發(fā)性分區(qū)

將得到的各項(xiàng)回歸系數(shù)代入公式（2）和公式（3），在GIS軟件中通過(guò)柵格計(jì)算工具將各影響因子圖層疊加得到每個(gè)柵格單元地質(zhì)災(zāi)害發(fā)生概率（P），其范圍為0.009 94～0.999 48。按自然斷點(diǎn)法分為4個(gè)區(qū)間：（0.009 94，0.196 21）（低易發(fā)區(qū)）、[0.0196 21，0.479 49）（中易發(fā)區(qū)）、[0.479 49，0.782 17）（高易發(fā)區(qū)）、[0.782 17，0.999 48]（極高易發(fā)區(qū)），各分區(qū)面積占比分別為56.86%、21.17%、12.27%、9.70%（圖3）。極高和高易發(fā)區(qū)主要集中于主城區(qū)、居民區(qū)和交通干道沿線，中易發(fā)區(qū)多與礦區(qū)分布相重合，低易發(fā)區(qū)分布在古交市四周人口密度低、人類(lèi)活動(dòng)較少的山地地區(qū)。

3.2 基本農(nóng)田地質(zhì)環(huán)境安全性評(píng)價(jià)

依據(jù)前文得到的古交市地質(zhì)災(zāi)害易發(fā)性分區(qū)，對(duì)古交市現(xiàn)有基本農(nóng)田地質(zhì)環(huán)境安全性進(jìn)行評(píng)價(jià)。將地質(zhì)災(zāi)害易發(fā)性等級(jí)分區(qū)圖和基本農(nóng)田分布圖進(jìn)行疊加分析，地質(zhì)災(zāi)害易發(fā)性等級(jí)和基本農(nóng)田安全性等級(jí)（圖4）的對(duì)應(yīng)關(guān)系為極高-極低、高-低、中-中、低-高。從古交市基本農(nóng)田地質(zhì)環(huán)境安全性評(píng)價(jià)等級(jí)圖（圖4）和安全性分級(jí)統(tǒng)計(jì)表（表3）來(lái)看，古交市中部地區(qū)的基本農(nóng)田地質(zhì)環(huán)境安全性最低，并沿主要交通干線呈射線狀向四周輻射。位于古交市中部的鎮(zhèn)城底鎮(zhèn)、梭峪鄉(xiāng)、西曲街辦、屯蘭街辦、桃園街辦和東曲街辦，其范圍內(nèi)有22.43%的基本農(nóng)田安全性等級(jí)為極低，30.91%的基本農(nóng)田安全性等級(jí)為低。在不采取任何措施的情況下，古交市有近30%的基本農(nóng)田將受到地質(zhì)災(zāi)害的破壞。

將基本農(nóng)田地質(zhì)環(huán)境安全性評(píng)價(jià)結(jié)果與道路、城鎮(zhèn)村居民區(qū)、礦區(qū)疊加展示，如圖5所示，可發(fā)現(xiàn)道路兩側(cè)和居民區(qū)周邊的基本農(nóng)田地質(zhì)環(huán)境安全性等級(jí)多為極低和低，經(jīng)統(tǒng)計(jì)，在道路和居民區(qū)周邊100 m范圍內(nèi)，有基本農(nóng)田397.35 hm2，其中地質(zhì)環(huán)境安全性等級(jí)為極低的基本農(nóng)田有258.20 hm2，占64.98%，安全性等級(jí)為低的占30.68%;礦區(qū)范圍內(nèi)的基本農(nóng)田共有3 247.43 hm2，安全性等級(jí)為極低的基本農(nóng)田占24.40%，安全性等級(jí)為低的基本農(nóng)田占29.51%。根據(jù)表3可知，有90.10%的極低安全等級(jí)、67.96%低安全等級(jí)基本農(nóng)田是分布在人類(lèi)工程活動(dòng)影響范圍內(nèi)的。經(jīng)實(shí)地考查發(fā)現(xiàn)，古交城市中心和交通干線附近，仍有部分工程活動(dòng)遺留的高陡邊坡，且缺少必要的防護(hù)措施。由于古交市多黃土梁峁分布，當(dāng)?shù)卮迕襁x址建房時(shí)難以遠(yuǎn)離邊坡，施工時(shí)若對(duì)周?chē)履_及巖土體結(jié)構(gòu)造成破壞，易誘發(fā)地質(zhì)災(zāi)害，增加對(duì)周邊基本農(nóng)田的威脅。長(zhǎng)期的礦產(chǎn)資源開(kāi)發(fā)更是對(duì)地表的強(qiáng)烈改造，生產(chǎn)過(guò)程中不僅會(huì)對(duì)土壤造成污染，還會(huì)引發(fā)地裂縫、地面塌陷等問(wèn)題，形成規(guī)模較大的塌坑和積水洼地，加劇水土流失速率，導(dǎo)致土壤肥力下降，使基本農(nóng)田質(zhì)量嚴(yán)重下降甚至無(wú)法耕種。而分布在研究區(qū)四周山地地區(qū)的基本農(nóng)田安全性等級(jí)較高，這是由于該區(qū)域高程相對(duì)較高，又無(wú)大量礦產(chǎn)資源，人口密度小，人類(lèi)工程活動(dòng)對(duì)地表的擾動(dòng)少，地表植被相對(duì)茂密，可鞏固水土、阻滯降雨對(duì)地表的沖刷，從而在一定程度上提高了該區(qū)域基本農(nóng)田的安全性等級(jí)。

4 結(jié)論與展望

本研究利用信息量耦合二元Logistic回歸模型對(duì)古交市基本農(nóng)田地質(zhì)環(huán)境安全性進(jìn)行了評(píng)價(jià)，得到以下結(jié)論：（1）古交市地質(zhì)災(zāi)害易發(fā)性等級(jí)分為極高易發(fā)區(qū)、高易發(fā)區(qū)、中易發(fā)區(qū)和低易發(fā)區(qū)，面積占比分別為9.70%、12.27%、21.17%、56.86%，易受地質(zhì)災(zāi)害影響的范圍大。極高和高易發(fā)區(qū)主要集中于主城區(qū)、居民區(qū)和交通干道沿線，中易發(fā)區(qū)多與礦區(qū)分布相重合，低易發(fā)區(qū)分布在古交市四周人口密度低、人類(lèi)活動(dòng)較少的山地地區(qū)。（2）以地質(zhì)災(zāi)害易發(fā)性為基礎(chǔ)開(kāi)展古交市基本農(nóng)田地質(zhì)環(huán)境安全性評(píng)價(jià)研究，基本農(nóng)田安全性等級(jí)為極低、低、中、高所對(duì)應(yīng)的面積占比分別為11.82%、17.02%、26.67%、44.49%，古交市基本農(nóng)田安全正面臨地質(zhì)環(huán)境的嚴(yán)重威脅，各項(xiàng)人類(lèi)活動(dòng)是導(dǎo)致地質(zhì)環(huán)境對(duì)基本農(nóng)田產(chǎn)生負(fù)面影響的主要因素，積極采取相應(yīng)措施保護(hù)基本農(nóng)田數(shù)量、質(zhì)量不受破壞已刻不容緩。以古交市為參考，黃土高原地區(qū)將城鎮(zhèn)周邊、交通沿線易被占用的優(yōu)質(zhì)耕地劃為基本農(nóng)田或開(kāi)展基本農(nóng)田布局優(yōu)化時(shí)，應(yīng)深入考慮當(dāng)?shù)氐刭|(zhì)環(huán)境對(duì)基本農(nóng)田的約束作用，切實(shí)做好水土保持、防災(zāi)減災(zāi)等相應(yīng)工作，必要時(shí)須對(duì)基本農(nóng)田周邊地質(zhì)災(zāi)害隱患點(diǎn)進(jìn)行工程治理，對(duì)于現(xiàn)有地質(zhì)災(zāi)害易發(fā)區(qū)內(nèi)的基本農(nóng)田，可以通過(guò)土地置換的方式進(jìn)行調(diào)整，確保做到牢守耕地紅線，保護(hù)糧食安全。

參考文獻(xiàn)：

[1]周尚意，朱阿興，邱維理，等. 基于GIS的農(nóng)用地連片性分析及其在基本農(nóng)田保護(hù)規(guī)劃中的應(yīng)用[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2008，24（7）：72-77，313.

[2]奉 婷，張鳳榮，李 燦，等. 基于耕地質(zhì)量綜合評(píng)價(jià)的縣域基本農(nóng)田空間布局[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2014，30（1）：200-210，293.

[3]孫萍萍，張茂省，程秀娟，等. 黃土高原地質(zhì)災(zāi)害發(fā)生規(guī)律[J]. 山地學(xué)報(bào)，2019，37（5）：737-746.

[4]董光龍，趙 軒，劉金花，等. 基于耕地質(zhì)量評(píng)價(jià)與空間集聚特征的基本農(nóng)田劃定研究[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào)，2020，51（2）：133-142.

[5]Song Y F，Lu Y J，Liu T J，et.al. Variation of vegetation fractional coverage and its relationship with climate in a desert steppe：optimization of farmland layout in a farming-pastoral ecotone using the ecological suitability index[J]. Ecological Engineering，2020，150：105834.

[6]付國(guó)珍，擺萬(wàn)奇.耕地質(zhì)量評(píng)價(jià)研究進(jìn)展及發(fā)展趨勢(shì)[J]. 資源科學(xué)，2015，37（2）：226-236.

[7]Cheng Q W，Jiang P H，Cai L Y，et al. Delineation of a permanent basic farmland protection area around a city centre：case study of Changzhou City，China[J]. Land Use Policy，2017，60：73-89.

[8]劉霈珈，吳克寧，趙華甫，等. 基于耕地綜合質(zhì)量的基本農(nóng)田布局優(yōu)化——以河南省溫縣為例[J]. 中國(guó)土地科學(xué)，2015，29（2）：54-59，2.

[9]王玉軍，劉 瓊，歐名豪. 基于地質(zhì)環(huán)境適宜性的基本農(nóng)田布局優(yōu)化[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2019，35（14）：252-260.

[10]曹麗萍，羅志軍，冉鳳維，等. 基于耕地質(zhì)量和空間集聚格局的縣域基本農(nóng)田劃定[J]. 水土保持研究，2018，25（4）：349-355，364.

[11]周 明，王占岐.基于耕地質(zhì)量和空間聚類(lèi)的縣域基本農(nóng)田劃定——以湖北省團(tuán)風(fēng)縣為例[J]. 水土保持研究，2016，23（5）：316-321.

[12]吳彥澎. 一種基于網(wǎng)格的空間聚類(lèi)算法及在基本農(nóng)田規(guī)劃中的應(yīng)用[J]. 江蘇農(nóng)業(yè)科學(xué)，2020，48（16）：250-254.

[13]樊芷吟，茍曉峰，秦明月，等. 基于信息量模型與Logistic回歸模型耦合的地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)[J]. 工程地質(zhì)學(xué)報(bào)，2018，26（2）：340-347.

[14]廖蕓婧. 安徽省歙縣地質(zhì)災(zāi)害易發(fā)性區(qū)劃與氣象預(yù)警研究[D]. 成都：成都理工大學(xué)，2015：35-36.

[15]陳立華，李立豐，吳 福，等. 基于GIS與信息量法的北流市地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)[J]. 地球與環(huán)境，2020，48（4）：471-479.

[16]梁麗萍，劉延國(guó)，唐自豪，等. 基于加權(quán)信息量的地質(zhì)災(zāi)害易發(fā)性評(píng)價(jià)——以四川省瀘定縣為例[J]. 水土保持通報(bào)，2019，39（6）：176-182，321.