研究煤礦瓦斯和煤塵的監(jiān)測與控制11

1、關(guān)于煤礦瓦斯和煤塵監(jiān)測與控制的研究【摘要】 煤礦安全生產(chǎn)是我國目前亟待解決的問題之一，做好井下瓦斯和煤塵的監(jiān)測與控制是實現(xiàn)安全生產(chǎn)的關(guān)鍵環(huán)節(jié)。本題研究空氣中瓦斯的濃度，煤塵懸浮于空氣中的濃度，對煤礦安全的威脅。因此我們要找出引起瓦斯，煤塵濃度的因素。參照監(jiān)測表的數(shù)Excle,lingo,Matlab軟件進行求解。 在問題一中，需要對“低瓦斯礦井”還是“高瓦斯礦井”進行評估，根據(jù)煤礦安全規(guī)程第一百三十三條提供的分類標準。高瓦斯煤礦的標準瓦斯相對涌出量大于10m3/t或絕對涌出量大于40m3/min; 低瓦斯礦井相對涌出量小于10m3/t 且絕對涌出量小于或等于40m3/min，根據(jù)附表2的數(shù)據(jù)算

2、出總回風量的相對和絕對瓦斯涌出量具體公式（見模型一），然后算出三個工作點瓦斯絕對量之和，基本等于總回風巷絕對瓦斯涌出量?？偦仫L巷流出的相對瓦斯量應該與三個工作點瓦斯流出量基本相同。求得瓦斯相對涌出量為23.1832m3/t，絕對瓦斯涌出量為9.7535m3/min，因而該礦井為高瓦斯礦井。 在問題二中，需要對煤礦不安全的程度（即發(fā)生爆炸事故的可能性）進行概率計算，根據(jù)煤礦安全規(guī)程第一百六十八條的規(guī)定，并參照附表1，對可能發(fā)生爆炸的點進行統(tǒng)計為25個，而發(fā)生爆炸的量為1620個，即瓦斯、煤塵和瓦斯與煤塵混合物的濃度在三個工作面與三個回風巷在一個月三個班次可能發(fā)生爆炸事故的量，最后將各點可能發(fā)生爆

3、炸的概率求和，得出該煤礦不安全的程度（即發(fā)生爆炸事故的可能性）為1.5431%，所以該礦井發(fā)生爆炸的可能性較大。因此我們要采取相應措施減少爆炸機率。 在問題三中，需對主要風巷的風速、瓦斯?jié)舛?、煤塵濃度之間關(guān)系已最佳風通為目標建立優(yōu)化模型運用Matlab軟件擬出風速與瓦斯?jié)舛龋L速與煤塵濃度的圖形。進而確立風速分別與瓦斯?jié)舛?、煤塵的函數(shù)關(guān)系根據(jù)圖像進行擬合，并運用Lingo軟件，最小二乘法求出各函數(shù)的誤差。然后我們把最佳風速定義為：在滿足安全生產(chǎn)，高效生產(chǎn)的情況下的基礎上定為最佳風速。通過計算得到最佳風速：工作面一3.29m/s 工作面二2.37m/s，掘進工作面0.85m/s，回風巷一1.34

4、m/s，回風巷二2.4m/s，總回風巷6m/s，最佳風量1800m3/min、局部通風量為203.17m3/min。關(guān)鍵詞：Excel Matlab擬合 Lingo 相對瓦斯涌出量 絕對瓦斯涌出量 最小二乘法1、問題重述煤礦安全生產(chǎn)是我國目前亟待解決的問題之一，做好井下瓦斯和煤塵的監(jiān)測與控制是實現(xiàn)安全生產(chǎn)的關(guān)鍵環(huán)節(jié)（見附件1）。瓦斯是一種無毒、無色、無味的可燃氣體，其主要成分是甲烷，在礦井中它通常從煤巖裂縫中涌出。瓦斯爆炸需要三個條件：空氣中瓦斯達到一定的濃度；足夠的氧氣；一定溫度的引火源。煤塵是在煤炭開采過程中產(chǎn)生的可燃性粉塵。煤塵爆炸必須具備三個條件:煤塵本身具有爆炸性；煤塵懸浮于空氣中并

