Calculation of maximum ground movement and deformation caused by mining

Transcription

1 Tas. Nofeous Me. Soc. Chia 1(011) s5-s59 Calculaio of maximum goud moveme ad defomaio caused by miig LI Pei xia 1,, TAN Zhi xiag 1,, DENG Ka zhog 1, 1. Key Laboaoy fo Lad Eviome ad Disase Moioig of Sae Bueau of Suveyig ad Mappig, Chia Uivesiy of Miig ad Techology, Xuzhou 111, Chia;. Jiagsu Key Laboaoy of Resouces ad Eviomeal Ifomaio Egieeig, Chia Uivesiy of Miig ad Techology, Xuzhou 111, Chia Received 19 Jue 011; acceped 10 Novembe 011 Absac: I ode o kow he maximum miig goud defomaio of abiay suface poi, diecios of he maximum suface il, cuvaue, hoizoal displaceme ad hoizoal defomaio caused by muli wokig faces wee deduced based o pobabiliy iegal mehod. The disibuio foms of suface defomaios i all diecio φ [0, π] wee obaied ad also equaios of maximum defomaios wee give based o calculaio of miig goud defomaio i diecio φ. A miig subsidece aalysis sysem was developed wih VB.0. The sysem implemes he pobabiliy iegal miig subsidece pedicio wih diec iegaio mehod, ad i ca avoid eos of abiay shaped wokig face subdivisio of adiioal mehod, ad accuacy of miig subsidece pedicio ca be impoved wih he diec iegal mehod. The sysem implemes he coou cha ad pofile cha ad also daa aalysis auomaio by maipulaio of he SURFER keel fucio ad i complemes he defec of exisig subsidece pedicio sofwae. Calculaio of maximum defomaios caused by muli wokig faces miig i a coal mie of Xuzhou, Chia was show as applicaio example. Egieeig applicaio idicaes ha maximum defomaio of miig goud suface ca be calculaed ad aalyzed by he sysem. The eseach povides a heoeical basis ad calculaio ool fo miig subsidece pedicio ad aalysis. Key wods: miig subsidece; pobabiliy iegal mehod; moveme ad defomaio; sysem developme 1 Ioducio Pobabiliy iegaio is he mos widely used mehod i miig subsidece pedicio [1] ad oe of he ecommeded mehods i he egulaio of he Naioal Coal Bueau []. The defomaio of suface moveme is a value elaed o is diecio; hey ae diffee i diffee diecios eve he same poi. I has impoa sigificace o safey miig, defomaio esisa sucue desig ad poecio i miig aea ad miig damage ideificaio o deemie he maximum defomaio ad mai (o maximum) defomaio diecio. A pese, he mos of he miig subsidece pedicio sofwae ca be used o calculae goud moveme ad defomaio alog a specified diecio caused by ecagula mie wokig faces. Bu hee ae also he followig defecs [3 11]. 1) Oly defomaio ad moveme alog specified diecios ca be calculaed, bu he maximum defomaio ad is diecio kow as he mos impoa values ca o be; o oly a sigle miig wokig face subsidece ca be calculaed, which gealy limis he accuacy of miig subsidece pedicio ad age of applicaio. ) Ifomaio of wokig faces ad goud moveme paamees eed o ee maually whe calculaio of miig subsidece affecively, i is complex o pepae he daa files. 3) Geeally, o ecagula wokig faces ae divided io ecagula o appoximae ecagula oes, i ivolves issues such as coodiae asfomaio ad small wokig faces coecio which iceases he complexiy of calculaio ad educes he accuacy. ) Mos pogams have gea defec o visualizaio. I ode o obai accuae goud moveme ad defomaio, maximum value of abiay poi miig Foudaio iem: Poec (107173) suppoed by he Naioal Naual Sciece Foudaio of Chia; Poec (CXZZ11_0300) suppoed by Scieific Iovaio Reseach Pogam fo College Posgaduaes i Jiagsu Povice Coespodig auho: LI Pei xia, Tel: ; E mail: pxlicum@1.com

