本發(fā)明屬于煤礦采空區(qū)形變監(jiān)測(cè)，尤其涉及基于多源數(shù)據(jù)的煤礦采空區(qū)形變動(dòng)態(tài)耦合監(jiān)測(cè)方法及系統(tǒng)。

背景技術(shù)：

1、在礦業(yè)開采過(guò)程中，采空區(qū)的沉降監(jiān)測(cè)是確保開采過(guò)程安全的重要環(huán)節(jié)。若不對(duì)采空區(qū)的沉降進(jìn)行監(jiān)測(cè)，會(huì)引發(fā)采礦區(qū)災(zāi)害以及地表塌陷，為了保證開采過(guò)程的安全，故對(duì)采空區(qū)的沉降進(jìn)行監(jiān)測(cè)是十分必要的。

2、由于沉降監(jiān)測(cè)數(shù)據(jù)的準(zhǔn)確性和實(shí)時(shí)性對(duì)采空區(qū)的管理和安全保障至關(guān)重要；為了實(shí)時(shí)了解采空區(qū)的沉降情況，以便及時(shí)采取必要的措施，需定期對(duì)采空區(qū)的沉降進(jìn)行監(jiān)測(cè)，但目前對(duì)采空區(qū)的沉降進(jìn)行監(jiān)測(cè)時(shí)，監(jiān)測(cè)數(shù)據(jù)的采集頻率只能依靠預(yù)先設(shè)置的時(shí)間間隔進(jìn)行采集，但單一的數(shù)據(jù)采集頻率會(huì)影響到監(jiān)測(cè)的準(zhǔn)確性和實(shí)時(shí)性，當(dāng)數(shù)據(jù)采集頻率過(guò)高時(shí)，會(huì)導(dǎo)致數(shù)據(jù)冗余，增加數(shù)據(jù)處理成本；當(dāng)采集頻率過(guò)低時(shí)，可能會(huì)漏掉關(guān)鍵的沉降速度變化信息，從而影響監(jiān)測(cè)的準(zhǔn)確性，即如何根據(jù)礦山開采的實(shí)際情況動(dòng)態(tài)設(shè)置采集頻率，提高監(jiān)測(cè)數(shù)據(jù)的準(zhǔn)確性。

3、為解決上述問(wèn)題，發(fā)明專利一種采空區(qū)沉降監(jiān)測(cè)方法、裝置、電子設(shè)備及存儲(chǔ)介質(zhì)(公開號(hào)cn117870618a，公告日20240412)公開獲取沉降監(jiān)測(cè)設(shè)備對(duì)采空區(qū)進(jìn)行監(jiān)測(cè)的第一沉降速度數(shù)據(jù)，沉降監(jiān)測(cè)設(shè)備包括分層沉降監(jiān)測(cè)設(shè)備、次聲監(jiān)測(cè)設(shè)備、壓力傳感器設(shè)備、雷達(dá)測(cè)量設(shè)備以及地表沉降監(jiān)測(cè)設(shè)備中的任意一種沉降監(jiān)測(cè)設(shè)備；獲取采空區(qū)的第一采集頻率；對(duì)第一沉降速度數(shù)據(jù)進(jìn)行計(jì)算，得到沉降速度變化率；根據(jù)沉降速度變化率對(duì)第一采集頻率進(jìn)行調(diào)整，以便于沉降監(jiān)測(cè)設(shè)備對(duì)采空區(qū)的沉降進(jìn)行數(shù)據(jù)采集。

4、但現(xiàn)有專利在數(shù)據(jù)融合上，不能有效消除的大氣延遲誤差，對(duì)采空區(qū)的沉降精度仍有待提高。

