李華，顏坤，侍建國(guó)

(天津市水文水資源勘測(cè)管理中心，天津 300061)

0 引言

近年來(lái)，隨著計(jì)算機(jī)技術(shù)的日益發(fā)展，國(guó)內(nèi)的地下水模擬技術(shù)取得了較大的進(jìn)步，但缺乏大區(qū)域長(zhǎng)時(shí)間序列的高精度地下水?dāng)?shù)值模型評(píng)價(jià)地下水資源的相關(guān)研究。該文以天津市為例，全面介紹了利用Modflow建立大區(qū)域長(zhǎng)時(shí)間序列數(shù)值模型的過(guò)程，為相關(guān)技術(shù)研究提供技術(shù)支撐，為深層地下水資源評(píng)價(jià)提供科學(xué)依據(jù)。

在一次海洋基座式起重機(jī)回轉(zhuǎn)平臺(tái)物資質(zhì)量監(jiān)造中，勝利油田的設(shè)備監(jiān)造人員發(fā)現(xiàn)龍骨之間和封板角焊縫焊腳高度與圖紙要求存在1毫米偏差后，立即要求廠家對(duì)其焊縫打磨后重新焊接，否則不允許進(jìn)行下道工序。

天津市淺層地下水研究程度較高，而深層地下水研究程度相對(duì)較低，曾開(kāi)展過(guò)2次區(qū)域水資源綜合評(píng)價(jià)，第一次為1956—1984年，評(píng)價(jià)內(nèi)容包括淺層地下水資源量、可開(kāi)采量和深層地下水可開(kāi)采量；第二次為1980—2000年，評(píng)價(jià)內(nèi)容包括淺層地下水資源量和可開(kāi)采量。兩次評(píng)價(jià)的重點(diǎn)集中在淺層地下水資源，均未涉及深層地下水資源的構(gòu)成評(píng)價(jià)，影響了深層地下水資源的合理開(kāi)發(fā)利用。

最后，必須要定期的對(duì)學(xué)生進(jìn)行階段性考核，準(zhǔn)確有效地了解學(xué)生的學(xué)習(xí)進(jìn)展情況，然后依照學(xué)生的學(xué)習(xí)進(jìn)展對(duì)教學(xué)目標(biāo)不斷的進(jìn)行調(diào)整，以此來(lái)保證教學(xué)目標(biāo)的適用性，進(jìn)而保證高中化學(xué)分層教學(xué)模式能夠有效地發(fā)揮出應(yīng)有的作用。

傳統(tǒng)地下水資源的評(píng)價(jià)一般只評(píng)價(jià)地下水的資源量和可開(kāi)采量，難以進(jìn)行不同地下水開(kāi)發(fā)利用方案下的地下水不同要素的預(yù)測(cè)和分析。采用新技術(shù)和手段進(jìn)行水資源調(diào)查評(píng)價(jià)，全面分析深層地下水資源構(gòu)成，揭示其戰(zhàn)略資源屬性特征，評(píng)價(jià)可開(kāi)采資源量，提高地下水資源管理效率，提高地下水資源管理水平成為亟待解決的問(wèn)題。

歷史模型初始水位采用1980年3月底水位值，深層地下水水位監(jiān)測(cè)井較少，水位監(jiān)測(cè)成果難以滿足模擬需求，同時(shí)此時(shí)間段內(nèi)深層地下水開(kāi)采量比較穩(wěn)定，水位變幅相對(duì)較小，因此在1980年3月底深層地下水長(zhǎng)觀孔監(jiān)測(cè)水位成果的基礎(chǔ)上，加入地下水統(tǒng)測(cè)孔水位數(shù)值。利用克里金插值法[1]，得到平原區(qū)第Ⅱ～Ⅵ承壓含水組的初始水位值，針對(duì)東南部區(qū)域第Ⅳ和第Ⅴ承壓含水組水位觀測(cè)孔更少的情況，利用已有的鉆孔標(biāo)高進(jìn)行修正，結(jié)合水位外推，確定Ⅳ和第Ⅴ承壓含水組初始水位。第Ⅵ承壓含水組及以下地下水無(wú)開(kāi)采量，水位動(dòng)態(tài)變化趨勢(shì)很小，采用1997年12月水位值作為初始水位，采用1997年12月31日水位值作為現(xiàn)狀模型初始水位。