2 LI Pei xia, e al/tas. Nofeous Me. Soc. Chia 1(011) s5 s59 subsidece caused by muli wokig face was deived accodig o basic piciples of pobabiliy iegal mehod; he accuae disibuio of suface moveme was aalyzed ad he pobabiliy iegal mehod miig subsidece pedicio was acualized wih a double iegal diecly. Miig subsidece pedicio ad aalysis sysem based o he pobabiliy iegal mehod was developed wih Acive X Auomaio echology i SURFER sofwae which combie he visualizaio daa pocessig ad esuls expessio ogehe. I ealized he visual expessio of miig subsidece calculaio ad esul aalysis. Fudameal heoy of pobabiliy iegal mehod Fudameal heoy of miig subsidece pedicio is pobabiliy iegal mehod based o he sochasic medium heoy. Accodig o he basic heoy of pobabiliy iegal mehod, he suface subsidece, goud il, cuvaue, hoizoal displaceme ad defomaio ca be calculaed wih double iegals as show i Eqs. (1 5). W ( x, y ) = W W ( x, y )d sd = = 1 D 0 e W 0 ( x s) + ( y ) (1) = 1 D W ( x, y) W ( x, y ) i( x, y, ϕ) = cosϕ + si ϕ = x y = 1 D π W 0 [( x s ) cos ϕ + ( x s) + ( y ) ( y )si ϕ ]exp π dsd i( x, y, ϕ) i( x, y, ϕ ) K ( x, y, ϕ) = cosϕ + si ϕ = x y {1 [( x s ) cos ϕ + π W 0 π = 1 D ( x s) + ( y ) ( y )si ϕ ] }exp π dsd π bw 0 U ( x, y, ϕ) = [( x s )cos 3 ϕ + = 1 D ( x s) + ( y ) ( y )si ϕ ]exp π dsd () (3) () π bw 0 π ε ( x, y, ϕ) = {1 [( x s ) cos ϕ + 3 = 1 D ( x s) + ( y ) ( y )si ϕ ] }exp π dsd s53 (5) whee W(x, y) is he miig subsidece of poi (x, y), is he umbe of wokig faces; W 0 is he maximum subsidece, W 0 =mqcosα; m is miig hickess; q is subsidece faco; α is dip agle of coal seam; D is he h miig aea, W e (x, y) is subsidece of ui miig of poi (x, y); is he mai ifluece adius of he h miig aea, =H /a β; H is miig deph of he h miig aea; a β is age of mai effec agle; b is displaceme faco; i(x, y, φ) is he suface il of poi (x, y) alog he diecio φ; K(x, y, φ) is he suface cuvaue; U(x, y, φ) is he suface hoizoal displaceme; ε(x, y, φ) is he suface defomaio. 3 Calculaio ad disibuio fom of maximum suface moveme ad defomaio caused by miig 3.1 Diecio of maximum defomaio Diecio of maximum suface il ad hoizoal displaceme defomaio Due o he popoioal of hoizoal displaceme ad suface il, hey have he same maximum defomaios diecio, ad he oly he suface il was show as example. Assume ha he diecio of maximum suface il (o hoizoal displaceme) of a abiay poi A(x, y) is φ i. Seekig paial deivaive of Eq. () ad le i equals o zeo, as show i Eq. (). i( x, y, ϕ ) ϕ ϕ= ϕ i = 0 Eq. (7) ca be obaied wih soluio of Eq. (). ϕ = aca i ( y )exp π dsd W 0 ( x s) + ( y ) = 1 D W 0 ( x s) + ( y ) = 1 D () ( x s)exp π dsd (7) whee φ i is he diecio of he maximum suface il (o maximum hoizoal displaceme). The maximum suface il of poi A ca be calculaed by puig he φ i of Eq. (7) io he Eq. (). I he simila way, he maximum hoizoal displaceme of poi A ca be calculaed by puig i io he Eq. () Diecio of maximum suface cuvaue ad hoizoal defomaio Due o he popoioal of suface cuvaue ad