5、在采空區(qū)形變監(jiān)測(cè)上，發(fā)明專利一種煤礦采空區(qū)形變敏感性評(píng)價(jià)方法(公開號(hào)cn117111057a，公開日20231124)，公開s1、建立煤礦采空區(qū)形變敏感性評(píng)價(jià)指標(biāo)庫(kù)；s2、利用sbas-insar技術(shù)獲取煤礦采空區(qū)地表形變速率；s3、在研究區(qū)內(nèi)均勻選取樣本點(diǎn)，以煤礦采空區(qū)形變敏感性評(píng)價(jià)指標(biāo)最為輸入變量，以相應(yīng)地表形變速率作為標(biāo)簽，構(gòu)建模型訓(xùn)練樣本數(shù)據(jù)集，并對(duì)其按照7:3的比例進(jìn)行隨機(jī)分組，分別作為訓(xùn)練集和測(cè)試集。該煤礦采空區(qū)形變敏感性評(píng)價(jià)方法，通過(guò)聯(lián)合地表形變探測(cè)技術(shù)、遙感技術(shù)和機(jī)器學(xué)習(xí)方法，可在一定程度上彌補(bǔ)常規(guī)基于調(diào)查統(tǒng)計(jì)資料的煤礦采空區(qū)形變敏感性評(píng)價(jià)方法的局限性。

6、但現(xiàn)有專利在采礦區(qū)的整體和高梯度的變形監(jiān)測(cè)和預(yù)測(cè)上精度有待進(jìn)一步提高。

7、針對(duì)形變耦合監(jiān)測(cè)，發(fā)明專利一種煤巖多場(chǎng)耦合監(jiān)測(cè)試驗(yàn)裝置(公開號(hào)cn107727508a，公開日20180223)公開通過(guò)對(duì)試樣持續(xù)穩(wěn)定的施加圍壓、軸向壓力、孔隙流壓以及對(duì)試樣進(jìn)行加熱，能夠模擬開采過(guò)程中長(zhǎng)期處在高地應(yīng)力、高滲透壓以及高地溫耦合作用下的煤巖體的受力情況，并且該監(jiān)測(cè)試驗(yàn)裝置能夠采集試樣整個(gè)受載過(guò)程中產(chǎn)生的各向應(yīng)變信號(hào)、聲發(fā)射信號(hào)以及電信號(hào)，根據(jù)各個(gè)信號(hào)可以得到試樣長(zhǎng)期演化破壞過(guò)程中的軸向和徑向應(yīng)變、試樣破壞過(guò)程中的聲發(fā)射信號(hào)和電荷信號(hào)的時(shí)空演化規(guī)律，更有利于建立實(shí)際開采過(guò)程中受到多場(chǎng)耦合作用的煤巖體的力學(xué)行為演化與失穩(wěn)災(zāi)變的作用機(jī)制，揭示煤巖體動(dòng)力災(zāi)害本源機(jī)理，為煤巖體動(dòng)力災(zāi)害防治提供更加可靠的試驗(yàn)基礎(chǔ)，根據(jù)試驗(yàn)基礎(chǔ)建立相應(yīng)的煤巖體失穩(wěn)災(zāi)變判定及防治技術(shù)。但現(xiàn)有技術(shù)對(duì)于礦區(qū)整體的大梯度變形監(jiān)測(cè)精度上仍需要進(jìn)一步提高。

8、通過(guò)上述分析，現(xiàn)有技術(shù)存在的問(wèn)題及缺陷為：現(xiàn)有技術(shù)對(duì)煤礦采空區(qū)井上井下形變動(dòng)態(tài)耦合監(jiān)測(cè)精度低，對(duì)地表變形監(jiān)測(cè)效果差。

技術(shù)實(shí)現(xiàn)思路

1、為克服相關(guān)技術(shù)中存在的問(wèn)題，本發(fā)明公開實(shí)施例提供了基于多源數(shù)據(jù)的煤礦采空區(qū)形變動(dòng)態(tài)耦合監(jiān)測(cè)方法及系統(tǒng)，尤其涉及一種基于深度學(xué)習(xí)的三維采礦區(qū)建模軟件的物性參數(shù)厚度確定方法及系統(tǒng)。

2、所述技術(shù)方案如下：基于多源數(shù)據(jù)的煤礦采空區(qū)形變動(dòng)態(tài)耦合監(jiān)測(cè)方法，包括：

