近日，中科院測地所十三五“123”規劃重點(diǎn)培育方向、江利明研究員首席負責的“大地測量精密探測新理論與新技術(shù)”創(chuàng )新團隊，與中國水利水電科學(xué)研究院趙勇教授團隊合作，在利用衛星干涉雷達測量（InSAR）研究華北深層地下水儲量變化方面取得重要進(jìn)展，相關(guān)成果在水文水資源領(lǐng)域權威期刊Water Resources Research在線(xiàn)發(fā)表，并受邀將在2018年AGU年會(huì )上做口頭報告，論文題目為Combining InSAR and hydraulic head measurements to estimate aquifer parameters and storage variations of confined aquifer-system in Cangzhou, North China Plain（聯(lián)合InSAR及水頭觀(guān)測估算華北平原滄州地區含水層參數及水儲量變化）。
自20世紀70年代以來(lái)，華北平原深層地下水（承壓水）長(cháng)期超采形成世界最大的“地下水漏斗區”，導致大面積地面沉降、地下水資源枯竭和水環(huán)境惡化。因此，厘清深層地下水變化與地面沉降的耦合關(guān)系，查明華北平原深層地下水儲量變化，對于該地區地下水資源科學(xué)管理和調控，減緩與地下水超采相關(guān)的地質(zhì)環(huán)境問(wèn)題均具有重要意義。雖然目前重力衛星、地面水井監測及地下水建模用于估算地下水參數和儲量變化，但其分辨率、效率和精度不足以滿(mǎn)足當前區域尺度上水資源評價(jià)和調控的實(shí)際應用要求，而且上述方法均無(wú)法定量估算深層地下水儲量變化中的不可恢復部分。近年來(lái)，國際上利用InSAR觀(guān)測反演含水層參數方面取得了一些進(jìn)展，然而這些方法需要大量的水力學(xué)參數進(jìn)行約束（如前期預固結水頭等），難以在缺乏長(cháng)期、連續的歷史水頭和沉降資料的華北平原推廣應用。
針對上述問(wèn)題，江利明團隊在深入研究深層地下水變化與地表形變耦合關(guān)系的基礎上，提出了一種基于時(shí)序InSAR觀(guān)測的承壓含水層釋水系數及水儲量變化估算的新方法，并成功用于華北深層地下水超采最嚴重的滄州中部地區。該方法的優(yōu)勢在于能夠獲取高空間分辨率的承壓含水層骨架釋水系數，而且無(wú)需前期固結水頭等先驗知識，可有效分離深層地下水變化的可恢復部分（含水層彈性釋水）及不可恢復部分（弱透水層壓密釋水）。該研究不僅填補了華北平原地下水弱透水層壓密釋水研究方面的空白，也為定量探明華北平原深層地下水儲量動(dòng)態(tài)變化提供了新的手段。
研究發(fā)現，2003至2010年滄州地區呈現多個(gè)沉降漏斗（最大沉降速率達7.4 cm/yr），但有部分地區發(fā)生了明顯抬升（最大抬升速率約2.5 cm/yr）（圖1（a）），且季節性形變滯后水頭變化最大達175天（圖2）。此外，滄州中部地區深層地下水呈明顯的虧損趨勢（圖1b），年虧損率約為8.20±2.36×107 m3，其中大部分為弱透水層壓密釋水（圖1d），意味著(zhù)該時(shí)期地下水嚴重超采，甚至一些地區的含水層呈疏干狀態(tài)。值得注意的是，2005年以來(lái)滄州政府逐步控制深層地下水開(kāi)采并增加地表水供給，部分地區水位逐年回升，弱透水層壓密釋水虧損速率呈減緩趨勢（圖3），表明該地區地下水壓采措施對深層地下水超采起到了一定的減緩作用。
隨著(zhù)2014年12月南水北調中線(xiàn)工程正式向京津冀豫4省市實(shí)施生態(tài)補水，華北平原地下水預期會(huì )有一定程度的恢復，生態(tài)環(huán)境亦將改善，該團隊正在開(kāi)展聯(lián)合攻關(guān)，進(jìn)一步評估南水北調中線(xiàn)工程對該地區水資源和生態(tài)環(huán)境的影響。該項研究得到國家自然科學(xué)基金創(chuàng )新研究群體項目（課題）“華北平原含水層參數及地下水變化InSAR反演”和大地測量與地球動(dòng)力學(xué)國家重點(diǎn)實(shí)驗室重點(diǎn)項目“城市群地面沉降與地下水動(dòng)態(tài)耦合效應的大地測量集成監測研究”的聯(lián)合資助。
論文及支撐材料鏈接：https://agupubs.onlinelibrary.wiley.com/doi/abs/10.1029/2017WR022126
圖1 滄州中部2004至2010年年均形變速率（a）；深層地下水儲量虧損量（b）；可恢復的含水層彈性釋水（c）；不可恢復的弱透水層壓密釋水（d）。紫色實(shí)線(xiàn)為劃定的深層地下水虧損區域及恢復區域
圖2 典型井位季節性水頭與形變時(shí)間序列的相關(guān)性分析。地表形變滯后水頭變化為0至175天
圖3 滄州中部2004至2010年每年不可恢復水儲量虧損的等效水柱高空間分布。2005_2006表示2005年和2006年不可恢復水儲量虧損的平均值