1 深層地下水流—地面沉降耦合模型的建立

1.1 時(shí)空離散

1.1.1 空間剖分

原始洞窟壁畫(huà)中就有對(duì)動(dòng)物油脂、礦物質(zhì)粉末、巖石、木炭等材料的使用，而文藝復(fù)興時(shí)期藝術(shù)家所使用的蛋清、石灰、亞麻仁油等，在其所處的時(shí)代也無(wú)疑是一種“新鮮的材料”。材料使用的變化一直在促使藝術(shù)概念的變化，因?yàn)槊糠N材料因其與生俱來(lái)的物性不同，運(yùn)用的創(chuàng)作方法不同，必然導(dǎo)致其背后所蘊(yùn)含的信息也就不同。

紅色文化在高校校園文化中占比不足，紅色教育活動(dòng)開(kāi)展較少，教育形式也局限于征文比賽、演講比賽、合唱比賽、觀影活動(dòng)等。受時(shí)間、經(jīng)費(fèi)等因素的限制，去烈士陵園掃墓、參觀紅色教育基地等活動(dòng)機(jī)會(huì)更是少之又少，導(dǎo)致學(xué)生接受紅色教育途徑少，校園紅色教育實(shí)效性差。除了線下活動(dòng)，據(jù)調(diào)查，各大高校少有建立專門(mén)的紅色教育網(wǎng)站或在官網(wǎng)上有紅色專題宣傳板塊，即便已建立的紅色教育網(wǎng)站也存在內(nèi)容陳舊、更新慢、吸引力不強(qiáng)等問(wèn)題，點(diǎn)擊量少，不能充分發(fā)揮網(wǎng)絡(luò)紅色教育作用，教育體系不夠完善[4]。

研究區(qū)域?yàn)樘旖蚴姓麄€(gè)平原區(qū)，面積約1.07×104km2?？臻g剖分采用有限差分法，使用GMS軟件中的Modflow模塊進(jìn)行矩形剖分，單元格500m×500m，共剖分7層346行248列，活動(dòng)單元格總數(shù)276111個(gè)，其中II層活動(dòng)單元格42628個(gè)，Ⅲ～Ⅵ層每層活動(dòng)單元格38171個(gè)(Ⅵ層以下概化為Ⅶ層，作為越流邊界層)。

1.1.2 時(shí)間離散

模擬期為1980—2013年，考慮到模擬期長(zhǎng)達(dá)34年，模擬初始時(shí)間段和截止時(shí)間段內(nèi)資料的詳細(xì)程度有較大差異，為了便于調(diào)整驗(yàn)證模型，將1980—2013年分為2個(gè)模擬期。1980年3月—1998年3月，由于歷史資料相對(duì)缺乏，以每3個(gè)月為一個(gè)應(yīng)力期，共72個(gè)應(yīng)力期，稱為歷史模型；1998年1月—2013年12月，資料較為豐富，以1個(gè)月為應(yīng)力期，共192個(gè)應(yīng)力期，稱為現(xiàn)狀模型。

1.2 數(shù)據(jù)處理

1.2.1 定解條件處理

(1)初始水位處理

秧床整地作畦要求上細(xì)下粗、上實(shí)下松。在播種前施肥并耕翻 1～2次，耕翻應(yīng)耕深15cm以上，做成畦寬1.2～1.4m、溝寬50cm、溝深20～30cm的畦，苗床四周開(kāi)好圍溝，以保證排水通暢。畦整平后用2.0kg壯秧劑加干細(xì)土20kg充分混拌可撒施25m2苗床，再反復(fù)耙勻于0～5cm深的土層內(nèi)，用木板壓平，澆透水即可播種。