3、s1，基于北斗gnss與d-insar融合技術(shù)的礦區(qū)沉降監(jiān)測(cè)，在井上采礦區(qū)建立gnss連續(xù)觀測(cè)站網(wǎng)點(diǎn)，構(gòu)建區(qū)域水汽模型和大氣層延遲誤差改正模型，完成時(shí)間域與空間域的融合模型和算法融合處理；采用雙插雙估計(jì)didp算法實(shí)現(xiàn)對(duì)gnss與d-insar監(jiān)測(cè)數(shù)據(jù)的加密，削弱大氣延遲誤差影響；

4、s2，針對(duì)獲取的削弱大氣延遲誤差影響后的礦區(qū)沉降監(jiān)測(cè)信息，對(duì)于采礦區(qū)沉降區(qū)域發(fā)生的大梯度形變產(chǎn)生的失相干，利用基于對(duì)數(shù)邏輯模型的sbas-insar與gnss耦合技術(shù)的采礦區(qū)變形監(jiān)測(cè)方法進(jìn)行計(jì)算采礦活動(dòng)在任何時(shí)刻產(chǎn)生的變形信息；

5、s3，分別部署電法系統(tǒng)和微震監(jiān)測(cè)系統(tǒng)進(jìn)行井下電阻率與微震聯(lián)合應(yīng)力變化監(jiān)測(cè)，實(shí)現(xiàn)井下采礦區(qū)大梯度異常體的圈定；結(jié)合井上采礦區(qū)沉降形變監(jiān)測(cè)信息，聯(lián)合井下電阻率與微震聯(lián)合應(yīng)力變化監(jiān)測(cè)，進(jìn)行協(xié)同監(jiān)測(cè)，完成對(duì)采空區(qū)形變過(guò)程協(xié)同動(dòng)態(tài)耦合監(jiān)測(cè)。

6、在步驟s1中，構(gòu)建區(qū)域水汽模型，包括：

7、以區(qū)域水汽飽和網(wǎng)格體每層的平均溫度、氣壓、水汽壓作為因變量，確定每層的水汽飽和分層因子，表達(dá)式為：

8、

9、式中，my(a)為區(qū)域水汽飽和網(wǎng)格體分層因子，a1為平均溫度的權(quán)重因子，a2為平均氣壓的權(quán)重因子，a3為平均水汽壓的權(quán)重因子，ta為溫度平均值，ωa為氣壓平均值，sa為水汽壓平均值，tm為第m層大氣溫度值，ωm為第m層大氣氣壓值，sm為第m層大氣水汽壓值。

10、在步驟s1中，構(gòu)建大氣層延遲誤差改正模型利用區(qū)域水汽飽和網(wǎng)格體每層傳輸通道間的相關(guān)性，通過(guò)矩形預(yù)測(cè)分別對(duì)大氣層圖像三個(gè)因變量每層傳輸通道進(jìn)行延遲像素值預(yù)測(cè)，組成預(yù)測(cè)延遲誤差組，在三維預(yù)測(cè)誤差直方圖中通過(guò)修改預(yù)測(cè)延遲誤差組的映射關(guān)系，舍棄嵌入稀疏變異大的映射，重建新映射，實(shí)現(xiàn)區(qū)域水汽飽和網(wǎng)格體的嵌入和區(qū)域水汽飽和網(wǎng)格體的提取。

11、進(jìn)一步，區(qū)域水汽飽和網(wǎng)格體的嵌入，包括：

12、步驟一，對(duì)原始大氣層圖像按照因變量每層傳輸通道進(jìn)行分離，對(duì)原大氣層圖像按照溫度、氣壓、水汽壓因變量每層傳輸通道分離三張灰度圖像，分別記為溫度層傳輸通道、氣壓層傳輸通道和水汽壓層傳輸通道；

13、步驟二，大氣層圖像預(yù)處理；