3 s5 LI Pei xia, e al/tas. Nofeous Me. Soc. Chia 1(011) s5 s59 hoizoal defomaio, he maximum defomaios have he same diecio, ad he oly he cuvaue was show as example. Assume ha he diecio of maximum cuvaue (o hoizoal defomaio) of a abiay poi A(x, y) is φ K. Seekig paial deivaive of Eq. (3) ad le i equals o zeo, as show i Eq. (8). K( x, y, ϕ ) ϕ ϕ = ϕ k = 0 Eq. (9) ca be obaied wih soluio of Eq. (8). 1 W aca 0 ϕ k = ( x s )( y ) = 1 D ( x s) + ( y ) exp π dsd / = 1 D W 0 [( x s) ( y ) ] ( x s) + ( y ) exp π dsd (8) ( 9 ) whee φ K is he diecio of he maximum suface cuvaue (o maximum hoizoal defomaio). The maximum cuvaue of poi A ca be calculaed by puig φ K of Eq. (9) io he Eq. (3). I he simila way, he maximum hoizoal defomaio of poi A ca be calculaed by puig i io Eq. (5). 3. Calculaio of maximum suface moveme ad defomaio I chape 3.1, hough supeposiio calculaio of muliple wokig faces, he diecios of maximum goud moveme ad defomaio of abiay poi caused by muli faces miig wee give usig highe mahemaics kowledge. I pocess of poecio of miig subsidece, defomaio esisa sucue desig ad poecio i miig aea ad miig damage ideificaio, he maximum value of defomaio is he maily coce value bu o is diecio. Accodig o fudameal heoy of pobabiliy iegal mehod, some meaigful coclusios wee go by fuhe calculaio of Eqs. ( 5) Calculaio of suface il ad hoizoal displaceme Fo he abiay poi A(x, y) i he affeced aea of muliple wokig face miig subsidece, is coodiae is deemied. The il of poi A alog he abiay diecio φ ca be expessed as Eq. (10). 1 i( x, y, ϕ) = π W0 [( x s ) cos ϕ + = 1 D i ( x s) + ( y ) ( y )si ϕ ]exp π dsd π W 0 i( ϕ) = cos ϕ ( x s)dsd + If = 1 D π W 0 y s (10) = 1 D si ϕ ( )d d π W 0 ( )d d, = 1 D A = x s s = 1 D π W 0 B = ( y )dsd, ad he Eq. (10) ca be expessed as Eq. (11). i( ϕ) = Acosϕ + B si ϕ = A + B cos θ = si θ = A B ( cosϕ + si ϕ ) A + B A + B Agle θ mees he equiemes of Eqs. (1 13). A A A + B B + B Eq. (10) ca be expessed as: i( ϕ) = Acosϕ + B si ϕ = A + B (cosθ cosϕ + siθ si ϕ ) = (11) (1) (13) A + B cos( ϕ θ ) (1) As kow i Eq. (1), age of suface il of abiay poi alog he abiay diecio is bewee [ A + B, A + B ] ; hee ae wo equal ad opposie maximum suface il A + B. The values associaed wih he locaio of he poi oly bu ohig elaioship wih diecios. The maximum value of hoizoal displaceme is simila, ad i which he maximum hoizoal displaceme of A(x, y) is show i Eq. (15). = 1 D π bw 0 ( x s)dsd + 3 π bw 0 ( x s) ( y ) 3 1 D = + ( y )exp π dsd (15) The values associaed wih he locaio of he poi oly ad have ohig elaioship wih diecios. Wih calculaio mehod above, if he diecio of 1/