5、達到一定的濃度；存在引爆的高溫熱源。試驗表明，一般情況下煤塵的爆炸濃度是30 2000g/m3，而當?shù)V井空氣中瓦斯?jié)舛仍黾訒r，會使煤塵爆炸下限降低，結(jié)果如附表1所示。國家煤礦安全規(guī)程給出了煤礦預防瓦斯爆炸的措施和操作規(guī)程，以及相應的專業(yè)標準 (見附件2)。規(guī)程要求煤礦必須安裝完善的通風系統(tǒng)和瓦斯自動監(jiān)控系統(tǒng)，所有的采煤工作面、掘進面和回風巷都要安裝甲烷傳感器，每個傳感器都與地面控制中心相連，當井下瓦斯?jié)舛瘸瑯藭r，控制中心將自動切斷電源，停止采煤作業(yè)，人員撤離采煤現(xiàn)場。具體內(nèi)容見附件2的第二章和第三章。附圖1是有兩個采煤工作面和一個掘進工作面的礦井通風系統(tǒng)示意圖，請你結(jié)合附表2的監(jiān)測數(shù)據(jù)，按照煤

6、礦開采的實際情況研究下列問題： （1）根據(jù)煤礦安全規(guī)程第一百三十三條的分類標準 (見附件2)，鑒別該礦是屬于“低瓦斯礦井”還是“高瓦斯礦井”。（2）根據(jù)煤礦安全規(guī)程第一百六十八條的規(guī)定，并參照附表1，判斷該煤礦不安全的程度（即發(fā)生爆炸事故的可能性）有多大？ （3）為了保障安全生產(chǎn)，利用兩個可控風門調(diào)節(jié)各采煤工作面的風量，通過一個局部通風機和風筒實現(xiàn)掘進巷的通風（見下面的注）。根據(jù)附圖1所示各井巷風量的分流情況、對各井巷中風速的要求（見煤礦安全規(guī)程第一百零一條），以及瓦斯和煤塵等因素的影響，確定該煤礦所需要的最佳（總）通風量，以及兩個采煤工作面所需要的風量和局部通風機的額定風量（實際中，井巷可能

7、會出現(xiàn)漏風現(xiàn)象）。注：掘進巷需要安裝局部通風機，其額定風量一般為150400 m3/min。局部通風機所在的巷道中至少需要有15%的余裕風量（新鮮風）才能保證風在巷道中的正常流動，否則可能會出現(xiàn)負壓導致乏風逆流，即局部通風機將乏風吸入并送至掘進工作面。名詞解釋：（1）采煤工作面：礦井中進行開采的煤壁 (采煤現(xiàn)場)。（2）掘進巷：用爆破或機械等方法開鑿出的地下巷道，用以準備新的采煤區(qū)和采煤工作面。（3）掘進工作面：掘進巷盡頭的開掘現(xiàn)場。（4）新鮮風：不含瓦斯和煤塵等有害物質(zhì)的風流。（5）乏風：含有一定濃度的瓦斯和煤塵等有害物質(zhì)的風流。2、問題分析煤礦安全生產(chǎn)是目前社會重點關(guān)注的熱點問題之一，尤其

8、是在能源緊張，對煤碳的需求量不斷增加的情況下，煤礦的安全生產(chǎn)問題更是值得我們關(guān)注，這也是建設平安和諧社會的重要組成部分。大部分事故的罪魁禍首都是瓦斯或煤塵爆炸，瓦斯在煤礦的開采中是不可避免的。因此，礦井下的瓦斯和煤塵對煤礦的安全生產(chǎn)構(gòu)成了重大威脅，做好煤礦井下瓦斯和煤塵的監(jiān)測與控制是保證煤礦安全生產(chǎn)的關(guān)鍵所在。2.1對問題一分析對題目進行反復研究，我們從題目中不難得出，如何判斷此煤礦的瓦斯類別，主要是根據(jù)礦井相對瓦斯涌出量和礦井絕對瓦斯涌出量的范圍來進行判斷，解決問題一主要是從以下幾個方面來考慮:第一個方面根據(jù)煤礦通風系統(tǒng)示意圖，確定瓦斯產(chǎn)出的部位，反復讀題和研究煤礦通風系統(tǒng)示意圖，可以得出結(jié)