14、步驟三，計(jì)算每張灰度圖的預(yù)測(cè)延遲誤差組，對(duì)每張灰度圖采用矩形預(yù)測(cè)方案獲取模糊節(jié)點(diǎn)的預(yù)測(cè)值，跟原始延遲像素值相減，獲得模糊節(jié)點(diǎn)的預(yù)測(cè)誤差序列(et1，et2…etm)，(eω1，eω2…eωm)，(es1，es2…esm)，將三張灰度圖像模糊節(jié)點(diǎn)的預(yù)測(cè)延遲誤差組成預(yù)測(cè)延遲誤差組em＝(etm，eωm，esm)；其中，e為延遲誤差預(yù)測(cè)值；

15、步驟四，選擇用于嵌入?yún)^(qū)域水汽飽和網(wǎng)格體的預(yù)測(cè)延遲誤差組，計(jì)算每一個(gè)預(yù)測(cè)延遲誤差組的局部復(fù)雜度kn，找到預(yù)測(cè)延遲誤差組的臨界值c，被替換的lsbs、壓縮后的位圖和區(qū)域水汽飽和網(wǎng)格體消息，修改滿足kn≤c的預(yù)測(cè)延遲誤差組的映射來(lái)嵌入?yún)^(qū)域水汽飽和網(wǎng)格體；

16、步驟五，模糊節(jié)點(diǎn)集區(qū)域水汽飽和網(wǎng)格體嵌入，在模糊節(jié)點(diǎn)的預(yù)測(cè)誤差直方圖中嵌入?yún)^(qū)域水汽飽和網(wǎng)格體信息，稱為模糊節(jié)點(diǎn)集區(qū)域水汽飽和網(wǎng)格體嵌入；

17、步驟六，大氣層清晰圖像節(jié)點(diǎn)集區(qū)域水汽飽和網(wǎng)格體嵌入，在清晰點(diǎn)的預(yù)測(cè)誤差直方圖中嵌入?yún)^(qū)域水汽飽和網(wǎng)格體信息，實(shí)現(xiàn)大氣層清晰圖像節(jié)點(diǎn)集區(qū)域水汽飽和網(wǎng)格體嵌入，模糊節(jié)點(diǎn)集嵌入?yún)^(qū)域水汽飽和網(wǎng)格體信息之后，模糊延遲像素點(diǎn)yi，j被修改為y′i，j；用嵌入?yún)^(qū)域水汽飽和網(wǎng)格體之后的模糊節(jié)點(diǎn)的延遲像素預(yù)測(cè)清晰點(diǎn)的延遲像素值，得到三張灰度圖像清晰點(diǎn)的預(yù)測(cè)誤差直方圖，嵌入?yún)^(qū)域水汽飽和網(wǎng)格體信息；

18、步驟七，將嵌入?yún)^(qū)域水汽飽和網(wǎng)格體之后的三張灰度圖合并，得到嵌入?yún)^(qū)域水汽飽和網(wǎng)格體后大氣層圖像，得到基于大氣層延遲誤差改正模型的大氣層圖像。

19、在步驟二中，大氣層圖像預(yù)處理，包括：

20、第一步，對(duì)大氣層圖像三個(gè)因變量每層傳輸通道的延遲像素進(jìn)行區(qū)域水汽飽和網(wǎng)格體預(yù)處理，調(diào)整延遲像素值等于邊界值0或255的延遲像素點(diǎn)到合理范圍內(nèi)；

21、第二步，構(gòu)造相對(duì)應(yīng)的三個(gè)位圖分別記錄三個(gè)因變量每層傳輸通道中被修改的延遲像素位置，被修改過(guò)的位圖中標(biāo)記為0，沒(méi)有被修改過(guò)的標(biāo)為1，無(wú)損壓縮三個(gè)位圖；

22、第三步，清空大氣層圖像溫度層傳輸通道的第一行延遲像素，使用lsb算法嵌入三個(gè)參數(shù)：預(yù)測(cè)延遲誤差組的臨界值c，壓縮后的位圖尺寸和區(qū)域水汽飽和網(wǎng)格體消息大小制造空間，被替換的lsbs和壓縮后的位圖作為負(fù)載的一部分；

23、在步驟五中，模糊節(jié)點(diǎn)集區(qū)域水汽飽和網(wǎng)格體嵌入，包括：