4 LI Pei xia, e al/tas. Nofeous Me. Soc. Chia 1(011) s5 s59 he maximum suface il (o hoizoal displaceme) eeds o o kow, maximum value of suface il (o hoizoal displaceme) ca be calculaed diecly. Take φ [0, π] i Eq. (1), disibuio figue ca be go as show i Fig. 1 (dashed deoes he egaive values). The suface il (o hoizoal displaceme) of abiay poi i miig affecig aea all has he disibuio fom show as Fig. 1. = 1 D π W 0 C = ( x s)( y ) ( x s) + ( y ), The, Eq. (1) ca be chaged io: K ( ϕ) = A B( cos ϕ 1) C siϕ cos ϕ = s55 Fig. 1 Disibuio schemaic diagam of suface il (o hoizoal displaceme) 3.. Calculaio of suface cuvaue ad hoizoal defomaio I a simila way, accodig o Eq. (3), he cuvaue of poi A(x, y) alog he diecio φ ca be expessed as Eq. (1). π W 0 π = 1 D K ( ϕ ) = {1 [( x s) + ( y ) ]} Le ( x s) + ( y ) m π W 0 ϕ x s y = 1 D 1 (cos 1) [( ) ( ) ] ( x s) + ( y ) siϕ cos ϕ = 1 D π W 0 ( x s)( y ) ( x s) + ( y ) π W 0 π = 1 D A = {1 [( x s) + ( y ) ]} ( x s) + ( y ), = 1 D 1 π W 0 B = [( x s) ( y ) ] ( x s) + ( y ), (1) A B cos( ϕ) C si( ϕ ) = B A B + C cos ϕ + B + C cos β = si β = C si ϕ (17) B + C Agle β mees he equieme of Eqs. (18 19). B B + C C B + C The Eq. (17) ca be expessed as: (18) (19) K ( ϕ) = A B + C cos( ϕ β ) (0) Le φ [0, π], hee ae five siuaios of he disibuio foms of cuvaue of poi A. Fisly, if A > B + C, he K(φ)>0 i age φ [0, π]. I his case, he suface cuvaues value of poi A is a posiive value i all diecio, i ages [ A B + C, A + B + C ], ad he maximum value of cuvaue is A + B + C. Accodig o similaiy, he hoizoal displaceme is a posiive value, ad poi A subec o esile defomaio i all diecio. Ude he cicumsaces disibuio of cuvaue (o hoizoal displaceme) is show as Fig. (a). Secodly, if A < B + C, he K(φ)>0 i age φ [0, π]. I his case, he suface cuvaues value of poi A is a egaive value i all diecio, i ages [ A B + C, A + B + C ], ad he maximum value of cuvaue is A B + C ; i is a egaive value. Accodig o similaiy, he hoizoal displaceme is a egaive value, ad poi A subec o compessio defomaio i all diecio. Ude he cicumsaces disibuio of cuvaue (o hoizoal displaceme) is show as Fig. (b). Thidly, if B + C < A < B + C, he he age of K(φ)>0 is [ A B +C, A+ B +C ] whe φ [0, π]. The suface cuvaue has wo ohogoal picipal diecios, he cuvaue value ae opposie sig

5 s5 LI Pei xia, e al/tas. Nofeous Me. Soc. Chia 1(011) s5 s59 alog wo mai diecios. Thee is a posiive cuvaue value alog a piciple diecio ad a egaive oe alog aohe. Similaiy, he hoizoal displaceme of poi A subecs o esile defomaio i oe piciple diecio ad compessive defomaio i aohe oe. Ude he cicumsaces disibuio of cuvaue (o hoizoal displaceme) is show as Fig (c). Fouhly, if A = B + C, he K(φ)>0 i age φ [0, π]. I his case, he suface cuvaues value of poi A is a posiive value i all diecio, i ages [0, A + B + C ], ad he maximum value of cuvaue is A + B + C. Ude he cicumsaces, he suface cuvaue value is A + B + C i a piciple diecio ad he value is zeo i aohe piciple diecio. Accodig o similaiy, he hoizoal displaceme is a opposie value. Poi A subecs o compessio defomaio i all diecios. Ude he cicumsaces disibuio of cuvaue (o hoizoal displaceme) is show as Fig. (d). Fifhly, if A = B + C, he K(φ)>0i age φ [0, π]. I his case, he suface cuvaues value of poi A is a egaive value i all diecios, i ages [ A B + C,0], ad he maximum value of cuvaue is A B + C. Ude he cicumsaces, he suface cuvaue value is A B + C i a piciple diecio ad he value is zeo i aohe piciple diecio. Accodig o similaiy, he hoizoal displaceme is a egaive value. Poi A subecs o compessio defomaio i all diecios. Ude he cicumsaces disibuio of cuvaue (o hoizoal displaceme) is show as Fig. (e). Ay whee i egio of miig subsidece, disibuios foms of cuvaue (o hoizoal defomaio) alog φ [0, π] ca oly appea oe of five foms show i Fig.. Developme of miig subsidece pedicio ad aalysis sysem.1 Composiio ad sucue of sysem The sysem was developed accodig o basic piciples of pobabiliy iegal mehod ad maximum value calculaio mehod based o he sucue ad modula desig philosophy. Is oveall implemeaio is a secoday developme based o AUTOCAD ad SURFER, which miig subsidece pedicio pa is based o AUTOCAD ad miig subsidece aalysis ad visualizaio pa based o he SURFER. The oveall idea is o pepae basic daa of miig subsidece pedicio by ey miig ifomaio hough daa maageme secio, file geeaio, daabase maageme; ad he goud moveme ad defomaio ca be pediced by miig subsidece pedicio secio; a las, saisical aalysis ad visualizaio expessio ca be implemeed by maipulaig he keel pogam of SURFER wih daa pediced. The flow cha of abiay shaped muliple wokig faces miig subsidece pedicio ad aalysis sysem is show i Fig. 3.. Implemeaio of abiay shaped muliple wokig faces miig subsidece pedicio Geeally, abiay shaped wokig faces should be divided io muliple ecagula faces o calculae miig subsidece. Subdivisio mehod is simple o impleme, bu i eeds o modify he pedicio paamees fo small ecagula wokig faces divided, ad calculaio esuls is appoximae fo subdivisio; ad i icease he hadess o pepae he daa. As kow i Eqs. (1 5), fo a abiay shaped miig aea D, miig subsidece pedicio ca be implemeed wih a double iegal. I ode o obai uppe ad lowe limi of iegal, miig egio D should subdivide io umbe of iagles ad pobabiliy iegal mehod miig subsidece pedicio ca be implemeed by double iegal diecly. Tiagulaio mehod ca divide abiay shaped wokig faces accuaely ad avoid he defecs of ecagle subdivisio. Basic flow cha of miig subsidece pedicio calculaio is show i Fig.. Fig. Disibuio schemaic diagam of suface cuvaue (o hoizoal defomaio).3 Miig subsidece aalysis usig SURFER The miig subsidece aalysis sysem used Acive X Auomaio echology ad he pogam ieface o cusomize he SURFER. SURFER ca be used o daw he coou, pase diagams, images ad veco gaphic easily, ad also be used o hee dimesio space aalysis such as ed suface aalysis, volume ad aea calculaio, slope aalysis, pofile aalysis e al. 0

