張政 金娟 郭興強(qiáng) 蔡保德

摘 要：磷石膏滲濾液對(duì)碳酸鹽巖及巖溶地下水產(chǎn)生重要影響，但其影響機(jī)制尚不清楚。本文開(kāi)展了磷石膏滲濾液對(duì)碳酸鹽巖的溶蝕試驗(yàn)，對(duì)比分析了不同pH值、水動(dòng)力條件及溶蝕時(shí)間對(duì)研究區(qū)不同碳酸鹽巖樣品的溶蝕效果，探討各試驗(yàn)條件對(duì)溶蝕作用的影響機(jī)制。結(jié)果表明：巖樣失重率呈現(xiàn)對(duì)數(shù)型增長(zhǎng)趨勢(shì);滲濾液pH值為1.86的試驗(yàn)組中巖樣失重率大小依次為婁山關(guān)組、法郎組、青巖組，pH值為2.71的試驗(yàn)組中巖樣失重率大小依次為：法郎組、青巖組、婁山關(guān)組;動(dòng)態(tài)試驗(yàn)失重率明顯高于靜態(tài)試驗(yàn)。滲濾液pH值和水動(dòng)力條件與其他變量的交互作用均表現(xiàn)強(qiáng)烈，滲濾液pH值、水動(dòng)力條件及二者間交互作用對(duì)巖樣溶蝕的影響在試驗(yàn)前期作用較大而在后期減弱，巖石層組類(lèi)別在試驗(yàn)后期作用凸顯。整體上看，試驗(yàn)變量對(duì)巖樣失重率的影響作用大小依次為滲濾液pH值、溶蝕時(shí)間、水動(dòng)力條件、巖石層組類(lèi)別。

關(guān)鍵詞：磷石膏滲濾液;碳酸鹽巖;溶蝕作用;室內(nèi)模擬試驗(yàn);巖溶地下水

中圖分類(lèi)號(hào)：X523?? 文獻(xiàn)標(biāo)識(shí)碼： A

貴州省地處喀斯特地區(qū)，以磷化工產(chǎn)業(yè)為重要支柱產(chǎn)業(yè)，因而關(guān)注和研究磷化工現(xiàn)存污染物對(duì)環(huán)境的影響和潛在隱患具有必要性和地區(qū)獨(dú)特性[1-2]。福泉市地處黔中腹地，號(hào)稱“亞洲磷都”，是我國(guó)大磷礦生產(chǎn)基地之一。

磷石膏滲濾液由工業(yè)副產(chǎn)物磷石膏堆存產(chǎn)生[3]，若防滲措施不完善，就會(huì)導(dǎo)致含磷、氟的滲濾液污染周邊地表水和地下水，破壞水生態(tài)系統(tǒng)的穩(wěn)定。福泉市擺紀(jì)堆場(chǎng)由于之前防滲不力已造成污染，發(fā)財(cái)洞為污染地下水主要出口。

國(guó)內(nèi)學(xué)者通過(guò)各種水文地質(zhì)手段對(duì)擺紀(jì)堆場(chǎng)的滲漏污染進(jìn)行了大量調(diào)查和研究確定滲漏污染的方式和途徑。前人研究測(cè)算磷石膏層滲透系數(shù)[4]，提出濕法磷石膏堆場(chǎng)的滲濾液產(chǎn)量計(jì)算公式[5]，建立擺紀(jì)磷石膏堆場(chǎng)巖溶滲漏污染預(yù)測(cè)模型[6]，預(yù)測(cè)堆場(chǎng)連續(xù)滲漏20年污染物濃度的空間分布[7]。磷石膏滲濾液從堆場(chǎng)至發(fā)財(cái)洞的具體滲漏途徑還未曾確定，滲濾液污染物質(zhì)在地下水中的運(yùn)移和水化學(xué)特征成為探究滲漏途徑的重要的研究方向和關(guān)注點(diǎn)，連通示蹤試驗(yàn)已確定發(fā)財(cái)洞方向?yàn)榱资喽褕?chǎng)滲漏主方向[8]，污染物的運(yùn)移路