24、第1步，在嵌入稀疏變異小于等于2的映射中嵌入?yún)^(qū)域水汽飽和網(wǎng)格體，舍棄嵌入稀疏變異大于2的映射(1,1,1)→(0,0,0)，重建(0,0,0)的新映射；

25、第2步，對(duì)于符合要求的映射進(jìn)行分類，分為5類，第2類為預(yù)測(cè)延遲誤差組(0,0,0)的新建映射：

26、第1類：

27、第2類：

28、第3類：

29、第4類：

30、第5類：(n1，n2，n3){(n1-1，n2-1，n3-1)

31、式中，n1為橫向映射的延遲誤差變異值，n2為縱向映射的延遲誤差變異值，n3為空間方向映射的延遲誤差變異值；

32、第3步，將預(yù)測(cè)延遲誤差組中的每個(gè)成員和它預(yù)測(cè)值相加，得到處理后的延遲像素值。

33、進(jìn)一步，區(qū)域水汽飽和網(wǎng)格體的提取，包括：

34、步驟1，對(duì)嵌入?yún)^(qū)域水汽飽和網(wǎng)格體后的大氣層圖像，按照溫度、氣壓、水汽壓因變量每層傳輸通道分離三張灰度；

35、步驟2，讀取大氣層圖像溫度層傳輸通道第一行延遲像素的lsbs，取出三個(gè)參數(shù)的值；

36、步驟3，提取大氣層清晰圖像節(jié)點(diǎn)集區(qū)域水汽飽和網(wǎng)格體；

37、步驟4，提取模糊節(jié)點(diǎn)集區(qū)域水汽飽和網(wǎng)格體，根據(jù)恢復(fù)的大氣層清晰圖像節(jié)點(diǎn)集的延遲像素值，使用相同方法獲取三個(gè)每層傳輸通道模糊節(jié)點(diǎn)的預(yù)測(cè)誤差序列，提取出嵌入負(fù)載，恢復(fù)原延遲像素值；

38、步驟5，恢復(fù)原圖。

39、在步驟3中，提取大氣層清晰圖像節(jié)點(diǎn)集區(qū)域水汽飽和網(wǎng)格體，包括：

40、第i步，使用相同的掃描順序和預(yù)測(cè)方法獲取三個(gè)每層傳輸通道清晰點(diǎn)的預(yù)測(cè)誤差序列，組成預(yù)測(cè)延遲誤差組；

41、第ii步，對(duì)每一個(gè)預(yù)測(cè)延遲誤差組，計(jì)算局部復(fù)雜度kn；

42、第iii步，處理滿足kn≤c的預(yù)測(cè)延遲誤差組，通過(guò)嵌入時(shí)的映射方案的逆映射恢復(fù)出原來(lái)的延遲像素值，根據(jù)映射時(shí)不同的類型恢復(fù)出負(fù)載；

43、第iv步，提取出嵌入負(fù)載后，獲取位圖和被替換的lsbs，根據(jù)位圖，處理標(biāo)記為0的延遲像素，恢復(fù)原值。

44、在步驟s2中，基于對(duì)數(shù)邏輯模型的sbas-insar與gnss耦合技術(shù)的采礦區(qū)變形監(jiān)測(cè)方法進(jìn)行計(jì)算采礦活動(dòng)在任何時(shí)刻產(chǎn)生的變形信息，包括：

45、采用采礦區(qū)大梯度形變重力異常鑒別方法，在一致的時(shí)間間隔內(nèi)集成gnss和sbas-insar數(shù)據(jù)；根據(jù)采礦區(qū)的塌陷特征，得到礦區(qū)整體的大梯度變形模式；將對(duì)數(shù)邏輯log-logistic模型與單個(gè)insar對(duì)相結(jié)合，計(jì)算采礦活動(dòng)在任何時(shí)刻產(chǎn)生的變形信息；

46、礦區(qū)大梯度形變重力異常鑒別方法采用采礦區(qū)大梯度形變距離賦閾值環(huán)形均勻布設(shè)運(yùn)動(dòng)軌跡路線插值分割方法，均勻選取環(huán)形均勻布設(shè)運(yùn)動(dòng)軌跡路線上的計(jì)算點(diǎn)；計(jì)算點(diǎn)的值基于其周邊的點(diǎn)值采礦區(qū)大梯度形變距離賦閾值計(jì)算得出。