(2)邊界條件處理

這種理論在比較法上的支持是德國(guó)的實(shí)踐。早在1902年，德國(guó)帝國(guó)法院即在“Duotal Lagend”案中確認(rèn)了對(duì)于根據(jù)專利方法(一種生產(chǎn)制造專利方法)直接獲得的產(chǎn)品，在權(quán)利人投入市場(chǎng)后，在其上的專利保護(hù)的效力就用盡了[注]Entscheidungen des Reichsgerichts in Zivilsachen, Erster Band, Der ganzen Reihe einundfünfzigster Band., Verlag von Beit und Comp. 1903, 139.。

1.2.5 地下水開(kāi)采量數(shù)據(jù)的處理

系統(tǒng)上邊界確定為潛水含水層自由水面，潛水通過(guò)系統(tǒng)上邊界與周邊環(huán)境產(chǎn)生垂向上的越流交換，各含水組之間也存在垂向上的越流交換。底界與Ⅵ含水組有水量交換，定義為水位邊界。

粘性土層概化的合理性直接關(guān)系地面沉降計(jì)算結(jié)果。將天津市平原區(qū)粘性土層厚度小于或者等于1.5m的定義為非滯后粘性土層，力學(xué)參數(shù)包括非彈性骨架儲(chǔ)水系數(shù)和彈性骨架儲(chǔ)水系數(shù)。粘性土層厚度大于1.5m的定義為滯后性粘性土層，力學(xué)參數(shù)包括垂向滲透系數(shù)、非彈性骨架儲(chǔ)水系數(shù)、彈性骨架儲(chǔ)水系數(shù)、粘土土層的有效厚度及比例因子。將非滯后性粘性土層和不同承壓含水層組的砂層均概化為非滯后壓縮層，在同一承壓含水層組中的滯后性粘性土層處理為滯后性壓縮層。

將初始沉降量設(shè)為0來(lái)滿足模型的擬合需要。考慮到粘性土的彈性以及塑性變形，引入前期固結(jié)壓力概念。當(dāng)前期最低水位高于當(dāng)前地下水水位時(shí)，粘性土體將產(chǎn)生塑性變形變化，當(dāng)前期最低水位低于當(dāng)前地下水水位時(shí)，粘性土體將產(chǎn)生彈性變形變化，通過(guò)地下水水位觀測(cè)孔監(jiān)測(cè)到的水位值得到各承壓含水組的最低水位。

針對(duì)偏遠(yuǎn)、交通不便的村民辦證和投訴不方便等情況，依托在鄉(xiāng)鎮(zhèn)、村委設(shè)立“便民服務(wù)點(diǎn)”，發(fā)放《便民服務(wù)聯(lián)系卡》，實(shí)行上門(mén)服務(wù)。

1.2.2 源匯項(xiàng)處理

利用GIS軟件的空間分析，利用程序?qū)懭氲较鄬?duì)應(yīng)的Modflow格式文件中，詳細(xì)情況見(jiàn)表1。

雜志里講的正是一個(gè)中年男人的情感故事，他邊看邊和自己做著比較，正看得入神，門(mén)忽然開(kāi)了，走進(jìn)幾個(gè)身穿制服的人，夾著公文包，面色冷峻，其中一個(gè)說(shuō)道：我們是專案組的，現(xiàn)在決定對(duì)你實(shí)行刑事拘留。兩個(gè)年輕力壯的人遂上前把他一夾，兩只手腕間多了一副錚亮的手銬，閃著耀眼的白光。他在白光中睜開(kāi)眼，天色大亮，原來(lái)是一場(chǎng)夢(mèng)。他渾身濕漉漉的，用手胡亂抹著汗水。妻子那邊的被子掀開(kāi)，衣服和挎包都不見(jiàn)了，枕頭下也沒(méi)有手機(jī)?；叵雺?mèng)中的情節(jié)，尹愛(ài)群的脊背又是陣陣寒氣。