徑呈現(xiàn)出狹窄的條帶狀，主要沿地下水徑流方向運(yùn)移[9]。堆場(chǎng)周邊已污染的地下水有從Ca2+-HCO-3型向Ca2+-SO2-4型演化的趨勢(shì)[10]，巖溶管道系統(tǒng)對(duì)污染物的吸附降解能力低下，滲濾液中的PO3-4和TFe在巖溶管道遷移過(guò)程中可能形成Fe（OH）3、羥磷灰石和藍(lán)鐵礦沉淀析出[11]。前人的研究普遍將重點(diǎn)放在了磷石膏滲濾液對(duì)地下水的污染上，還未涉及磷石膏滲濾液對(duì)巖溶含水介質(zhì)的溶蝕影響，但酸性磷石膏滲濾液進(jìn)入到巖溶管道與裂隙中后，會(huì)加速溶蝕碳酸鹽巖地層，擴(kuò)張巖溶含水介質(zhì)，甚至可能產(chǎn)生新的滲漏途徑，磷石膏滲濾液對(duì)碳酸鹽巖的溶蝕作用不容小覷。本研究選取三種處于滲漏途徑起點(diǎn)和終點(diǎn)附近地層的碳酸鹽巖，開(kāi)展室內(nèi)模擬試驗(yàn)，以pH值、水動(dòng)力條件、巖樣層組、溶蝕時(shí)間為試驗(yàn)變量，通過(guò)多因素方差分析比較各因素對(duì)溶蝕作用的影響，為磷石膏滲濾液的巖溶作用提供重要試驗(yàn)依據(jù)，為深入研究巖溶含水介質(zhì)的擴(kuò)張和新滲漏途徑的產(chǎn)生提供參考。

1 研究區(qū)概況

研究區(qū)磷石膏堆場(chǎng)坐落于貴州省黔南州的福泉市馬場(chǎng)坪鎮(zhèn)北東方向約2 km處。堆場(chǎng)位于馬場(chǎng)坪向斜西翼，出露地層為三疊系中統(tǒng)，堆場(chǎng)及周邊出露巖層青巖組二段（T2q2）、青巖組三段（T2q3）、法郎組一段（T2f1-1 和T2f1-2）、法郎組二段（T2f2），如圖1所示。研究區(qū)發(fā)育有擺郎發(fā)財(cái)洞巖溶管道、老落凼巖溶管道、鴨草壩巖溶管道、擺紀(jì)豆腐橋巖溶管道[7]。

法郎組一段下亞段（T2f1-1）和青巖組三段（T2q3）為相對(duì)隔水層，法郎組二段（T2f2）和青巖組二段（T2q2）地層為強(qiáng)含水層組，溶蝕裂隙、落水洞和巖溶管道發(fā)育較好，法郎組一段上亞段（T2f1-2）為中等含水層組，是發(fā)育出露泉水最多的地層，占研究區(qū)總?cè)c(diǎn)數(shù)的30%[12]。由渣場(chǎng)至發(fā)財(cái)洞A-A地質(zhì)剖面圖（圖2）可知，T2f2多表現(xiàn)為高臺(tái)地，T2f1-2表現(xiàn)為陡坡和斜坡更接近擺紀(jì)堆場(chǎng)，溶蝕路徑會(huì)由于地下水位和構(gòu)造線的控制向地層較深處延伸，堆場(chǎng)至發(fā)財(cái)洞的巖溶通道可能按法郎組一段上亞段（T2f1-2）的灰?guī)r地層走向發(fā)育。發(fā)財(cái)洞周邊區(qū)域地層為寒武系中上統(tǒng)婁山關(guān)組（∈2-3ls）白云巖，屬于弱含水層組，發(fā)財(cái)洞為滲漏出水口。

2 材料與方法

2.1 樣品采集方法