47、進(jìn)一步，采礦區(qū)大梯度形變距離賦閾值環(huán)形均勻布設(shè)運(yùn)動(dòng)軌跡路線插值分割方法，包括：

48、設(shè)定重力異常場(chǎng)h(h，o)由區(qū)域場(chǎng)s(h，o)和局部場(chǎng)f(h，o)構(gòu)成；b(h，o)代表與中心點(diǎn)相距為d的某點(diǎn)的重力場(chǎng)平均值，表達(dá)式為：

49、

50、式中，m為重力場(chǎng)總節(jié)點(diǎn)數(shù)量，i為第i個(gè)重力場(chǎng)總節(jié)點(diǎn)hi(d)為以(h，o)點(diǎn)為圓心，直徑為d的環(huán)形均勻布設(shè)運(yùn)動(dòng)軌跡路線上點(diǎn)；

51、引入閾值系數(shù)f，利用重力異常h(h，o)與重力場(chǎng)平均值b(h，o)的差值來(lái)定義切割直徑d下的局部異常場(chǎng)，表達(dá)式為：

52、f(h，o)＝f[h(h，o)-b(h，o)]

53、式中，f(h，o)為切割直徑d下的局部異常場(chǎng)；

54、通過(guò)對(duì)梯度層位置場(chǎng)網(wǎng)格數(shù)據(jù)采取采礦區(qū)大梯度形變距離賦閾值法插值計(jì)算獲得，選取i點(diǎn)處距離小于單個(gè)網(wǎng)格點(diǎn)距的j個(gè)點(diǎn)參與計(jì)算，表達(dá)式為：

55、

56、

57、式中，h(d)為直徑為d的環(huán)形均勻布設(shè)運(yùn)動(dòng)軌跡路線，j為第j個(gè)網(wǎng)格點(diǎn)，vj為第j個(gè)網(wǎng)格點(diǎn)距離，dj為第j個(gè)網(wǎng)格點(diǎn)的直徑，m為重力場(chǎng)總節(jié)點(diǎn)數(shù)量；

58、此時(shí)，區(qū)域場(chǎng)s(h，o)和局部場(chǎng)異常場(chǎng)f(h，o)分別為切割直徑d所對(duì)應(yīng)的區(qū)域場(chǎng)和局部場(chǎng)，閾值系數(shù)f由圓心處、環(huán)形均勻布設(shè)運(yùn)動(dòng)軌跡路線上的梯度層位置場(chǎng)值共同確定，由下式計(jì)算所得：

59、

60、

61、式中，ui(d)為直徑為d的和hi(d)的差值比值，hi(d)為以(h，o)點(diǎn)為圓心，直徑為d的環(huán)形均勻布設(shè)運(yùn)動(dòng)軌跡路線上點(diǎn)的重力場(chǎng)值，為hi(d)在切割直徑為d的環(huán)形均勻布設(shè)運(yùn)動(dòng)軌跡路線上的對(duì)稱點(diǎn)；

62、ui(d)與和hi(d)的差值相關(guān)，具有對(duì)稱性，因此閾值系數(shù)f的計(jì)算簡(jiǎn)化為：

63、

64、則區(qū)域場(chǎng)s(h，o)為重力異常場(chǎng)h(h，o)與重力場(chǎng)平均值b(h，o)的賦閾值平均，表達(dá)式為：

65、s(h，o)＝(1-f)h(h，o)+fb(h，o)。

66、本發(fā)明的另一目的在于提供一種基于多源數(shù)據(jù)的煤礦采空區(qū)形變動(dòng)態(tài)耦合監(jiān)測(cè)系統(tǒng)，該系統(tǒng)實(shí)施基于多源數(shù)據(jù)的煤礦采空區(qū)形變動(dòng)態(tài)耦合監(jiān)測(cè)系統(tǒng)，該系統(tǒng)包括：