9、論，瓦斯主要是在采煤過程當中產(chǎn)生，最后瓦斯的涌出主要是通過總回風巷，跟回風巷和回風巷無關(guān)，所以我們確定總回風巷的瓦斯涌出量為礦井瓦斯涌出量；第二個方面是如何對礦井相對瓦斯涌出量和礦井絕對瓦斯涌出量進行定義，研究題目，根據(jù)題目所給的礦井相對瓦斯涌出量的單位，可以確定礦井相對瓦斯涌出量的定義為產(chǎn)一個單位煤時涌出多少瓦斯，而礦井絕對涌出量是瓦斯涌出與工作時間的關(guān)系，可以表示為在一個工作時間內(nèi)涌出多少瓦斯。根據(jù)以上分析的結(jié)果，我們可以建立一個計算標準，把得到的結(jié)果與標準進行對比來判斷該礦井的瓦斯類別。2.2對問題二分析如何確定該礦發(fā)生安全事故的可能性，明顯需要求解的是概率問題，根據(jù)附件2煤礦安全規(guī)程第

10、一百六十八條，煤礦產(chǎn)生爆炸的原因主要有三種情況，一個是瓦斯?jié)舛瘸瑯?，第二個是煤塵濃度超標，最后是他們的中和濃度屬于可爆炸范圍，這些情況都會引起煤礦爆炸的可能性發(fā)生，認真分析本題后兩種情況在該礦中都不會發(fā)生，對于問題二我們只需分析甲烷傳感器的報警濃度范圍，當甲烷濃度大于等于百分之一時，甲烷傳感器就會報警，礦井下面就有可能發(fā)生安全事故；分析附表2當中甲烷在每個班次每個檢測點的瓦斯?jié)舛?，把等于或者超過百分之一的選擇中出來，確定瓦斯不達標準的頻率，再與總的檢測次數(shù)進行對比，我們可以確定煤礦發(fā)生瓦斯爆炸的概率。2.3對問題三分析問題三是根據(jù)附圖1研究各井巷風量的分流情況、對各井巷中風速的要求，以及瓦斯和

11、煤塵等因素的影響，確定最佳通風量，以及兩個采煤工作面所需要的風量和局部通風機的額定風量。實際中，井巷可能會出現(xiàn)漏風現(xiàn)象，由于漏風是很難估計的，與很多因素有關(guān)。而且相對于總風量，漏風只是少數(shù)可以忽略不計。所以假設井巷不漏風。在煤礦中，降低瓦斯?jié)舛鹊奈ㄒ晦k法是通過通風系統(tǒng)將瓦斯排除到井外，從降低瓦斯?jié)舛鹊慕嵌瓤紤]的話，巷道中風速越大越好。但是風速越大，煤層又會增加，而煤層同樣會引起爆炸。所以要通過控制一定量的通風量來使得瓦斯?jié)舛群兔簤m達到一定的平衡而使其不發(fā)生爆炸。在達到這個條件下總風量越小越好，所以以總風量最小為目標函數(shù)，再根據(jù)題目當中的條件確定目標函數(shù)的約束條件，再用lingo編程求解，從而解