堆場(chǎng)滲漏點(diǎn)所在地層尚不明朗 [4，7，10]，根據(jù)實(shí)地調(diào)研，選取在堆場(chǎng)至發(fā)財(cái)洞之間分布的青巖組二段（T2q2）、法郎組一段上亞段（T2f1-2）和婁山關(guān)組（∈2-3ls）地層作為受試巖樣層組。于擺紀(jì)堆場(chǎng)東北側(cè)和擺紀(jì)坡腳分別采集法郎組和青巖組灰?guī)r巖樣，在發(fā)財(cái)洞西側(cè)巖壁采集婁山關(guān)組白云巖巖樣。采樣過(guò)程中選擇體積較大的巖石露頭，采集斷面新鮮、均勻，無(wú)侵蝕風(fēng)化痕跡。將采集巖石切割打磨成直徑38 mm，厚度10 mm的圓形薄片。各受試巖樣組的偏光顯微鏡鑒定情況見(jiàn)表1（下文中青巖組代表青巖組亮晶在巖;法郎組代表法郎組泥晶在巖;婁山關(guān)組代表婁山關(guān)組中晶白云巖），經(jīng)檢測(cè)青巖組和法郎組巖樣主要礦物組成為方解石和少量石英，婁山關(guān)組巖樣主要礦物組成為方解石和白云石。

磷石膏滲濾液的采樣點(diǎn)為磷石膏堆場(chǎng)下的滲濾液收集池，發(fā)財(cái)洞水樣為發(fā)財(cái)洞口出水，磷石膏滲濾液pH值在1.62～1.86之間（表2），滲濾液中PO3-4、SO2-4、F-的濃度較高，發(fā)財(cái)洞水樣可能受到地下水匯流的影響，其水溶性離子濃度與之相比明顯降低。Cl-基本不參與磷石膏滲濾液在遷移過(guò)程中的地球化學(xué)作用，因此根據(jù)Cl-在滲漏源頭（磷石膏滲濾液）和終點(diǎn)（發(fā)財(cái)洞水樣）中的濃度差距推斷地下水對(duì)磷石膏滲濾液的稀釋倍數(shù)在15～18倍左右。

2.2 室內(nèi)試驗(yàn)設(shè)計(jì)

以滲濾液初始pH值、水動(dòng)力條件、巖石層組類(lèi)別和溶蝕時(shí)間為試驗(yàn)變量設(shè)計(jì)試驗(yàn)方案（表3）。

試驗(yàn)用滲濾液通過(guò)稀釋磷石膏滲濾液原液獲得，并根據(jù)滲濾液原液與水稀釋比的差異控制試驗(yàn)用滲濾液初始pH值，分別模擬未與碳酸鹽巖接觸反應(yīng)、只考慮地下水匯流稀釋的磷石膏滲濾液經(jīng)近距離和遠(yuǎn)距離遷移后的狀態(tài)。第一組試驗(yàn)的滲濾液為磷石膏滲濾液原液和去離子水1∶1稀釋?zhuān)瑴y(cè)得pH值為1.86;第二組試驗(yàn)的稀釋比為1∶19，測(cè)得pH值為2.71。

兩組試驗(yàn)中再設(shè)置動(dòng)態(tài)試驗(yàn)組和靜態(tài)試驗(yàn)組。動(dòng)態(tài)試驗(yàn)組裝置采用自主設(shè)計(jì)的亞克力沖刷槽，3個(gè)層組的巖樣分置在不同沖刷槽內(nèi)，每槽設(shè)計(jì)容納28個(gè)巖樣，巖樣交替插入固定在沖刷反應(yīng)區(qū)的長(zhǎng)條方孔板上。按照每個(gè)巖樣以4 L滲濾液進(jìn)行溶蝕的反應(yīng)配比，每個(gè)沖刷槽配給滲濾液112 L。進(jìn)水流量為5.5～6 L/min，斷面平均流速約3～3.3 cm/s。靜態(tài)試驗(yàn)組容器為容量5 L的聚乙烯桶，向桶中注入4 L溶液，并利用小三腳架支撐巖樣，高度距離桶底約3 cm，以保證巖樣與滲濾液有充分的反應(yīng)接觸面積。

