賈星存
(山西西山晉興能源有限責任公司 斜溝煤礦,山西 呂梁 033602)
近年來,成熟的礦井通風輔助決策型軟件已經幫助不少礦井完成了優化改造[1,2]. 借助計算機強大的計算能力,再復雜的風網也能夠快速完成解算任務[3,4]. 這使得人們能夠在計算機上對通風系統調節方案進行“預演”,找到最優的優化方案[5,6]. 本文運用Ventsim三維通風仿真軟件,對斜溝煤礦15采區進行了三維通風仿真模型的建立和優化調節方案的仿真設計.
斜溝煤礦15采區通風由斜溝回風井負擔,由副平硐、主帶式輸送機井進風,斜溝回風井同時給+700水平提供約1 000 m3/min 風量。15采區共有18503工作面和16502工作面兩個生產工作面;有18505材料巷和18505帶式輸送機巷兩個煤巷掘進生產工作面;其他用風地點有井下火藥庫、水倉水泵房、15采區變電所、15采區消防材料庫。
優化前斜溝回風井主扇負壓:1 600 Pa,風葉角度:-3°/-3°,排風量:8 460 m3/min,采區目前通風阻力分布合理,符合規程要求,但考慮到日后加大產量,需降低阻力,提高風量,以保障生產需求,需要對15采區通風系統進行優化改造。
斯考得亨斯雷算法是通風網絡解算的通用算法,經過大量實踐驗證,穩定高效,且能夠獲得高精度解[7-9]. 求解的基本方程有:
1) 質量守恒。
網絡中流入和流出節點的質量守恒,網絡m個節點,可建立m-1個方程。
(1)
式中:
bij—基本關聯矩陣的元素;
qi—第i條分支的流量。
式(1)用矩陣表達為:
BQ=0
(2)
式中:
B—基本關聯矩陣;
Q—分支流量列矩陣。
2) 能量守恒。
對于一條巷道,根據伯努利方程可得:
(3)
hi=ΔP+ΔD+ΔW=Rq2
(4)
對任一個穩態回路:
∑hi-hz-hf=0
(5)
對所有的回路,表示成矩陣形式為:
CH=0
(6)
風網有n條分支,m個節點,這樣又得到n-m+1個方程。式中:C為基本回路矩陣;H為分支阻力列矩陣;R為巷道的風阻;hi、hz、hf分別為分支阻力、自然風壓和機械風壓。
根據質量守恒和能量守恒共n個方程n個未知量,可以聯立這些方程求解。聯立上述方程得到一組非線性方程組:
Fn(Qn)=0
(7)
因為是非線性方程組,用迭代法求解:
(8)
(9)
(10)
(11)
這樣便得到了一次迭代修正后的風量值,經過反復迭代運算,ΔQ小到滿足精度的要求,便可完成解算,此時的Q便是需要的解。
具體到通風網絡解算中,斯考德亨斯雷解法是將J矩陣簡化為主對角矩陣,見式(11),以簡化矩陣求逆操作,減少計算量,但同時其收斂效果受到影響,這種解法需要將解算的量進行適當排序方可有較好的收斂效果,但同時這種解法對初始值的要求放寬了。
Ventsim三維通風仿真系統提供了一個用于分析風流模擬、熱模擬、污染物模擬和通風經濟性的集成工具箱。Ventsim能夠模擬與記錄生產礦井風流及風壓,模擬新掘和廢棄井巷后風網系統的變化,輔助進行短期和長期通風系統規劃,輔助礦井通風系統進行風機選型,輔助進行礦井通風經濟性分析,模擬煙霧、粉塵、有害氣體擴散路徑和濃度,輔助災害和緊急情況處理預案制定,可對煤礦井下熱源、濕源和冷源進行分析,支持深井空氣壓縮分析、主要風路經濟斷面選型、風網通風經濟性和通風能力分析,通過在礦井中的不同情況,基于動態時間分析污染物、氣體、柴油顆粒和熱擴散、串聯通風和循環風預測,可對井下柴油機排放顆粒物濃度進行模擬[10-11].
1) 建立模型。通過在采掘工程平面圖上的導線點標高,提取巷道中心線,導入到Ventsim中,通過轉化成巷道功能將中心線轉化成三維巷道模型。模型見圖1.