表1 源匯項(xiàng)處理對(duì)照

1980—2013年平均排泄量為6.90×108m3/a，開(kāi)采量為6.28×108m3/a，占排泄量的90%～93%；側(cè)向流出量為0.31×108m3/a，占排泄量的4.44%；越流量為0.13×108m3/a，占排泄量的1.90%；壓縮釋水量為0.19×108m3/a，占排泄量的2.72%(表2)。

依據(jù)平原區(qū)各項(xiàng)水文地質(zhì)參數(shù)的不同分區(qū)，將降水入滲系數(shù)、各承壓含水組滲透系數(shù)、田間灌溉入滲補(bǔ)給系數(shù)、潛水含水層給水度和承壓含水組儲(chǔ)水率等不同水文地質(zhì)參數(shù)的分區(qū)初始值輸入到建立的模型中，通過(guò)擬合地下水不同承壓含水層組典型水位觀測(cè)孔和流場(chǎng)動(dòng)態(tài)變化過(guò)程線，識(shí)別不同承壓含水組的水文地質(zhì)參數(shù)，確定各水文地質(zhì)參數(shù)的分區(qū)數(shù)值。

1.2.4 粘性土層概化及壓縮參數(shù)處理

(3)沉降的定解條件處理

寶坻斷裂南部潛水含水層多含咸水，基本處于不開(kāi)采狀態(tài)，地下水水位年內(nèi)變化幅度較??；斷裂北部潛水含水層處于地下水系統(tǒng)的補(bǔ)給區(qū)，水位恢復(fù)較快且變幅較小，故不考慮潛水層的壓縮釋水。上層含水層的越流補(bǔ)給是各含水組的主要補(bǔ)給來(lái)源，如地下水開(kāi)采量較大，則水位年內(nèi)變幅也較大，且呈現(xiàn)多年持續(xù)下降趨勢(shì)，因此將Ⅱ～Ⅵ含水組作為壓縮釋水層。

根據(jù)不同含水層的砂層及厚度、粘性土層的層數(shù)及厚度，計(jì)算出每一個(gè)地下水水位觀測(cè)孔的滯后壓縮層和非滯后壓縮層對(duì)應(yīng)的參數(shù)，利用克里金插值法，得到每一個(gè)模型網(wǎng)格的參數(shù)數(shù)據(jù)，應(yīng)用VBA生成Modflow2005中SUB模塊所需要的.sub文件。

北部山區(qū)和平原區(qū)分界線、環(huán)渤海海岸帶和寶坻斷裂為自然邊界，其余均為人為邊界。流量邊界確定為山區(qū)和平原區(qū)分界線，渤海灣海岸帶以水頭值為0的水頭邊界處理，平原區(qū)北部的寶坻斷裂上下部有水量交換關(guān)系，并非嚴(yán)格意義上的不透水邊界，和其他人為邊界統(tǒng)一處理為三類邊界。

建模采用的核心模塊為Modflow-2005，但在處理大區(qū)域模型時(shí)，存在運(yùn)行時(shí)間長(zhǎng)、開(kāi)采量只能處理成點(diǎn)狀等缺點(diǎn)。通過(guò)改進(jìn)RCH子程序包，開(kāi)發(fā)出新的子程序包RAW(rechargeand well，RAW)，實(shí)現(xiàn)了垂向上對(duì)多層單元格補(bǔ)排量的表達(dá)，有效壓縮了源匯項(xiàng)文件，減少了大區(qū)域數(shù)值模型所需要讀取的數(shù)據(jù)量，提高了模型模擬的工作效率。

1.3 模型識(shí)別驗(yàn)證

模型模擬的地面沉降[10-19]與實(shí)測(cè)的地面沉降增長(zhǎng)趨勢(shì)基本一致，其中塘沽分層標(biāo)、西青分層標(biāo)和軍糧城分層標(biāo)監(jiān)測(cè)點(diǎn)地面沉降模擬值與實(shí)測(cè)值相差較小，大直沽天津堿廠分層標(biāo)模擬值與實(shí)測(cè)值相差較大，總體均表現(xiàn)為地面沉降值為增長(zhǎng)的趨勢(shì)(圖5)。