6 LI Pei xia, e al/tas. Nofeous Me. Soc. Chia 1(011) s5 s59 s57 Fig. 3 Geeal flow cha of abiay shape muliple wokig faces miig subsidece pedicio ad aalysis sysem diffee kids of Acive obecs embedded i SURFER ca be called by exeal applicaios, ad almos all fucios ca be accomplished by cusomize developme [11 15]. Suface moveme ad defomaio calculaed by miig subsidece module ca be used o make vaiey maps ad daa aalysis wih keel fucio i SURFER. Seps ae as follows: 1) Thee dimesioal suface moveme ad defomaio daa should be go fisly by miig subsidece pedicio module, poi coodiae ad defomaio value ae icluded. ) SURFER coecio esablished, ad daa pediced should be gid pocessed by SURFER keel daa pocessig fucio. 3) Wih expessio fucio of SURFER gaphic, vaious maps ad daa ca be made, ad do he ovelay aalysis wih wokig faces. ) Vaiey of saisical aalysis ad daa pocessig ca be doe wih SURFER daa aalysis fucio. Fig. Flow cha of miig subsidece pedicio. Sysem fucios ad feaues Compaed wih simila miig subsidece pedicio sysem i Chia, he miig subsidece sysem has he followig feaues. 1) Goud subsidece, il, cuvaue, hoizoal displaceme ad hoizoal defomaio caused by muliple abiay shaped wokig faces ca be calculaed wih he sysem; also he diecios of