動(dòng)態(tài)試驗(yàn)巖樣的取樣時(shí)間點(diǎn)共有8個(gè)。靜態(tài)試驗(yàn)巖樣的取樣時(shí)間點(diǎn)有4個(gè)，與動(dòng)態(tài)試驗(yàn)對(duì)比分析。每個(gè)時(shí)間點(diǎn)取3個(gè)平行樣。

2.3 分析方法

巖樣溶蝕前用去離子水超聲清洗2次，每次10 min，洗去樣品表面的粉末和浮塵。105 ℃烘干10 h后，置于干燥皿中冷卻至室溫，利用電子天平稱重，記錄巖樣初始質(zhì)量m0。溶蝕后，先將其表面不溶物刷洗，再按照相同程序清洗、烘干、稱重，記錄巖樣溶蝕后質(zhì)量m1，計(jì)算巖樣失重率，并根據(jù)3個(gè)平行巖樣的失重率計(jì)算平均失重率。由于每片巖樣的質(zhì)量有輕微的差距，用失重率這一計(jì)算方法能更好的表現(xiàn)巖樣的被溶蝕程度，失重率計(jì)算如式（1）：

失重率=巖樣初始質(zhì)量m0-巖樣溶蝕后質(zhì)量m1巖樣初始質(zhì)量m0。（1）

3 結(jié)果與討論

3.1 巖樣失重率變化趨勢(shì)

從整體上看，失重率的增長(zhǎng)曲線隨時(shí)間的推移逐漸平緩，即巖樣的溶蝕速率逐漸變小。失重率的主要增長(zhǎng)區(qū)間在前15 d，溶蝕前期三種巖樣試片的失重率差距不大，15 d后失重率的增長(zhǎng)趨勢(shì)減緩，三個(gè)層組巖樣間失重率的差距逐漸擴(kuò)大（圖3）。

同組試驗(yàn)中，三種巖樣的失重率時(shí)間曲線的高低分布在動(dòng)態(tài)與靜態(tài)試驗(yàn)組中表現(xiàn)一致。第一組（pH=1.86）的巖樣失重率從大到小依次為：婁山關(guān)組、法郎組、青巖組，第二組（pH=2.71）的巖樣失重率從大到小依次為：法郎組、青巖組、婁山關(guān)組。法郎組巖樣的失重率整體大于青巖組，婁山關(guān)組巖樣在第一組中失重率最高，在第二組中失重率最低。第一組動(dòng)態(tài)組巖樣的失重率時(shí)間曲線兩兩間差異顯著（p<0.01），第一組靜態(tài)組中，法郎組與青巖組失重率時(shí)間曲線未達(dá)到顯著差異，而婁山關(guān)組與法郎組（p<0.05）、婁山關(guān)組與青巖組（p<0.01）間存在顯著差異。由于第二組巖樣失重率整體偏低，不同層組巖樣兩兩間皆未達(dá)到顯著性差異。

溶蝕試驗(yàn)60 d后，第一組中青巖組、法郎組和婁山關(guān)組的動(dòng)態(tài)試驗(yàn)和靜態(tài)試驗(yàn)的最大失重率分別為42.4%、46.8%、54.4%和36.6%、38.5%、41.2%，第二組中青巖組、法郎組和婁山關(guān)組的最大失重率依次為：動(dòng)態(tài)試驗(yàn)8.44%、8.01%和6.67%;靜態(tài)試驗(yàn)5.65%、6.33%和4.95%。其他試驗(yàn)變量相同、pH值不同時(shí)，第一組失重率是第二組的4.4～10.1倍;其他試驗(yàn)變量相同、水動(dòng)力條件不同時(shí)，動(dòng)態(tài)試驗(yàn)組與靜態(tài)試驗(yàn)組失重率之比為1.2～1.8。通過(guò)對(duì)比可知，室內(nèi)試驗(yàn)條件中滲濾液pH值對(duì)巖樣溶蝕作用的影響大于水動(dòng)力條件。

3.2 試驗(yàn)變量對(duì)失重率的影響分析