1.3.1 地下水流場(chǎng)

模型擬合校正過(guò)程要考慮2個(gè)重要過(guò)程，一是校正不同水文地質(zhì)參數(shù)，二是推算壓縮釋水量和反演地下水開(kāi)采量。通過(guò)不斷調(diào)整水文地質(zhì)參數(shù)，進(jìn)行地下水流場(chǎng)[7-9]、觀測(cè)孔過(guò)程線與實(shí)測(cè)數(shù)據(jù)的擬合對(duì)比。識(shí)別后的流場(chǎng)擬合整體效果較好，各含水組實(shí)測(cè)水位和模擬水位總體變化趨勢(shì)一致，單觀測(cè)井實(shí)測(cè)水位和模擬水位變化趨勢(shì)也基本一致(圖1—圖4)。

圖1 1997年12月Ⅳ含水組流場(chǎng)擬合

圖2 1997年12月Ⅴ含水組流場(chǎng)擬合

圖3 1980-1998年寶坻區(qū)BD18觀測(cè)井(Ⅱ含水組)水位過(guò)程線擬合

圖4 1998-2013年靜海區(qū)JH36觀測(cè)井(Ⅲ含水組)水位過(guò)程線擬合

1.3.2 沉降量的擬合

對(duì)建立的大區(qū)域、多層結(jié)構(gòu)、長(zhǎng)時(shí)間序列地下水流-地面沉降數(shù)值模型，應(yīng)用試估-校正法進(jìn)行識(shí)別和檢驗(yàn)[2-6]，總體原則如下：①模型識(shí)別的水文地質(zhì)參數(shù)均要符合實(shí)際的水文地質(zhì)條件；②模擬的地下水流場(chǎng)與實(shí)際的地下水流場(chǎng)要保持基本一致；③模擬的地下水均衡變化情況要基本符合研究區(qū)域?qū)嶋H的水均衡變化情況；④模擬的地下水水位動(dòng)態(tài)變化趨勢(shì)與實(shí)測(cè)的水位動(dòng)態(tài)變化趨勢(shì)基本相似。

圖5 2011—2012年地面沉降量擬合

2 深層地下水資源評(píng)價(jià)

2.1 地下水均衡項(xiàng)

補(bǔ)給項(xiàng)[20]包括通用水頭邊界流入量、給定水頭邊界流入量、山前側(cè)向補(bǔ)給量、越流補(bǔ)給量及壓縮釋水量，排泄項(xiàng)包括通用水頭邊界流出量、給定水頭邊界流出量、人工開(kāi)采量、越流排泄量及壓縮釋水量。

2.2 地下水均衡分析

1980—2013年平均補(bǔ)給量為6.46×108m3/a，其中Ⅱ和Ⅶ含水組越流補(bǔ)給量為2.72×108m3/a，占補(bǔ)給量的33.82%；粘性土壓縮釋水量為1.48×108m3/a，占補(bǔ)給量的18.36%；彈性釋水量為0.09×108m3/a，占補(bǔ)給量的1.09%；北部邊界及邊界的側(cè)向補(bǔ)給量為2.17×108m3/a，占補(bǔ)給量的26.97%。層間越流、側(cè)向流入、壓縮釋水已經(jīng)成為深層地下水的主要補(bǔ)給來(lái)源。

1.2.3 水文地質(zhì)參數(shù)處理

總之，運(yùn)用多媒體輔助語(yǔ)文教學(xué)，可以說(shuō)是有利有弊，關(guān)鍵在于教師的把握，在于教師的思想態(tài)度。用好多媒體可以促進(jìn)教學(xué)，提高教學(xué)效果，既有利于師，也有利于生；如果使用不當(dāng)多媒體，可能會(huì)導(dǎo)致教學(xué)效果差，學(xué)生反感，教師尷尬，先進(jìn)的設(shè)備變成害生的毒藥。