7 s58 LI Pei xia, e al/tas. Nofeous Me. Soc. Chia 1(011) s5 s59 maximum defomaio ad is values ca be go, i s he mai iovaios i he developme of his sysem. ) Wokig faces daa maageme secio. The sysem ca exac basic ifomaio of wokig faces i AUTOCAD; basic daa ca be acquied ad maaged by he sysem diecly which impove he efficiecy. 3) Esablishme of miig subsidece paamees daabase maageme pogam. Miig subsidece daa ca be maaged ad used effecively wih i. ) Goud moveme ad defomaio caused by abiay shaped muliple wokig faces ca be pediced by he miig pedicio sysem. Goud subsidece, il, cuvaue, hoizoal displaceme, hoizoal defomaio ad is maximum values ca be calculaed by opeaios of he sysem. The pedicio esuls ca be oupu i he daa files. Goud defomaio caused by abiay shaped ad abiay umbes of wokig faces ca be go wih he sysem. 5) Miig subsidece aalysis sysem. The miig subsidece aalysis sysem is a subpogam of miig subsidece pedicio developed by secodly developme wih Acive X Auomaio echology i SURFER. Miig subsidece daa pediced ca be pocessed o visualizaio expessio, map poducio, ad auomaio daa aalysis hough he opeaio of he sysem; also base suface map of wokig faces, moveme ad defomaio coou, pofile, hee dimesioal sufaces, ovelay cha ad elaed aalysis ca be poduced. pedicio ad aalysis sysem descibed above. As space limied, oly coou of maximum goud il ad cuvaue caused by No.701 ad No.731 wok faces ae show i Fig. 5 ad Fig.. Fig. 5 Coou cha of maximum suface il (mm/m) Fig. Coou cha of maximum suface cuvaue (mm/m ) 5 Applicaio ad discussio 5.1 Egieeig applicaio I ode o mie coal pilla esouces ude he villages i a coal mie i Xuzhou, Chia, wo wokig faces No.701 ad No.731 wee plaed o mie. Sufe above he wokig faces is fla, aveage elevaio is 3 37 m, desig legh of No.713 is m, desig widh m, miig he 7h coal seam, aveage deep is 597 m, aveage hickess of coal seam is 7 m, seam dip agle is 0 1 ad aveage agle is. No.701 wokig face is locaed i oh of No.713 face; desig legh is 5 m, desig widh m, miig he 7h coal seam oo, aveage miig deph is 50 m, aveage hickess of coal seam is 7 m, ad aveage coal seam dip agle is 8. Buildigs i ea villages would be damaged afe wo wok faces miig. I ode o kow he goud moveme ad defomaio afe miig ad povide a basis fo mie decisio makig. Miig subsidece ad goud defomaio such as subsidece, il, cuvaue, hoizoal displaceme, hoizoal defomaio values alog a specified diecio also he maximum goud defomaios wee calculaed by he miig subsidece Maximum values ad is diecio of abiay poi goud moveme ad defomaio ca be calculaed accuaely;compae wih adiioal mehod, i aalysis he goud defomaio wih he maximum values of all diecios isead of alog he specified diecio, much moe accuae buildig damage degee ad damage age ca be go wih he ew mehod above ad i povides a much moe accuae piciple fo goud moveme ad buildigs damage aalysis; he mehod has impoa pacical sigificaces o impove miig subsidece pedicio aalysis, poec he goud buildig ad coal mie safey. 5. Discussio Pobabiliy iegal mehod is a mos widely used mehod i Chia miig subsidece pedicio. No ecagula wokig faces geeally should be divided io ecagula o appoximae ecagula oes adiioally; i ivolves issues such as coodiae asfomaio ad small wokig faces coecio ad iceases he complexiy of calculaio ad educes he calculaio accuacy. O he ohe had, fo muliple wok faces miig, i will uable o deemie he maximum defomaio diecio ad is value afe