僅憑失重率比值的分析，對(duì)影響失重率的主要因素及貢獻(xiàn)程度的探討不夠深入全面，未包含其他試驗(yàn)變量和交互作用的影響，故利用多因素方差分析法對(duì)試驗(yàn)變量進(jìn)行全面分析。以試驗(yàn)變量滲濾液初始pH值（pH）、水動(dòng)力條件（Hyd）、巖石層組類(lèi)別（Lyr）和溶蝕時(shí)間（Dt）作為4個(gè)分類(lèi)變量，以溶蝕時(shí)間7、15、30和60 d的巖樣失重率數(shù)據(jù)為樣本總體，樣本總量共144個(gè)。

多次進(jìn)行多因素方差分析并逐步剔除不顯著因素后，決策系數(shù)R2為0.991，說(shuō)明已較為全面的囊括了對(duì)失重率有影響的因素（表4）。4個(gè)主效應(yīng)均呈現(xiàn)出顯著性（p<0.01），根據(jù)效應(yīng)量偏η2可發(fā)現(xiàn)主效應(yīng)對(duì)于巖樣失重率的影響作用大小依次為pH、Dt、Hyd、Lyr。交互作用中pH×Hyd、pH×Lyr、pH×Dt、Hyd×Lyr、pH×Hyd×Lyr對(duì)巖樣的失重率存在顯著影響（p<0.01），其中pH×Hyd、pH×Lyr、pH×Dt的效應(yīng)量較大，分別達(dá)到0.546、0.442、0.655，Hyd×Lyr、pH×Hyd×Lyr效應(yīng)量較小，分別為0.185、0.195。Dt在主效應(yīng)中的影響占第二位，但它與其他主效應(yīng)的交互效應(yīng)僅pH×Dt具有顯著性，說(shuō)明Dt因素對(duì)于溶蝕過(guò)程的影響以主效應(yīng)為主，pH在這一交互效應(yīng)中占據(jù)主導(dǎo)地位。

通過(guò)觀察不同時(shí)間階段效應(yīng)量的階段性變化，可以分析溶蝕過(guò)程中各因素對(duì)溶蝕作用的影響程度的變化狀況（圖4）。試驗(yàn)后期pH、Hyd、pH×Hyd效應(yīng)量占比減小，Lyr、pH×Lyr效應(yīng)量占比逐漸增加，Hyd×Lyr、pH×Hyd×Lyr則基本持平，說(shuō)明pH、Hyd及二者間交互作用對(duì)巖樣溶蝕的影響在試驗(yàn)前期發(fā)揮較大作用，隨試驗(yàn)的進(jìn)行影響減弱，而Lyr在試驗(yàn)后期作用逐漸凸顯。

pH和Hyd的交互作用表現(xiàn)為pH值較低的滲濾液會(huì)使得溶蝕效果更突出，經(jīng)水流沖刷剝離巖樣被溶蝕的表層后，會(huì)形成更大的比表面積，從而加大反應(yīng)接觸面積，進(jìn)一步加快溶蝕反應(yīng)速度，拉開(kāi)失重率差距，二者是協(xié)同作用。對(duì)于pH與Dt的交互作用（圖5（a）），可以看出pH對(duì)曲線的斜率即溶蝕反應(yīng)速度起決定性作用;隨著Dt的增加，第一組失重率邊際平均值的增幅減小，證明Dt對(duì)于溶蝕過(guò)程的貢獻(xiàn)逐漸減弱（圖5（b））。pH×Lyr（圖5（c））和Hyd ×Lyr（圖5（d））的共同點(diǎn)是代表青巖組與法郎組的線段幾乎平行，說(shuō)明二者間可能不存在這兩種交互作用，而婁山關(guān)組與二者迥異，說(shuō)明這兩種交互作用可能與白云巖的溶蝕特性有關(guān)，白云巖的溶蝕作用受pH和Hyd的影響較灰?guī)r更為敏感。現(xiàn)有研究表明，鹽酸中白云巖的溶蝕速率由溶蝕界面流速和pH值決定[13]，酸的種類(lèi)數(shù)和濃度[14]、同離子效應(yīng)[15-17]、晶格能[18]、酸不溶物[19]等因素均可能造成白云巖與灰?guī)r溶蝕速率的差異。