深層地下水均衡差為-0.44×108m3/a，深層地下水系統(tǒng)處于負(fù)均衡狀態(tài)。

表2 1980—2013年平均深層地下水資源組成(104m3/a)

2.3 可開(kāi)采資源量評(píng)價(jià)

(1)評(píng)價(jià)原則

南水北調(diào)通水之前，天津深層地下水開(kāi)采量遠(yuǎn)超過(guò)其獲得的補(bǔ)給量，從水均衡的角度將外界補(bǔ)給量作為深層地下水的可開(kāi)采資源量，結(jié)果相對(duì)偏小，也不符合實(shí)際情況。因此，深層地下水可開(kāi)采資源量的評(píng)價(jià)原則是：一方面要不產(chǎn)生嚴(yán)重環(huán)境地質(zhì)問(wèn)題，另一方面要結(jié)合天津深層地下水開(kāi)發(fā)利用現(xiàn)狀及歷史，人為給定環(huán)境約束條件。

(2)評(píng)價(jià)方法

主要考慮不同區(qū)域深層地下水的不同開(kāi)發(fā)利用程度[21]。開(kāi)發(fā)利用程度相對(duì)較高的區(qū)域，產(chǎn)生地質(zhì)災(zāi)害的可能性就越大，禁止利用深層地下水釋水量，可開(kāi)采資源量為側(cè)向徑流補(bǔ)給量和越流補(bǔ)給量?jī)身?xiàng)之和。相反，開(kāi)發(fā)利用程度相對(duì)較低的區(qū)域，產(chǎn)生地質(zhì)災(zāi)害的可能性就越小，可開(kāi)采資源量在考慮側(cè)向徑流補(bǔ)給量和越流補(bǔ)給量的基礎(chǔ)上，結(jié)合現(xiàn)狀開(kāi)采條件，要考慮可動(dòng)用的釋水量。

深層地下水可利用總量為3.08×108m3/a，其中Ⅱ含水組為2.20×108m3/a，Ⅲ含水組為0.52×108m3/a，Ⅳ含水組為0.26×108m3/a，Ⅴ含水組為0.10×108m3/a。

3 南水北調(diào)通水后地下水流場(chǎng)和地面沉降變化趨勢(shì)研究

3.1 地下水關(guān)鍵控制水位

利用數(shù)值模擬進(jìn)行地下水資源評(píng)價(jià)，通常分2個(gè)階段進(jìn)行，一是解逆問(wèn)題(也稱反演)，二是解正問(wèn)題(也稱正演)。解逆問(wèn)題是根據(jù)水位、開(kāi)采量等資料進(jìn)行數(shù)值計(jì)算，求出水文地質(zhì)參數(shù)，驗(yàn)證數(shù)值模型。解正問(wèn)題是利用解逆問(wèn)題的參數(shù)，計(jì)算不同開(kāi)采方案時(shí)的水位值，進(jìn)行方案對(duì)比和水位預(yù)報(bào)。

通過(guò)建立研究區(qū)地下水流-地面沉降耦合數(shù)值模型，將地下水開(kāi)采量作為約束條件代入模型，確定地下水關(guān)鍵性控制水位。以天津市為例，根據(jù)天津市地下水中長(zhǎng)期規(guī)劃，在壓采條件下求解出不同目標(biāo)年下的關(guān)鍵控制水位，計(jì)算出2020年和2030年平均紅線水位分別為-42.44m和-37.86m，用數(shù)值模擬法確定關(guān)鍵控制水位，是一種全面考慮水文地質(zhì)因素，計(jì)算精度較高的方法。

3.2 地下水流場(chǎng)變化趨勢(shì)