圖1 15采區三維通風系統模型圖
2) 錄入數據。雙擊三維巷道,在提供的屬性對話框中錄入巷道的名稱、形狀、斷面尺寸等信息,特別是要錄入阻力測定的巷道摩擦阻力系數,另外長度如果和Ventsim依據坐標計算的有較大差別,也應該錄入。對于有調節設施的巷道,應盡量測得兩段壓差,如若不能,則可暫時設置為固定風量分支,錄入漏風量即可。斜溝風機將代表風硐的巷道設置為風機分支,然后錄入風機特性曲線。與其他采區連通的巷道,在負壓變化不大的情況下,暫時設置為固定風量分支。15采區阻力測定部分數據見表1.

表1 15采區阻力測定部分結果匯總表
3) 模擬解算。數據錄入完成后,進行試解算,然后將解算結果與實測風量進行對比,模擬數值與實測值的對比結果見表2.誤差在5%以內,此模型可用于進一步優化方案仿真。

表2 實際風量與模擬風量對比表
經過對15采區通風系統的分析,提出兩種優化方案,并分別進行了模擬,均達到了理想的優化效果。
1) 優化方案一。
主要考慮供給6#煤層系統巷的風可以同時作為16502工作面進風,增加16502工作面進風的同時,減少系統總的通風阻力。
目前,供給6#煤層的風一部分供給16502工作面,一部分供給帶式輸送機和輔運下山,實際經過系統巷后的風可以再通往工作面。具體措施是從16502工作面進風巷道打一條短巷到集中回風下山延長段,工作面原來進風變成回風,新鮮風由6#煤帶式輸送機下山和輔運下山經材料巷到工作面。同時,去掉與材料巷連接的系統巷風窗,16502通風系統優化方案見圖2.
綜合分析發現系統巷耗費了大量新鮮風,應適當加固風門,調節風窗減少系統巷供風。方案一對15采區回風上山調節措施增大阻力;集中帶式輸送機下山延長段18605工作面進風后,風門進行加固。
優化前后15采區工作面風量變化對比見表3,另外工作面優化前斜溝風機風量8 473 m3/min,風壓1 623.5 Pa;優化后風量8 635 m3/min,風壓1 453.8 Pa.優化調節效果明顯。

表3 方案一優化前后風量對比表

A—新掘巷道 B—加風墻(風門)封堵 C—拆除風窗 D—加風墻(風門)封堵圖2 優化方案一示意圖
2) 優化調節方案二。
由表2可知,15采區集中回風下山延長段風量為7 440 m3/min,風速為5.2 m/s,風量過于集中,增加了礦井通風阻力,增大了通風成本。方案二設計增加一條專用回風道,見圖3. 圖中A為新增加的15采區專用回風巷,設計斷面為20 m2.
利用VentSim對方案進行模擬,其中新掘巷道摩擦阻力系數依據經驗值,按照錨噴支護,普通爆破斜巷取0.010 7 Ns2/m8,優化前后用風地點風量對比見表4.

表4 方案二優化前后風量對比表

A—新掘巷道 B—加風墻(風門)封堵 C—加風墻(風門)封堵圖3 優化方案二示意圖
優化前斜溝風機風量8 473 m3/min,風壓1 623.5 Pa;優化后風量8 622 m3/min,風壓1 467 Pa.
利用Ventsim對斜溝煤礦15采區通風系統進行了優化模擬仿真,提出了兩種通風系統優化方案,并進行了仿真實驗。由仿真結果可知,兩種方案在優化后系統總阻力均有所下降,分別是170 Pa、156 Pa,風量均有所增加,總風量均分別增加了162 m3/min、149 m3/min.兩種方案優化效果相當,在實際改造中可根據經濟性進行選擇。