4 結(jié)論

（1）隨時(shí)間的推移，巖樣的溶蝕速率由大變小，失重率的變化表現(xiàn)為對(duì)數(shù)型增長(zhǎng)曲線。在pH值為1.86的試驗(yàn)組中婁山關(guān)組的巖樣失重率大于法郎組、青巖組，pH值為2.71的試驗(yàn)組中法郎組、青巖組的失重率大于婁山關(guān)組。法郎組失重率整體上大于青巖組。

（2）試驗(yàn)變量對(duì)巖樣失重率的影響作用大小依次為：滲濾液初始pH值、溶蝕時(shí)間、水動(dòng)力條件、巖石層組類(lèi)別。滲濾液pH值和水動(dòng)力條件與其他變量的交互作用均表現(xiàn)強(qiáng)烈，溶蝕時(shí)間則以主效應(yīng)影響為主。滲濾液pH值、水動(dòng)力條件及二者間交互作用對(duì)巖樣溶蝕的影響在試驗(yàn)前期發(fā)揮較大作用并隨著試驗(yàn)的進(jìn)行影響減弱，而巖石層組類(lèi)別在試驗(yàn)后期作用逐漸凸顯。婁山關(guān)組在不同試驗(yàn)組中表現(xiàn)出較大的失重率差異可能與其溶蝕特性有關(guān)，pH值和水流沖刷作用對(duì)婁山關(guān)組的影響可能比法郎組更大。

參考文獻(xiàn)：

[1]

湛亞禮. 從環(huán)境地質(zhì)學(xué)及環(huán)境地貌學(xué)的發(fā)展淺析貴州喀斯特環(huán)境問(wèn)題[J]. 科技風(fēng)， 2018（28）： 125.

[2]孫志立. 貴州省硫磷化工產(chǎn)業(yè)回顧與“十三五”展望[J]. 硫酸工業(yè)， 2017（1）： 8-16.

[3]李光明， 李霞， 賈磊， 等. 國(guó)內(nèi)外磷石膏處理和處置概況[J]. 無(wú)機(jī)鹽工業(yè)， 2012（10）： 11-13.

[4]褚學(xué)偉. 擺紀(jì)磷石膏堆場(chǎng)滲漏及子壩穩(wěn)定性分析[D]. 貴陽(yáng)： 貴州大學(xué)， 2006.

[5]蒙明富. 磷石膏堆場(chǎng)滲濾液產(chǎn)量預(yù)測(cè)研究[D]. 貴陽(yáng)： 貴州大學(xué)， 2016.

[6]綦婭. 擺紀(jì)磷石膏堆場(chǎng)巖溶滲漏污染預(yù)測(cè)[D]. 貴陽(yáng)： 貴州大學(xué)， 2006.

[7]張科正， 陳舟， 陳長(zhǎng)生， 等. 巖溶區(qū)某磷石膏堆場(chǎng)地下水滲流與污染物運(yùn)移數(shù)值模擬[J]. 中國(guó)煤炭地質(zhì)， 2018， 30（5）： 46-52.

[8]彭展翔， 褚學(xué)偉. 擺紀(jì)磷石膏堆場(chǎng)滲漏分析[J]. 地下水， 2012， 34（5）： 14-15， 87.

[9]鐘春里. 巖溶管道系統(tǒng)中污染物擴(kuò)散及地下水自凈能力研究： 貴州某磷石膏堆場(chǎng)為例[J]. 貴州地質(zhì)， 2011， 28（2）： 126-130.

[10]覃應(yīng)機(jī). 磷石膏堆場(chǎng)滲漏對(duì)巖溶地下水環(huán)境影響研究[D]. 貴陽(yáng)：貴州大學(xué)， 2017.