12、決第三個問題。2、 模型假設1、 假設在開采過程中只產(chǎn)生瓦斯和煤塵。2、 假設該礦井檢測的數(shù)據(jù)是絕對可信準確。3、 假設該礦井在開采過程中嚴格遵守安全規(guī)程。4、 假設礦井通風示意圖嚴格反應出該礦井內(nèi)部結(jié)構(gòu)。5、 假設礦井通風示意圖嚴格反應出該礦井內(nèi)部結(jié)構(gòu)。6、 假設該礦井在檢測過程中監(jiān)控和傳感器裝備不存在故障。 4、符號說明符號描述說明井礦絕對瓦斯涌出量井礦相對瓦斯涌出量第天第個班次的風速第天第個班次的瓦斯第的日產(chǎn)量為時間常數(shù)（480min）第個工作點的風速第個工作點瓦斯與風速的關(guān)聯(lián)系數(shù)第個工作點煤塵與風速的關(guān)聯(lián)系數(shù)5、模型建立與求解本題共三問，前兩問由問題分析可知，主要是通過理解單位定義，利

13、用數(shù)學理解來建立關(guān)系式來進行求解。第三問是一個相互牽制與約束關(guān)系，一個方面要求出瓦斯?jié)舛扰c風速之間的函數(shù)關(guān)系，根據(jù)函數(shù)關(guān)系分析誤差范圍，然后求出總風量最小，為求得目標可以通過建立一個優(yōu)化模型來求解。5.1判斷煤礦類別根據(jù)題意，如何判斷煤礦的類別，主要是根據(jù)原題當中的附件2煤礦安全規(guī)程第一百三十三條來判斷，由問題一的分析可知，我們主要是明確以下幾個方面：第一個分析礦井結(jié)構(gòu)，第二個如何定義相對涌出量，絕對涌出量的表達意義，第三個是確定關(guān)系式如何計算相對涌出量，絕對涌出量，最后是計算出結(jié)果與判斷標準進行比較。解決問題一我們主要通過以下四步來進行。圖一：判斷煤礦類別流程圖1）如何確定礦井瓦斯涌出量對題

14、目進行分析，礦井瓦斯產(chǎn)生主要是從煤層裂縫中涌出，可以確定產(chǎn)生瓦斯的部位主要是挖煤部位，他們分別是工作面，工作面，掘進工作面。他們涌出的瓦斯通過回風巷，回風巷最后從總回風巷涌出，由此而分析得知，總回風巷涌出的瓦斯具有一定的代表說明意義。根據(jù)如何計算出流體流量的定義，我們確定礦井瓦斯涌出量可以表示為：瓦斯的濃度與風速和總回風巷截面積之間的關(guān)系，具體表示式如下：2）根據(jù)單位理清相對涌出量與絕對涌出量意義分析題目，相對涌出量的單位表示為立方米每噸，說明的是產(chǎn)出的煤與瓦斯之間的關(guān)系，可具體表示為產(chǎn)一噸煤的同時涌出多少瓦斯量；絕對涌出量的單位表示為立方米每分鐘，說明的是產(chǎn)出的煤與時間之間的關(guān)系，可具體表示

15、為在一個工作班次內(nèi)，工作一分鐘涌出多少瓦斯量。3）建立表達式與計算根據(jù)以上兩步分析，我們確定了相對涌出量與絕對涌出量的計算關(guān)系式，建立了如下表達式：4）計算結(jié)果與判斷 我們把附表2相關(guān)數(shù)據(jù)整理成excel表格，把以上關(guān)系式轉(zhuǎn)化為公式運用到excel表格當中，當計算絕對涌出量時我們把總回風巷計算出來的結(jié)果與工作面，工作面，掘進工作面計算出來的結(jié)果之和進行對比。得出的數(shù)據(jù)分別是9.753立方米每分鐘，9.4305立方米每分鐘，考慮到井下漏風現(xiàn)象，結(jié)果相差不大，不影響判斷結(jié)果，可以作為標準判斷的依據(jù)；計算相對涌出量時得出的結(jié)果為23.183立方米每噸，整理結(jié)果制出下表： 表1：判斷依據(jù)表總回風巷三工

