【正文】 ....... 46 結(jié)論 ..................................................................... 51 參考文獻(xiàn)： ............................................................... 52 致 謝 ................................................................... 53 安徽理工大學(xué)畢業(yè)論文 1 緒論 開采沉陷的研究對國民經(jīng)濟(jì)的發(fā)展和人民生活水平的提高有著重要的意義，然而隨著開采沉陷研究的規(guī)模越大，有關(guān)開采沉陷的信息量也相應(yīng)的增加，對這些信息的存儲和管理也就相當(dāng)繁瑣。 生成的曲線圖更 有利于直觀現(xiàn)時地分析開采對地表的影響。 本科畢業(yè)論文 基于 CAD 和 VBA 的地表移動觀測 站 數(shù)據(jù)處理 系統(tǒng) DATA PROCESSING SYSTEM ON MINING AREA SURFACE MOVEMENT STATION BASED ON AUTOCAD VBA 學(xué)院（部）： 測繪學(xué)院 專業(yè)班級： 測繪 091 班 學(xué)生姓名： 劉石磊 指導(dǎo)教師： 王磊講師 20xx 年 5 月 31 日 安徽理工大學(xué)畢業(yè)論文 II 基于 CAD 和 VBA 的地表移動觀測 站 數(shù)據(jù)處理 系統(tǒng) 摘 要 采用 AutoCAD VBA 二次開發(fā)技術(shù)，實(shí)現(xiàn)了地表移動觀測數(shù)據(jù) 處理、成圖自動一體化 。且圖中按原比例繪出對應(yīng)地下煤層開采情況，可很方便地從圖中直接量取各角量參數(shù)，距離參數(shù)。 surface movement and deformation。如何有效地對開采沉陷信息管理和利用起來，如何使得開采沉陷的信息能夠在更大范圍內(nèi)安全地進(jìn)行傳輸、共享和管理，一直是很多企業(yè)單位急需解決的問題。 地下煤層被開采出來以后，巖體內(nèi)部形成一個采空區(qū)，導(dǎo)致周圍巖體的原有應(yīng)力平衡狀態(tài)受到破壞，引起應(yīng)力的重新分布，使巖體產(chǎn)生移動和破壞，直到達(dá)到新的平衡。 建立地表移動觀測站是研究開采沉陷和建筑物破壞規(guī)律的最可靠手段 。目前是最直接有效的方法安徽理工大學(xué)畢業(yè)論文 2 是建立地表移動觀測站，定期對觀測站上的監(jiān)測點(diǎn)及控制點(diǎn)進(jìn)行觀測。在輸出下沉、傾斜、曲率、水平移動、水平變形計(jì)算結(jié)果的 同時，將其各曲線以一定比例在 CAD 中成圖，可從曲線上量取任意點(diǎn)的對應(yīng)值。需要求取的地表移動參數(shù)，可在圖中直接精確量取和標(biāo)注。 1885年法國礦長 總結(jié)了開采沉陷時間的矛盾意見 [3]。本世紀(jì) 20 至 30年代 ,斯奇米茨、凱因斯特和巴斯研究了開采影響的作用面積及分帶 ,可作為影響函數(shù)法的初級階段 。學(xué)者 改進(jìn)了沉陷觀測的方法，并對開采沉陷的的理論方法也做出了巨大貢獻(xiàn) [5]。在西方發(fā)達(dá)國家 ,利用信息技術(shù)開發(fā)出了一批適用于礦山的應(yīng)用軟件如加拿大的 LYNX 系統(tǒng) ,用于地質(zhì)學(xué)研究從勘探資料分析地質(zhì)特征 。 隨著計(jì)算機(jī)和 3S 技術(shù)的廣泛應(yīng)用，將新技術(shù)應(yīng)用于開采沉陷觀測和數(shù)據(jù)處理也成為國際上普遍研究的的一個重要方面。 國內(nèi)研究現(xiàn)狀 我國是以煤炭為能源的國家，自新中國成立以來就開始開采沉陷方面的研究， 80、90年代至今，我國開采沉陷理論和實(shí)踐研究出現(xiàn)了日新月異的發(fā)展。王金莊等建立了用于預(yù)計(jì)主斷面地表下沉的典型曲線法。李增琪建立了巖層移動的三維層狀模型。吳立新、王金莊提出了條帶開采覆巖破壞的托板理論。 陳宜金，黃紹東對開采沉陷信息處理軟件系統(tǒng)各模塊間關(guān)系和設(shè)計(jì)原則進(jìn)行了探討 [12]。 郝長勝等研究 者利用 VB 作為系統(tǒng)開發(fā)前臺，基于 GIS 桌面軟件 MapInfo及 MapX 組件來進(jìn)行系統(tǒng)二次開發(fā)設(shè)計(jì)：根據(jù)模塊各功能應(yīng)盡量共享、數(shù)據(jù)相互調(diào)用、減少資源內(nèi)存儲量等原則，該系統(tǒng)煤礦開采沉陷預(yù)計(jì)分析整體流程分為數(shù)據(jù)管理、變形預(yù)計(jì)參數(shù)計(jì)算、圖形輸出、三維顯示及預(yù)計(jì)結(jié)果分析五大模塊。例如賈小敏、余學(xué)祥等討論了應(yīng)用Visual Basic 和 Auto CAD 二次開發(fā)技術(shù)進(jìn)行礦區(qū)地表移動變形值計(jì)算和變形曲線繪制的方法，實(shí)現(xiàn)了在此環(huán)境下變形值的計(jì)算和曲線的 研究內(nèi)容 VBA ( visual basic for application) 是由 Microsoft 創(chuàng)建的，用來自動執(zhí)行任務(wù)的一個編程環(huán)境，它提供了一些用來創(chuàng)建圖形用戶界面的可拖拉工具和用來與AutoCAD 對象交換的編程語言。需要求取的地表移動參數(shù)，可在圖中直接精確量取和標(biāo)注 擬采用的研究思路與方法 研究思路與方法 以礦區(qū)實(shí)測數(shù)據(jù)為例 ,應(yīng)用 VBA 與 AutoCAD 二次開發(fā)技術(shù)進(jìn)行礦區(qū)地表移動變形值安徽理工大學(xué)畢業(yè)論文 5 的解 算 , 分析繪制地表移動變形曲線圖的 理論和方法 , 并以 VBA 為基礎(chǔ) , 實(shí)現(xiàn)了礦區(qū)地表移動變形值的快速求解 , 從而達(dá)到了對地下開采造成的地表移動變形的直觀性描述。 N為觀測線上的一觀測點(diǎn)， N’ 為 N受地下開采影響后的點(diǎn)。方位角 α12 、 α ， S 的求取在此就不給出了。cosβ， β=α12—α 同理可求出 n’ 的坐標(biāo)。移動和變形的計(jì)算主要包括：各測點(diǎn)的下沉和水平移動，相鄰兩測點(diǎn)間的傾斜和水平變形，相鄰三點(diǎn)的曲率變形 [24]。 水平移動 下沉盆地內(nèi)，某點(diǎn)沿某一水平方向的位移叫水平移動，用 U 表示，單位 mm，則 m次觀測 n 點(diǎn)的水平移動為 Un= Ln m — Ln 0 （ 26） 式中， Ln m， Ln 0——分別為 m 次觀測時和首次觀測時 n 號點(diǎn)至觀測線控制點(diǎn)間的水平距離，直接用坐標(biāo)求得。其最大特點(diǎn)在于兩個方面：第一，與 VB 有著幾乎相同的開發(fā)環(huán)境和語法，具備功能強(qiáng)大和易于掌握的特點(diǎn)；第二，在于它的 for Application功能，針對性強(qiáng)，它駐留在主程序的內(nèi)部，使其結(jié)構(gòu)精簡，且代碼運(yùn)行效率高，使語言的引擎在技術(shù)上與開發(fā)環(huán)境分離 [6]。 確定煤層采動信息 為反應(yīng)地下開采與地表移動和變形規(guī)律，在進(jìn)行各項(xiàng)計(jì)算之前應(yīng)當(dāng)首先確定觀測數(shù)據(jù)所對應(yīng)的地下開采情況。讀入基點(diǎn)坐標(biāo)，計(jì)算規(guī)劃方向 JX1 = Val() JY1 = Val() JX2 = Val() JY2 = Val() a = Atn((JY2 JY1) / (JX2 JX1)) 安徽理工大學(xué)畢業(yè)論文 11 讀入 觀測數(shù)據(jù)進(jìn)行坐標(biāo)規(guī)劃 數(shù)據(jù)文件編輯成 txt 文檔，首次觀測平面坐標(biāo)及高程在前，第 m次觀測在后。 Chr(13) amp。 煤層繪制的部分代碼如下： 1) 煤層走向信息繪制 Private Sub OptionButton1_Click() Dim layer1 As AcadLayer 安徽理工大學(xué)畢業(yè)論文 13 Set layer1 = (坐標(biāo) ) = 255 = layer1 Call MC Call CKQTC End Sub 2）煤層傾向信息繪制 Private Sub OptionButton2_Click() Dim layer1 As AcadLayer Set layer1 = (坐標(biāo) ) = 255 = layer1 Call mcq Call QXTC End Sub Public Sub MC() Dim HZ, ZD, KX, KY, JX1, JY1, SKQ, CH 39。煤層中間點(diǎn)采深，平均采厚 39。 在 模型空間中創(chuàng)建實(shí)體對象 Set solidObj = (point1, point2, point3, point4) ZoomAll End Sub 煤層傾向代碼如下： Public Sub mcq() Dim HZ, D, QBX, QBY, JX1, JY1, Q, CH, j 39。煤層中間點(diǎn)采深，平均采厚 39。 在模型空間中創(chuàng)建實(shí)體對象 Set solidObj = (point1, point2, point3, point4) ZoomAll End Sub 3）曲線繪制代碼 ： Private Sub CommandButton5_Click() Dim layer1 As AcadLayer Set layer1 = (下沉曲線 ) = acRed = layer1 Call ZS 安徽理工大學(xué)畢業(yè)論文 18 Dim wn() As Double Dim AQ() As Double Dim ds() As Double Dim gx() As Double Dim gy() As Double Dim n As Double n = Val() JX1 = Val() JY1 = Val() JX2 = Val() JY2 = Val() a = Atn((JY2 JY1) / (JX2 JX1)) Dim ss As String ss = Set fso = CreateObject() Set ttfile = (ss, 1) read = read = Replace(read, , ) Do While InStr(1, read, String(2, vbTab)) 0 read = Replace(read, String(2, vbTab), vbTab) Loop 39。 ds(i) amp。 AQ(3 * k) amp。 AQ(3 * k + 2) amp。 wi(j) amp。 AQ(3 * k + 1) amp。MsgBox read brr = Split(read, vbCrLf) ReDim ws(UBound(brr)) ReDim AQ(3 * UBound(brr) 1) ReDim ds(UBound(brr)) ReDim hs(UBound(brr)) ReDim gx(UBound(brr)) ReDim gy(UBound(brr)) ReDim hx(UBound(brr)) ReDim hy(UBound(brr)) ReDim ws(UBound(brr)) For i = LBound(brr) To UBound(brr) If brr(i) Then arr = Split(brr(i), vbTab) If UBound(arr) 5 Then dn = arr(0) gx(i) = Format(Sqr((CDbl(arr(1)) JX1) ^ 2 + (CDbl(arr(2)) JY1) ^ 2) * Cos(a) + JX1, ) gy(i) = Format(Sqr((CDbl(arr(1)) JX1) ^ 2 + (CDbl(arr(2)) JY1) ^ 2) * Sin(a) + JY1, ) ds(i) = Format(Sqr((gx(i) JX1) ^ 2 + (gy(i) JY1) ^ 2), ) hx(i) = Format(Sqr((CDbl(arr(4)) JX1) ^ 2 + (CDbl(arr(5)) JY1) ^ 2) * Cos(a) + JX1, ) 安徽理工大學(xué)畢業(yè)論文 27 hy(i) = Format(Sqr((CDbl(arr(4)) JX1) ^ 2 + (CDbl(arr(5)) JY1) ^ 2) * Sin(a) + JY1, ) hs(i) = Format(Sqr((hx(i) JX1) ^ 2 + (hy(i) JY1) ^ 2), ) ws(i) = hs(i) ds(i) End If End If Next For k = 0 To UBound(brr) 1 AQ(3 * k) = ds(k): AQ(3 * k + 1) = ws(k) / n: AQ(3 * k + 2) = 0 mystr = mystr amp。 , amp。mystr = mystr amp。 , amp。 vbCrLf Next = mystr Dim splineObj As AcadSpline Dim startTan(0 To 2) As Double Dim endTan(0 To 2) As Double startTan(0) = 0: startTan(1) = 0: startTan(2) = 0 endTan(0) = 1500: endTan(1) = 0: endTan(2) = 0 Set splineObj = (AQ, startTan, endTan) End Sub 之后，將 單位一致的 變形 量繪在同一坐標(biāo)下 ，