8 LI Pei xia, e al/tas. Nofeous Me. Soc. Chia 1(011) s5 s59 complex coodiae asfomaios, which limis he aalysis accuacy of pobabiliy iegal mehod gealy. Wok faces pofile eos ca be educed by usig diecly iegal mehod o impleme he pobabiliy iegal mehod; miig subsidece pedicio accuacy ca be impoved. All calculaio is compleed i he same coodiae wih he diecly iegal mehod. I does o eed coodiaes asfomaios, ad i educes he complex of he calculaio; maximum value of goud defomaio ad is diecio ca be deemied wih he mehod, ad i ca impoves he accuacy of he miig subsidece pedicio gealy. Coclusios 1) Wih summay of poblems i cue miig subsidece pedicio mehods, he maximum goud defomaio diecios ad is values ae eseached ad is calculaio mehods ae go based o pobabiliy iegal mehod. ) Rages ad disibuios siuaio of abiay poi goud moveme ad defomaio alog he diecio [0, π) ae aalyzed wih igoous mahemaical deivaio of pobabiliy iegal mehod. I shows ha abiay shaped muliple wok faces miig, abiay poi i he miig subsidece affecig aea, have maximum il ad hoizoal displaceme, he wo maximum values ae equal ad opposie i diecio; he maximum cuvaue ad hoizoal defomaio occu i fou diecios pepedicula o each ohe, ad hee ae five disibuio foms of cuvaue ad hoizoal displaceme i diecio [0, π). 3) Abiay shaped muliple wok faces miig subsidece pedicio ad aalysis sysem is developed. The sysem implemes he pobabiliy iegal mehod wih diec iegal mehod calculaio, avoidig he eos of ecagula subdivisio of o ecagula, ad impoved he accuacy of miig subsidece pedicio. Miig subsidece visualizaio, map poducio ad daa aalysis visualizaio ae implemeed by igh couplig of VB ad SURFER based o Acive X echology. Egieeig applicaio shows ha he sysem ca calculae ad aalyze he maximum value of goud defomaio; i has impoa pacical sigificace o impove he accuacy of miig subsidece. Refeeces s59 [1] PENG S S. Suface subsidece egieeig [M]. Coloado: Sociey fo Miig, Meallugy ad Exploaio, Ic, 199: [] Fomulaed by Naioal Coal Bueau of P.R.C. Regulaio of miig ad pilla leavig ude buildig, wae body, ailway ad mai udegoud egiee [M]. Beiig: Coal Idusy Pess, 000: (i Chiese) [3] WU Ka, ZHOU Mig. Miig subsidece pedicio sysem [M]. Xuzhou: Chia Uivesiy of Miig ad Techology Pess, 1999: 1 3. (i Chiese) [] HABIB A. Back of evelope miig subsidece esimaio [J]. Ausalia Geomechaics Joual, 009, (1): 9 3. [5] TANG Fu qua. Miig subsidece moioig usig he mehod of combiig ISAR ad GPS echology [J]. Joual of Coal Sciece ad Egieeig, 011, 17(): [] HUANG Le ig. Reseach o defomaio velociy i dyamic suface subsidece defomaio [J]. Tasacios of Nofeous Meals Sociey of Chia, 005, 15(S1): [7] DIAZ FERNANDEZ M E, ALVAREZ FERNANDEZ M I, ALVAREZ VIGIL A E. Compuaio of ifluece fucios fo auomaic miig subsidece pedicio [J]. Compuaioal Geoscieces, 010, 1(1): [8] TRCKOVA J. Expeimeal 3 D modelig of suface subsidece affeced by udegoud miig aciviies [J]. Joual of he Souh Afica Isiue of Miig ad Meallugy, 009, 109(1): [9] LIIVEIER D, HARLALKA A. Numeical sudy of he soil sucue ieacio wihi miig subsidece aeas [J]. Compues ad Geoechics, 010, 37(): [10] HAN Bao mi, OU Ji ku, CHAI Ya u, LU Xiu sha. Mehod fo pocessig daa obseved fom GPS fo subsidece suveyig i miig aea [J]. The Chiese Joual of Nofeous Meals, 1999, 9(): (i Chiese) [11] DING De xi, ZHANG Zhi u, BI Zhog wei. New appoach o pedicig miig iduced suface subsidece [J]. Joual of Ceal Souh Uivesiy of Techology, 00, 13(): 38. [1] SI Zheg, LI Shou yi, HUANG Lig zhi, CHEN Yao log. Visualizaio pogammig fo bach pocessig of coou maps based o VB ad sufe sofwae [J]. Advaces i Egieeig Sofwae, 010, 1(7 8): [13] MA Yu feg, GAO Chu xiag. Auomaic dawig figues of meeoological eleme by combied pogammig fo VB ad sufe [J]. Meeoology Joual of Ie Mogolia, 010, 010(): (i Chiese) [1] LUO Liag, ZENG Tao, LIANG Feg, Applicaio of ue hee dimesio display o subsidece foecas i miig aea [J]. Suvey, 009, 3(5): 0 0. (i Chiese) [15] TANG Fu qua, WANG Gui sheg, DAI Ju peg. Applicaio of Sufe8.0 i visualizaio expessio ad mappig of miig subsidece daa [J]. Idusy ad Mie Auomaio, 010, 010(10): 7 9. (i Chiese) (Edied by ZHAO Ju)