[11]查學(xué)芳， 覃應(yīng)機(jī)， 吳攀， 等. 磷石膏堆場(chǎng)滲漏影響下巖溶地下水地球化學(xué)過(guò)程[J]. 生態(tài)學(xué)雜志， 2018， 37（6）： 1708-1715.

[12]陳舟， 趙貴清， 王志光， 等. 巖溶區(qū)某磷石膏堆放場(chǎng)滲漏特征分析[J]. 水文地質(zhì)工程地質(zhì)， 2017， 44（2）： 144-150.

[13]GAUTELIER M， OELKERS E H， SCHOTT J. An experimental study of dolomite dissolution rates as a function of pH from 0.5 to 5 and temperature from 25 to 80 ℃[J]. Chemical Geology， 1999， 157（1-2）：13-26.

[14]莫春夢(mèng). 桂林市毛村流域碳酸鹽巖混合溶蝕實(shí)驗(yàn)研究[D]. 北京：中國(guó)地質(zhì)大學(xué)（北京）， 2019.

[15]WANG X， CHOU I M， HU W， et al. Kinetic inhibition of dolomite precipitation： insights from Raman spectroscopy of Mg2+-SO2-4 ion pairing in MgSO4-MgCl2-NaCl solutions at temperatures of 25 to 200 ℃ [J]. Chemical Geology， 2016， 435： 10-21.

[16]黃思靜， 楊俊杰， 張文正， 等. 石膏對(duì)白云巖溶解影響的試驗(yàn)?zāi)M研究[J]. 沉積學(xué)報(bào)， 1996， 14（1）： 103-109.

[17]閆志為. 硫酸根離子對(duì)方解石和白云石溶解度的影響[J]. 中國(guó)巖溶， 2008， 27（1）： 24-31.

[18]翁金桃. 方解石和白云石的差異溶蝕作用[J]. 中國(guó)巖溶， 1984（1）： 29-38.

[19]武鑫， 王藝霖， 黃敬軍， 等. 徐州地區(qū)碳酸鹽巖溶蝕特征及影響因素分析[J]. 地質(zhì)科技情報(bào)， 2019， 38（3）： 120-126.

（責(zé)任編輯：于慧梅）

Dissolution Experiment of Phosphogypsum Leachate on Carbonatite

ZHANG Zheng，? LI Jinjuan*，? GUO Xingqiang，? CAI Baode

（College of Resources and Environmental Engineering， Guizhou University， Guiyang 550025， China）

Abstract：

Phosphogypsum leachate has caused dissolution on carbonatite and led to karst groundwater pollution， however， its effect mechanism is not clear. In this study， the dissolution experiments of phosphogypsum leachate on Carbonatite were carried， dissolution effects of different carbonate rock samples under different pH values， hydrodynamic conditions and dissolution time were compared， and the influence mechanisms of these factors were discussed. The results indicate that the weight loss ratios of rock samples showed a logarithmic growth trend. The weight loss ratios from high to low was Loushanguan Formation， Falang Formation and Qingyan Formation at pH 1.86， and the sequence was Falang Formation， Qingyan Formation and Loushanguan Formation at pH 2.71. Otherwise dynamic weight-loss-ratio was higher than static obviously. Interactions between pH values of leachate， hydrodynamic conditions and other variables were strong. In the early stages of the dissolution， pH values of leachate， hydrodynamic conditions and their interactions played a major role on the impact of dissolution， but weakened in the later stages， while the effects of rock formation began to burgeon. As a whole， the effect of experimental variables on the weight loss rate of rock samples is as follows： pH value of leachate， dissolution time， hydrodynamic conditions， and rock group category.

Key words：

phosphogypsum leachate;carbonatite;dissolution;laboratory simulation experiment;karst groundwater

收稿日期：2020-04-25

基金項(xiàng)目：國(guó)家自然科學(xué)基金-貴州省喀斯特科學(xué)研究中心聯(lián)合資助項(xiàng)目（U1612442）

作者簡(jiǎn)介：張 政（1994-），男，碩士研究生，研究方向：環(huán)境工程，Email： zze531@outlook.com.

通訊作者：李金娟，Email： summy_lee@163.com.