16、作面瓦斯體積 （m3）421354.656 407398.003 月產(chǎn)量(t)18175/絕對瓦斯涌出(m3/min)9.754 9.431 相對瓦斯涌出(m3/t)23.183 /總時間（min）43200 43200該礦井瓦斯相對涌出量為23.183立方米每噸遠遠大于判斷標準的10立方米每噸，可以直接判斷該礦井為高瓦斯礦井。5.2判斷煤礦爆炸事故可能性如何確定煤礦發(fā)生爆炸事故的可能性，主要通過確定可能發(fā)生爆炸事故的概率來決定，發(fā)生爆炸事故的概率，主要是通過確定在總的檢測次數(shù)當中瓦斯不達標次數(shù)出現(xiàn)的比例，具體方法如下圖： 圖2：確定煤礦爆炸概率圖1） 確定瓦斯?jié)舛瘸瑯丝赡馨l(fā)生爆炸事故的班次與

17、部位表2：瓦斯超標的班次與對應工作點分析附表2的檢測數(shù)據(jù)，我們把瓦斯?jié)舛鹊扔诤统^百分之一的選擇出來，把一天的早、中、晚三個班次進行編號，30天的班次可以分為190個班次，再選出某個班次發(fā)生爆炸的相應部位，制成如下表格。班次工作面2 瓦斯?jié)舛劝啻位仫L巷2 瓦斯?jié)舛劝啻位仫L巷2 瓦斯?jié)舛劝啻位仫L巷2 瓦斯?jié)舛劝啻位仫L巷2 瓦斯?jié)舛?8191.01181.07491.05691.04251.01101241501.07751.03411.11151.03251.05511801.03501.03161.02371621851.06851171.01411.18681.02881.03 統(tǒng)計表格得到

18、數(shù)據(jù)可能會發(fā)生爆炸的次數(shù)總共有25次，對應的部位主要是工作面和回風巷。2）統(tǒng)計總檢測次數(shù)與計算概率根據(jù)附表2的統(tǒng)計情況，一天三個班次分別統(tǒng)計六個部位可以得到總的檢測次數(shù)為：天數(shù)*班次*檢測部位=540（次），產(chǎn)生爆炸的次數(shù)等于主要考慮瓦斯超標可能引發(fā)事故的次數(shù)為：25（次），我們建立如下關(guān)系式進行計算。 把統(tǒng)計的次數(shù)運用到以上關(guān)系式當中進行計算，得出P=4.629%。3） 再分析結(jié)果因為爆炸有三種情況他們分別是：一個是瓦斯?jié)舛瘸瑯?，第二個是煤塵濃度超標，最后是他們的中和濃度屬于可爆炸范圍，這些情況都會引起煤礦爆炸的可能性發(fā)生。我們以上所謂統(tǒng)計的總檢測次數(shù)，實際上是瓦斯的總檢測次數(shù)。單獨依據(jù)以上

19、算法是沒有任何依據(jù)，算出來的概率會偏大，不符合實際情況，所以總檢測次數(shù)=瓦斯的檢測次數(shù)*3=1620（次），我們再次修正，進行計算得到P=1.543%，最后確定該礦的不安全程度為1.543%。5.3確定所需最佳通風量 5.3.1擬合出瓦斯?jié)舛扰c風速的關(guān)系1）如何確定總風量最小，首先要通過觀察風速與瓦斯?jié)舛戎g的關(guān)系，整理附表2的相關(guān)風速和瓦斯?jié)舛鹊臄?shù)據(jù)，再運用matlab畫圖命令，執(zhí)行操作，可以得到瓦斯?jié)舛扰c風速之間的散點圖，觀察圖形，可以假設風速與瓦斯?jié)舛戎g的關(guān)系成一個函數(shù)關(guān)系，這個方法就叫擬合，擬合圖形如下：圖3：瓦斯?jié)舛扰c風速擬合圖 觀察瓦斯?jié)舛扰c風速之間的關(guān)系，我們可以看出瓦斯隨著風