67、基于北斗gnss與d-insar融合技術(shù)的礦區(qū)沉降監(jiān)測(cè)模塊，用于在井上采礦區(qū)建立gnss連續(xù)觀測(cè)站網(wǎng)點(diǎn)，構(gòu)建區(qū)域水汽模型和大氣層延遲誤差改正模型，完成時(shí)間域與空間域的融合模型和算法融合處理；采用雙插雙估計(jì)didp算法實(shí)現(xiàn)對(duì)gnss與d-insar監(jiān)測(cè)數(shù)據(jù)的加密，削弱大氣延遲誤差影響；

68、采礦活動(dòng)在任何時(shí)刻產(chǎn)生的變形信息計(jì)算模塊，用于針對(duì)獲取的削弱大氣延遲誤差影響后的礦區(qū)沉降監(jiān)測(cè)信息，對(duì)于采礦區(qū)沉降區(qū)域發(fā)生的大梯度形變產(chǎn)生的失相干，利用基于對(duì)數(shù)邏輯模型的sbas-insar與gnss耦合技術(shù)的采礦區(qū)變形監(jiān)測(cè)方法進(jìn)行計(jì)算采礦活動(dòng)在任何時(shí)刻產(chǎn)生的變形信息；

69、協(xié)同動(dòng)態(tài)耦合監(jiān)測(cè)模塊，用于分別部署電法系統(tǒng)和微震監(jiān)測(cè)系統(tǒng)進(jìn)行井下電阻率與微震聯(lián)合應(yīng)力變化監(jiān)測(cè)，實(shí)現(xiàn)井下采礦區(qū)大梯度異常體的圈定；結(jié)合井上采礦區(qū)沉降形變監(jiān)測(cè)信息，聯(lián)合井下電阻率與微震聯(lián)合應(yīng)力變化監(jiān)測(cè)，進(jìn)行協(xié)同監(jiān)測(cè)，完成對(duì)采空區(qū)形變過(guò)程協(xié)同動(dòng)態(tài)耦合監(jiān)測(cè)。

70、結(jié)合上述的所有技術(shù)方案，本發(fā)明所具備的有益效果為：針對(duì)某煤礦5-20303工作面的采空區(qū)，寬350m長(zhǎng)1000m范圍內(nèi)擬采用井下電阻率和微震聯(lián)合監(jiān)測(cè)，結(jié)合井上地表gnss，雷達(dá)衛(wèi)星insar技術(shù)，對(duì)采空區(qū)上三帶(彎曲下沉帶、裂隙帶以及冒落帶)形變過(guò)程協(xié)同動(dòng)態(tài)耦合監(jiān)測(cè)。發(fā)揮多手段的時(shí)空互補(bǔ)性和過(guò)程協(xié)同性，實(shí)現(xiàn)煤礦采空區(qū)井上井下形變動(dòng)態(tài)耦合監(jiān)測(cè)。

71、本發(fā)明根據(jù)需要，對(duì)某煤礦采空區(qū)，利用alos-2，radarsat-2，sentinel-1sar數(shù)據(jù)進(jìn)行連續(xù)4個(gè)月的insar采空區(qū)地表形變反演處理。此外，在現(xiàn)場(chǎng)合適位置布署北斗gnss，消除的大氣延遲誤差的影響，實(shí)現(xiàn)井上2種手段采空區(qū)形變高精度聯(lián)合監(jiān)測(cè)。為補(bǔ)充地表變形監(jiān)測(cè)方法的不足，分別部署電法系統(tǒng)和微震監(jiān)測(cè)系統(tǒng)進(jìn)行井下電阻率與微震聯(lián)合應(yīng)力變化監(jiān)測(cè)，實(shí)現(xiàn)井下采礦區(qū)大梯度異常體的圈定。結(jié)合井上井下4種方法協(xié)同監(jiān)測(cè)，完成對(duì)某采空區(qū)上三帶(彎曲下沉帶、裂隙帶以及冒落帶)形變過(guò)程協(xié)同動(dòng)態(tài)耦合監(jiān)測(cè)。