截至2017年12月27日，天津市累計(jì)受水23.7億m3。隨著最嚴(yán)格水資源管理制度的實(shí)施，地下水壓采進(jìn)程加快，2016年，深層地下水開(kāi)采量1.76億m3，比深層地下水可開(kāi)采資源量少1.32億m3。由于大幅度壓采深層地下水，深層地下水各承壓含水組水位將呈現(xiàn)整體上升趨勢(shì)。深層地下水系統(tǒng)將處于正均衡狀態(tài)，不同含水組儲(chǔ)存變化量為正值，彈性水在含水組中儲(chǔ)存了起來(lái)；壓縮釋水凈補(bǔ)給量為負(fù)值，多余的地下水儲(chǔ)存到土骨架中，儲(chǔ)存量處于逐漸恢復(fù)過(guò)程中。

ZHANG Ming-ming, ZHAO Pei, YU Yue-qing, ZHANG Cui-gai, GAO Wei, LI Fang

東麗區(qū)和津南區(qū)地下水水位有先下降后上升的變化趨勢(shì)，由于深層地下水開(kāi)采量仍稍大于各區(qū)可開(kāi)采量，在初始時(shí)間段內(nèi)，地下水水位有小幅下降，隨著周邊地區(qū)水位的持續(xù)上升，各區(qū)之間的內(nèi)部地下水交換量逐漸增大，水位也將逐漸上升。由于下層含水組水量主要由上層含水組越流補(bǔ)給，因此與上層相比，下層含水組得到的越流補(bǔ)給量相對(duì)越小，下層含水組水位的恢復(fù)速率比上層含水組也要小。

6.聘請(qǐng)張掖市螢火蟲(chóng)公益的成員為代理家長(zhǎng)。此公益組織主要是包含著個(gè)體老板和記者等多個(gè)主體，經(jīng)過(guò)合理的協(xié)調(diào)組合而成，他們一般是運(yùn)用周末或者是空閑時(shí)間，給村里的孩子們輔導(dǎo)課程，為其開(kāi)設(shè)講座等，讓孩子們擁有一個(gè)愉快的周末時(shí)光。對(duì)一些特別貧困的留守兒童進(jìn)行一對(duì)一資助，目前已經(jīng)有140多位留守兒童受到資助。

研究發(fā)現(xiàn)抑郁與服藥依從性密切相關(guān)[12]。腫瘤病人經(jīng)常出現(xiàn)情緒低落、悲觀、失望，特別是病情變化會(huì)加重病人恐懼、絕望心理，嚴(yán)重者甚至拒絕治療[13]。病人1：“我不看病了，看也白看，又治不好?！辈∪?：“我怎么會(huì)得這個(gè)病，治療這么長(zhǎng)時(shí)間，現(xiàn)在還是渾身沒(méi)勁。”病人8：“有這個(gè)病1年多了，這一年不是吃這藥就是吃那藥，都成藥罐子了，一說(shuō)吃藥就郁悶，就想發(fā)火”。

模型預(yù)測(cè)2018年各承壓含水組水位恢復(fù)較快，最大處水位上升4m左右，其余大部分地區(qū)水位上升2～3m。大部分區(qū)域地下水水位前期恢復(fù)較快，導(dǎo)致后期越流量減少，水位恢復(fù)呈現(xiàn)前快后慢，然后逐漸穩(wěn)定的趨勢(shì)。

第二，機(jī)械故障，導(dǎo)致這一故障出現(xiàn)的原因有很多，其中車輛的車門(mén)的相關(guān)部件由于擠壓等出現(xiàn)了變形，或者是在生產(chǎn)的過(guò)程中，尺寸的誤差比較大而使得硬件條件不符合相關(guān)要求是最重要的因素，同時(shí)在軌道車輛運(yùn)行的過(guò)程中，由于車輛車門(mén)的開(kāi)關(guān)頻率比較高，從而使得車門(mén)硬件條件發(fā)生變化也是導(dǎo)致這一類型故障出現(xiàn)的主要原因。在客流量比較大的時(shí)候，乘客在上下車的時(shí)候都會(huì)擠擠攘攘，很多乘客會(huì)被擠在車門(mén)的位置，長(zhǎng)此以往就會(huì)導(dǎo)致車門(mén)發(fā)生變形，這樣一來(lái)就會(huì)使得車門(mén)出現(xiàn)開(kāi)關(guān)遲鈍等常見(jiàn)故障，不利于軌道車輛的正常運(yùn)行。