20、速的增加濃度減小，可以假設一個擬合函數(shù)能反映他們之間的關(guān)系，擬合函數(shù)為：運用最小二乘法，算出最小誤差，建立如下求最小誤差模型： 把函數(shù)編寫成lingo程序見（附錄1），可以得到、，整理數(shù)據(jù)得到下表：表3：瓦斯?jié)舛扰c風速擬合表工作面1工作面2掘進面-0.26131.28750.0309-0.45681.83280.0833-0.80410.40990.4236回風巷1回風巷2總回風巷-0.23141.21110.033-0.63112.26090.085-0.11941.24620.03345.3.2擬合出煤塵濃度與風速的關(guān)系需要明確煤塵濃度與風速之間的關(guān)系，我們接著把題目當中表2的煤塵濃度與風

21、速的數(shù)據(jù)運用matlab擬合出圖形，擬合圖形如下： 圖4：煤塵濃度與風速擬合圖觀察擬合圖，可以得到煤塵濃度與風速之間的關(guān)系，我們可以看出煤塵隨著風速的增加濃度出現(xiàn)不穩(wěn)定的增大，可以假設一個擬合函數(shù)能反映他們之間的關(guān)系，擬合函數(shù)為：同樣把函數(shù)編寫成lingo程序見（附錄1），可以得到煤塵函數(shù)的、，整理數(shù)據(jù)得到下表：表4：煤塵濃度與風速擬合表工作面1工作面2掘進面abcabcabc1.98073.21660.0941.24855.01250.09141.94322.99320.1032回風巷1回風巷2總回風巷abcabcabc1.25464.83590.17650.81655.50910.1195

22、0.22975.86480.12045.3.3建立優(yōu)化模型求最佳通風量和兩個采煤工作面所需要的風量和局部通風機的額定風量仔細分析題目，要求得最佳通風量，這是一個多條件約束問題，主要有以下幾個方面：各井巷風量的分流情況、對各井巷中風速的要求，以及瓦斯和煤塵等因素的影響，確定該煤礦所需要的最佳（總）通風量，以及兩個采煤工作面所需要的風量和局部通風機的額定風量。（1） 確定總風量最小為目標要求得最佳通風量實際上，是需要求出在約束條件下，最小通風量，所以我們建立如下目標函數(shù)： （二）建立約束條件（1）根據(jù)附件表2，可知井巷中允許的風流速度，采煤工作面的允許風速為最低不得低于0.25m/s,最高不得高于

23、4m/s?；仫L巷的允許風速為最低不得低于0.25m/s,最高不得高于6m/s，建立如下約束條件：（2） 根據(jù)附件表2，可知礦井中允許的風流速度，采煤工作面和回風巷的風流速度不能超過8m/s，建立如下約束條件：（3） 工作面、工作面、掘進工作面，主巷道的風量之和等于總回風巷的風量，巷道中至少需要有15%的余裕風量（新鮮風）才能保證風在巷道中的正常流動，建立如下約束條件：（4） 如何保證該礦井是在安全的條件下工作，需要保證每個班次每個檢測點的瓦斯?jié)舛炔荒艹^1，建立如下約束條件：（5）根據(jù)表2：在保證每個班次每個檢測點的瓦斯?jié)舛炔荒艹^1%的前提下，煤塵的濃度不能超過15g/m3，建立如下約束條件

24、：（5） 根據(jù)掘進巷需要安裝局部通風機，其額定風量一般為150400 m3/min，否則可能會出現(xiàn)負壓導致乏風逆流，即局部通風機將乏風吸入并送至掘進工作面，建立如下約束條件：（6） 根據(jù)附表2第一百零一條相關(guān)規(guī)定，掘進工作面的風速需要小于回風巷的風速，工作面的風速小于回風巷的風速，根據(jù)這些規(guī)定，建立如下約束條件：（3） 模型確定與求解根據(jù)以上的目標函數(shù)和約束條件，我們建立模型，從模型中可以求解出最佳通風量，以及兩個采煤工作面所需要的風量和局部通風機的額定風量，具體模型如下： 把模型編寫成lingo程序（見附錄）運行程序得出結(jié)果（見附錄）得到每個檢測點的最佳風速，然后結(jié)合公式進行計算，得到最佳通