3.3 地面沉降變化趨勢(shì)

當(dāng)大量開(kāi)發(fā)利用地下水時(shí)，孔隙水壓力降低，直接導(dǎo)致各含水組砂層壓密，當(dāng)孔隙水壓力恢復(fù)時(shí)，各含水組砂層基本能恢復(fù)原狀，砂礫類巖土呈現(xiàn)彈性變形的狀態(tài)。當(dāng)?shù)叵滤笠?guī)模超采，粘性土壓密后，其結(jié)構(gòu)將發(fā)生不可逆的變化，即使孔隙水壓力恢復(fù)，粘性土層也基本上繼續(xù)保持壓密狀態(tài)，以塑性變形為主。故而當(dāng)?shù)叵滤畨翰蓵r(shí)，各含水組砂層有回彈趨勢(shì)，粘性土層仍然處于壓密狀態(tài)。

在現(xiàn)有的壓采方案下模型預(yù)測(cè)平原區(qū)地面沉降沒(méi)有繼續(xù)加大的趨勢(shì)，但也沒(méi)有出現(xiàn)全面回彈的現(xiàn)象，截至2018年，除薊州和寧河區(qū)地面沉降速率在2mm/a之外，其余各區(qū)地面沉降處于涵養(yǎng)回彈態(tài)勢(shì)，地面沉降整體開(kāi)始回彈，回彈平均速率為3.42mm/a。

4 結(jié)論

(1)在建立平原區(qū)長(zhǎng)時(shí)間序列地下水流-地面沉降數(shù)值模型時(shí)，可根據(jù)數(shù)據(jù)精度不同，建立不同的連貫?zāi)Ｐ?，?shù)據(jù)精度低的時(shí)間段內(nèi)，模擬期可適當(dāng)延長(zhǎng)，數(shù)據(jù)精度高的時(shí)間段內(nèi)，模擬期適當(dāng)縮短，前后模型的模擬期可重復(fù)，但不可缺失。

(2)用試估-校正法進(jìn)行識(shí)別和檢驗(yàn)地下水?dāng)?shù)值模型時(shí)，水文地質(zhì)參數(shù)、地下水流場(chǎng)、地下水均衡變化和地下水水位動(dòng)態(tài)變化的擬合情況要符合研究區(qū)實(shí)際的水文地質(zhì)條件。模型的擬合校正過(guò)程要考慮2個(gè)重要過(guò)程，一是校正不同水文地質(zhì)參數(shù)，二是推算壓縮釋水量和反演地下水開(kāi)采量。

(3)用地下水?dāng)?shù)值模型進(jìn)行平原區(qū)深層地下水可開(kāi)采資源量評(píng)價(jià)時(shí)，基本原則是既要不產(chǎn)生嚴(yán)重環(huán)境地質(zhì)問(wèn)題，又要結(jié)合研究區(qū)深層地下水開(kāi)發(fā)利用現(xiàn)狀及歷史，人為給定環(huán)境約束條件，評(píng)價(jià)方法要考慮研究區(qū)不同區(qū)域深層地下水不同開(kāi)發(fā)利用程度，開(kāi)發(fā)利用程度低的區(qū)域要考慮可動(dòng)用的深層地下水釋水量。

(4)模型預(yù)測(cè)2018年平原區(qū)地下水水位恢復(fù)較快，后期呈現(xiàn)前快后慢，然后逐漸穩(wěn)定的趨勢(shì)。在現(xiàn)有的壓采方案下模型預(yù)測(cè)平原區(qū)地面沉降沒(méi)有繼續(xù)加大的趨勢(shì)，但也沒(méi)有出現(xiàn)全面回彈的現(xiàn)象。