25、風量為1800.00，工作面所需風量為788.96，工作面所需風量為568.84，局部通風機額定風量為203.17，整理數(shù)據(jù)制出下表：表5：問題三最后結(jié)果工作點X( 1)X( 2)X( 3)X( 4)X( 5)X( 6)風速（m/s）3.29 2.37 0.85 1.34 2.40 6.00 最佳通風量（m3/min）1800.00 采煤工作面1（m3/min）788.96 采煤工作面2（m3/min）568.84 局部通風機（m3/min）203.17 6、模型的評價與推廣本文根據(jù)題中所給條件，建立了整數(shù)規(guī)劃優(yōu)化，由于題中限制條件過多，運用Lingo軟件快速求解，減少了運算量、Matlab軟

26、件擬合整理數(shù)據(jù)，形成相應的函數(shù)對應關(guān)系式，方便快捷。6.1模型優(yōu)點1）此模型應用面較廣，使用類似的的方法，集中生產(chǎn)等線性規(guī)劃問題都可適用。2）用Lingo軟件求解，方便快捷；Matlab軟件擬合，結(jié)果簡潔明了，參考價值高，利于問題分析求解。3）運用整數(shù)線性規(guī)劃，優(yōu)化模型，Matlab擬合等多種方法建立模型，解決了風速分別與瓦斯和煤塵濃度之間的函數(shù)關(guān)系，并求解出誤差范圍等。4）本模型建立煤礦的風速大小、瓦斯?jié)舛扰c煤塵濃度之間最佳工作數(shù)值范圍的問題研究，能夠推廣至其他方面。5）該模型整個過程思路清晰，結(jié)構(gòu)明了，可行性強，易于推廣。6.2模型缺點 1）Lingo、Matlab軟件運算求解時，計算機在

27、擬合和數(shù)據(jù)（y=a*x+b±3*）帶入求解存在一定的誤差。2）由于現(xiàn)實情況復雜，沒有考慮礦井中各個工作檢測點的環(huán)境存在其他特殊情況等。6.3模型推廣該問題運用優(yōu)化模型規(guī)劃方法，可以解決經(jīng)濟管理中大量的規(guī)劃問題。如工廠產(chǎn)品分配問題、部門工作人員的工作與休息時間的安排問題、工廠選址問題等?；镜姆椒ㄟ€是在滿足特定的需求下，求目標函數(shù)的最優(yōu)化值。 附錄1：sets:aa/1.6/:x,a,b,e,g,h,d;endsetsdata:a=-0.2613 -0.4568-0.8041 -0.2314 -0.6311-0.1194;b=1.2875 1.8328 0.4099 1.21112.2

28、609 1.2462;h=0.0309 0.0833 0.4236 0.0330.085 0.0334;e=1.9807 1.2485 1.9432 1.25460.8165 0.2297;g=3.2166 5.0125 2.9932 4.8359 5.5091 5.8648;d=0.094 0.0914 0.1032 0.1765 0.1195 0.1204;enddatamin=5*x(6)*60;for(aa(i)|i#le#3:x(i)>=0.25);for(aa(i)|i#le#3:x(i)<=4);for(aa(i)|i#ge#4#le#5:x(i)>=0.25)

29、;for(aa(i)|i#ge#4:x(i)<=6);for(aa(i):x(i)<=8);4*x(3)/0.85+4*(x(1)+x(2)+x(3)=5*x(6);for(aa(i):a(i)*x(i)+b(i)+3*h(i)<1);for(aa(i):e(i)*x(i)+g(i)+3*d(i)<15);240*x(3)>150;附錄2：Global optimal solution found. Objective value: 1360.243 Infeasibilities: 0.000000 Total solver iterations: 0 Variable Value Reduced Cost X( 1) 1.455033 0.000000 X( 2) 2.370184 0.000000 X( 3) 0.84653