楊錚園， 申立勇

(中國(guó)科學(xué)院大學(xué)數(shù)學(xué)科學(xué)學(xué)院，北京101408)

數(shù)據(jù)點(diǎn)擬合在計(jì)算機(jī)輔助幾何設(shè)計(jì)、計(jì)算機(jī)圖形學(xué)和數(shù)控加工等領(lǐng)域都有廣泛的應(yīng)用，擬合中較常使用B樣條曲線。B樣條曲線具有較好的造型能力[1-2]。在擬合算法及提高精確度和效率方面，前人已做了許多研究。通常把平方距離和作為擬合誤差，這樣擬合過程就可以描述成一個(gè)求最小擬合誤差的最優(yōu)化問題。在優(yōu)化迭代中，怎樣估計(jì)擬合誤差是一個(gè)關(guān)鍵的問題。目前已有很多估計(jì)擬合誤差的方法，包括PDM，TDM，SDM，GTDM和CDM[3-4]等，在擬合效果和速度上都各有優(yōu)缺點(diǎn)。Yang等[5]提出了隱式B樣條曲線擬合算法，適用于更一般的情形，比如幾條分離的閉曲線的擬合。Fl?ry[6]提出了在障礙存在下的B樣條曲線曲面擬合算法，該算法能夠進(jìn)行點(diǎn)云邊界的重構(gòu)，并推廣到避開一般真實(shí)障礙物的曲線曲面重構(gòu)。本文針對(duì)實(shí)際實(shí)驗(yàn)數(shù)據(jù)，采用多約束條件下的樣條擬合，處理油氣爆炸界限的估計(jì)問題，得到不同氧氣濃度下油氣濃度和爆炸強(qiáng)度曲線，應(yīng)用于燃爆預(yù)判。

近年來，工業(yè)生產(chǎn)過程危險(xiǎn)化學(xué)品事故頻繁發(fā)生，事故造成了生命和財(cái)產(chǎn)的重大損失。因此，危險(xiǎn)化學(xué)品事故的預(yù)測(cè)是眾多學(xué)者積極研究的主題[7-9]。Chang和Lin[10]整理了在2006年之前的40年間有記載的工廠事故，其中約74%的事故發(fā)生在石油精煉廠、石油港口和油庫(kù)。近年來這類事故的數(shù)量還在上升，包括2010年大連港油庫(kù)爆炸事故。油罐爆炸可能會(huì)導(dǎo)致多米諾骨牌效應(yīng)，加劇意外事故所造成的損害[11]。油罐爆炸主要是油氣燃爆引起的，因此對(duì)油氣爆炸范圍進(jìn)行分析非常必要，良好的分析結(jié)果將有助于防止油氣爆炸事故的發(fā)生。之前人們主要應(yīng)用不同分類器對(duì)某種油氣的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行分析，通過訓(xùn)練得到合適的分類器參數(shù)，從而可以對(duì)新數(shù)據(jù)進(jìn)行是否屬于可以燃爆范圍的預(yù)判。常用的分離器包括支持向量機(jī)(SVM)方法[12]、邏輯回歸(LR)方法[13]等。這類基于統(tǒng)計(jì)分析的分類器一般適合給出邏輯變量分類，即是否可燃爆。例如目前常用的燃?xì)馊莘e比作為因變量，而它的可燃性作為目標(biāo)變量。但是實(shí)際生產(chǎn)中人們不僅關(guān)注燃?xì)馐欠癖?，還關(guān)注其爆炸程度強(qiáng)弱，以便對(duì)后援工作進(jìn)行有效安排。因此，本文基于此實(shí)際需求，對(duì)某油罐燃?xì)饪扇夹砸约叭急瑥?qiáng)度進(jìn)行分析，其影響因素涉及其中的油氣濃度和氧氣濃度等。我們通過在特制的密閉管道內(nèi)進(jìn)行油氣爆炸試驗(yàn)，從而獲取在常溫常壓條件下不同油氣濃度、氧氣濃度預(yù)混氣體的爆炸極限及爆炸強(qiáng)度值(圖1)。然后分析不同氧氣濃度下的油氣燃爆界，以及可燃爆情況下的爆炸強(qiáng)度。

圖1 實(shí)驗(yàn)裝置示意圖(單位：mm)

本文在關(guān)注油氣數(shù)據(jù)擬合中，首先根據(jù)問題的物理化學(xué)背景和實(shí)驗(yàn)數(shù)據(jù)特征，對(duì)擬合曲線的預(yù)期性質(zhì)進(jìn)行估計(jì)。然后把這些性質(zhì)轉(zhuǎn)化為曲線擬合中的約束條件，對(duì)擬合曲線進(jìn)行帶約束的優(yōu)化。包括，不同氧氣濃度下的爆炸強(qiáng)度曲線應(yīng)該互不相交；爆炸強(qiáng)度曲線是凸曲線；爆炸強(qiáng)度曲線應(yīng)該為安全曲線等。這些要求在轉(zhuǎn)化為約束條件時(shí)又有不同的轉(zhuǎn)化和實(shí)現(xiàn)方式，需要不同的嘗試。本文首先回顧了B樣條曲線擬合的一般過程、SDM誤差項(xiàng)估計(jì)方法和Fl?ry的障礙約束方法，然后根據(jù)數(shù)據(jù)特征和問題背景導(dǎo)出曲線組擬合的約束條件，最后在實(shí)際實(shí)驗(yàn)數(shù)據(jù)上實(shí)現(xiàn)算法，在分析、比較擬合效果后選擇合適的方案。

1 實(shí)驗(yàn)數(shù)據(jù)和樣條擬合

1.1 實(shí)驗(yàn)數(shù)據(jù)采集

在實(shí)驗(yàn)裝置中(圖1)，首先測(cè)量實(shí)驗(yàn)油氣(HC)濃度、氧氣(O2)濃度和惰性氣體(CO2)濃度，惰性氣體主要是用來調(diào)整前兩種實(shí)驗(yàn)氣體濃度。然后通過點(diǎn)火裝置，再由傳感器測(cè)量并記錄引爆后的各自氣體濃度，包括油氣、氧氣、一氧化碳、二氧化碳。在3個(gè)位置安裝了濃度傳感器，將數(shù)據(jù)平均值作為最終實(shí)驗(yàn)數(shù)據(jù)進(jìn)行分析。另外，還安裝了5個(gè)壓強(qiáng)傳感器，測(cè)量爆炸瞬間的壓強(qiáng)，其平均值作為實(shí)驗(yàn)壓強(qiáng)值，該瞬時(shí)壓強(qiáng)稱之為燃?xì)獗◤?qiáng)度。而燃?xì)馊急缡侵腹潭ㄑ鯕鉂舛认掠蜌獍l(fā)生爆炸的濃度邊界，通常用上下界區(qū)間表示。表1和表2分別給出了引爆后發(fā)生爆炸和沒有爆炸的兩組實(shí)驗(yàn)數(shù)據(jù)；表3為爆炸瞬間壓強(qiáng)表。

整個(gè)實(shí)驗(yàn)過程中，共采集、整理了44組數(shù)據(jù)(圖2)。經(jīng)初步分析發(fā)現(xiàn)，隨著氧氣濃度的增大，油氣爆炸范圍增大、爆炸強(qiáng)度也增強(qiáng)。因此，將數(shù)據(jù)按氧氣濃度分類，然后對(duì)每一組數(shù)據(jù)進(jìn)行擬合。但是每組數(shù)據(jù)擬合不是獨(dú)立的，不同氧氣濃度下爆炸曲線有相互約束關(guān)系。

表1 爆炸實(shí)驗(yàn)數(shù)據(jù)(No.21)

表2 未爆炸實(shí)驗(yàn)數(shù)據(jù)(No.19)

表3 爆炸瞬間壓強(qiáng)(MPa)

圖2 實(shí)驗(yàn)數(shù)據(jù)

1.2 B樣條曲線擬合

設(shè)f:?→?2:t?f(t)是平面上的一條參數(shù)曲線。已知有n個(gè)控制點(diǎn)的B樣條曲線的一般表示為：

其中，Bk,i(t)∈?是k次B樣條基函數(shù)，是曲線的控制點(diǎn)。令：

則式(1)可以寫成矩陣表示：

給出平面上的一組無序數(shù)據(jù)點(diǎn)pk,k=1,…,m，用一條B樣條曲線來擬合這些點(diǎn)，使得這條曲線盡可能地描述這些點(diǎn)的軌跡形狀。這一過程可以描述為如下一個(gè)非線性優(yōu)化問題：已知點(diǎn)集求一組控制點(diǎn)c,i=1,…,n,∈?2，使i得目標(biāo)函數(shù)：

的值最小。其中|d(pk,f)|是pk到曲線f(t)的距離，ωfs是平滑項(xiàng)。

式(3)的最優(yōu)化問題可以用迭代過程解決。我們不妨設(shè)f(tk)是已知曲線f(t)上距離pk最近的點(diǎn)，pk到曲線f(t)的距離定義為|d(pk,f(tk))|。對(duì)控制點(diǎn)向量C做位移為δC的擾動(dòng)，則原始曲線f(t)變?yōu)椋?/p>

pk到曲線的距離未知，但是可以用合理的方法來估計(jì)。優(yōu)化問題的關(guān)鍵在于如何估計(jì)曲線的擬合誤差d2(pk,）。Wang等[3]提出的SDM方法給出的誤差項(xiàng)：

其中，Tk和Nk分別是f(tk)點(diǎn)的單位切向量和單位法向量，ρk是f(tk)點(diǎn)的曲率半徑。|dk|=||pk-f(tk)||，當(dāng)ρk與曲率中心K在曲線f(t)同側(cè)時(shí)dk＞0，異側(cè)時(shí)dk＜0。顯然dk＜ρk。式(5)可展開計(jì)算并求和，

其中：

另一方面，平滑項(xiàng)部分在實(shí)際應(yīng)用中常用最小化薄板能量，使用二階導(dǎo)數(shù)的平方[14]，因此fs也是C的二次函數(shù)：

于是式(3)的無約束最優(yōu)化問題等價(jià)于在現(xiàn)有值C的情況下，每次迭代中最小化如下函數(shù)：

根據(jù)以上記號(hào)，B樣條曲線擬合算法可描述如下：

輸入：需要擬合的一組二維無序數(shù)據(jù)點(diǎn)pk,k=1,…,m,∈?2。

輸出：一條B樣條擬合曲線f(t)。

(1) 根據(jù)所要擬合的數(shù)據(jù)點(diǎn)，合理設(shè)置初始控制點(diǎn),i=1,…,n，根據(jù)式(2)計(jì)算出初始的B樣條曲線f0(t)。

(2) 對(duì)每個(gè)數(shù)據(jù)點(diǎn)pk，計(jì)算其在當(dāng)前的曲線fr(t)上的最近點(diǎn)fr(tk)。

(3) 通過求解式(8)的最小值，計(jì)算擬合曲線控制點(diǎn)Cr的位移δCr。則新的擬合曲線為fr+1(t)=LT(t)(Cr+δCr)。如果新的擬合曲線滿足了事先設(shè)定的條件（比如迭代步數(shù)或者誤差閾值的控制），則轉(zhuǎn)到第4步；否則轉(zhuǎn)到第2步。

(4) 輸出擬合曲線fr+1(t)。

在樣條擬合的一般算法的基礎(chǔ)上，F(xiàn)l?ry[6]提出了在一般障礙物存在的情況下擬合點(diǎn)云的方法。

同上文類似的記號(hào)，并記f()處的向外單位法向量n)，則式(9)的約束等價(jià)于：

由于每個(gè)數(shù)據(jù)點(diǎn)對(duì)應(yīng)的與當(dāng)前的擬合曲線的控制點(diǎn)相關(guān)，故可寫成f(t)的控制點(diǎn)向量C的函數(shù)，并記式(10)左側(cè)的線性近似為′(C)，則：

因此，約束(10)可以近似寫成：

′(C)對(duì)每個(gè)控制點(diǎn)ci求偏導(dǎo)：

其中，I是二階單位矩陣，J是將切向量向曲線外旋轉(zhuǎn)π/2的矩陣。經(jīng)整理后可寫出梯度矩陣?Clk′(C)，式(12)的具體推導(dǎo)過程可參見文獻(xiàn)[6]。

2 多約束的B樣條曲線擬合

根據(jù)爆炸數(shù)據(jù)的物理化學(xué)特征，將不同的擬合約束轉(zhuǎn)化為合理的數(shù)學(xué)條件，并對(duì)條件做進(jìn)一步簡(jiǎn)化。

2.1 多條互不相交的曲線擬合

爆炸反應(yīng)的原理顯示，在相同的油氣濃度下，爆炸強(qiáng)度的大小跟氧氣濃度有關(guān)。不同氧氣濃度下的爆炸強(qiáng)度曲線應(yīng)該互不相交。因此在實(shí)際應(yīng)用于爆炸數(shù)據(jù)的擬合中，需要針對(duì)曲線性質(zhì)的特殊要求，對(duì)擬合曲線進(jìn)行帶約束的優(yōu)化。對(duì)于本文的實(shí)驗(yàn)數(shù)據(jù)，不能直接應(yīng)用上一節(jié)的方法。因?yàn)镕l?ry的算法處理的障礙物都是處于靜止?fàn)顟B(tài)下，擬合的多條曲線分別作為其他曲線的邊界約束在每一次迭代中會(huì)發(fā)生變化，因此約束曲線上取的樣點(diǎn)作為障礙點(diǎn)也在變化。

上述問題的擬合可以用兩種方法來完成。首先很直接地想到第一種方法：先擬合其中一條曲線，比如f2(t2)，然后擬合另一條曲線f1(t1)時(shí)，在f2(t2)上取足夠多的樣點(diǎn)作為障礙點(diǎn)集，即可用上一節(jié)的帶約束的算法來擬合f1(t1)。這樣的做法能夠滿足曲線互不相交，但是曲線擬合順序的不同可能導(dǎo)致產(chǎn)生不同的擬合效果，尤其在兩組點(diǎn)云重疊部分會(huì)有明顯的偏差，較先擬合的曲線對(duì)擬合結(jié)果的影響較大。

第二種方法旨在解決第一種方法中可能出現(xiàn)的問題，避免因曲線擬合的順序不同產(chǎn)生較大的偏差，對(duì){p1,k1}和{p2,k2}同時(shí)進(jìn)行擬合。把兩個(gè)控制點(diǎn)向量C1和C2寫成一個(gè)列向量，第r步迭代中，在當(dāng)前曲線f2(t2)上取足夠多的樣點(diǎn)f2(t2,k′),k′=1,…,m′。下一步迭代中由于曲線f2(t2)的變化，f2(t2,k′)也隨之移動(dòng)，記為(t2,k′)，要求(t2,k′)都在(t1)的內(nèi)側(cè)。于是可轉(zhuǎn)化為解問題：

記號(hào)定義同上文。曲線f2(t2)上的樣點(diǎn)f2(t2,k′)可以寫成C2的函數(shù)，故是C1和C2的函數(shù)。對(duì)屬于C1的控制點(diǎn)的偏導(dǎo)可參照式(12)，對(duì)屬于C2的控制點(diǎn)c2,i的偏導(dǎo)：

B2,i(t2)是f2(t2)的B樣條基函數(shù)。于是式(13)的約束可寫成：

多組數(shù)據(jù)的相互約束擬合可基于上述兩組數(shù)據(jù)擬合的約束條件。lk對(duì)不與當(dāng)前擬合曲線有邊界約束關(guān)系的曲線的控制點(diǎn)的偏導(dǎo)均為0，因此整理后也可寫成如(14)的線性不等式約束。

2.2 凸曲線的擬合

根據(jù)燃?xì)獗▽?shí)驗(yàn)的特點(diǎn)，在氧氣濃度固定的情況下，油氣濃度達(dá)到爆炸區(qū)間的某一值時(shí)，爆炸強(qiáng)度可達(dá)到最大，爆炸強(qiáng)度曲線應(yīng)是凸曲線。在上一小節(jié)的的約束擬合中加入新的約束用以保證擬合曲線的上凸性質(zhì)。不妨設(shè)擬合函數(shù)f(t)需要保持向外凸的性質(zhì)，故在曲線的任一點(diǎn)處的曲率中心均在曲線的內(nèi)側(cè)。這一條件也可敘述為f′′(t)始終與f(t)的向內(nèi)法向量同向。沿用上文的記號(hào)，J是將f′(t)向曲線外旋轉(zhuǎn)π/2的矩陣，則JT是將f′(t)向曲線內(nèi)旋轉(zhuǎn)π/2的矩陣。曲線f(t)=L(t)TC應(yīng)滿足：

寫成控制點(diǎn)向量C的約束：

上式的左側(cè)發(fā)生擾動(dòng)δC后在C處展開并線性化：

式(16)就是保持?jǐn)M合曲線外凸性質(zhì)的約束條件。在實(shí)際應(yīng)用中，可以在每一次迭代中在f(t)上取足夠多的樣點(diǎn)，用式(16)的約束求解子問題的最優(yōu)化問題。

上述條件實(shí)現(xiàn)中，算法需要在f(t)上取足夠多的樣點(diǎn)，因此存在如何控制所取樣點(diǎn)數(shù)目的問題。下面我們根據(jù)樣條的變差縮減性質(zhì)直接給出控制點(diǎn)約束，從而導(dǎo)出線性的簡(jiǎn)化條件。設(shè)擬合函數(shù)f(t)是外凸，由B樣條曲線變差縮減性，曲線的凸性受控于其控制多邊形凸性。如前文，J是將向曲線外旋轉(zhuǎn)π/2的矩陣，則JT=-J是向曲線內(nèi)旋轉(zhuǎn)π/2的矩陣。外凸曲線f(t)應(yīng)滿足：

上式的左側(cè)發(fā)生擾動(dòng)δC后在C處展開并線性近似，設(shè)ai=ci+1-ci,i=1,…,n-1，整理后得：

式(18)就是用變差縮減性質(zhì)保證曲線外凸的約束條件，只使用了控制點(diǎn)之間的關(guān)系，比直接用曲線上點(diǎn)的曲率構(gòu)成的式(16)簡(jiǎn)單得多。

2.3 外邊界約束的曲線擬合

油氣爆炸曲線擬合的目的在于對(duì)有爆炸危險(xiǎn)的濃度預(yù)警，在應(yīng)用中可能出現(xiàn)兩種錯(cuò)誤：一種是可燃爆數(shù)據(jù)被判定為安全；另一種則是不可燃爆數(shù)據(jù)被判定為可燃爆。在實(shí)際生產(chǎn)中，前一種錯(cuò)誤顯然更加危險(xiǎn)，一旦誘發(fā)燃爆，帶來的損失將遠(yuǎn)遠(yuǎn)大于后一種錯(cuò)誤。因此，爆炸強(qiáng)度曲線應(yīng)為安全曲線，即至少保證安全曲線之外沒有任何燃爆數(shù)據(jù)點(diǎn)。另一方面，還需考慮預(yù)警的準(zhǔn)確性，即希望擬合曲線盡量靠近實(shí)驗(yàn)數(shù)據(jù)。綜合兩個(gè)要求，需對(duì)爆炸數(shù)據(jù)作外邊界擬合。對(duì)于每一組數(shù)據(jù)的擬合曲線，把這組數(shù)據(jù)作為障礙點(diǎn)約束用于擬合數(shù)據(jù)點(diǎn)的外邊界。數(shù)據(jù)點(diǎn)的邊界擬合是Fl?ry算法的一種特殊情形，通過設(shè)置d(pk,f(t))＜0或d(pk,f(t))＞0即可擬合出的外邊界或內(nèi)邊界。約束條件應(yīng)為：

3 數(shù)據(jù)擬合結(jié)果與分析

通過油氣爆炸實(shí)驗(yàn)中獲取了四組數(shù)據(jù)點(diǎn)，分別記為{p20},{p18},{p16},{p15}，表示實(shí)驗(yàn)反應(yīng)在氧氣濃度為20%、18%、16%、15%采集到的爆炸壓強(qiáng)。用橫坐標(biāo)表示爆炸反應(yīng)前氣體中油氣的濃度，縱坐標(biāo)表示爆炸時(shí)測(cè)得的平均壓強(qiáng)，需要用四條B樣條曲線f20(t),f18(t),f16(t),f15(t)分別擬合這四組爆炸數(shù)據(jù)。首先用B樣條曲線擬合的一般方法分別擬合四組數(shù)據(jù)點(diǎn)，可以粗略地看一下爆炸曲線的走勢(shì)，得到的擬合結(jié)果如圖3。

圖3 一般B樣條擬合

從圖3可以看出，曲線無凸性且兩條曲線有交點(diǎn)，不符合燃?xì)馊急匦?。依?jù)曲線不相交約束的擬合算法，使用第一種方法，先對(duì)其中一組數(shù)據(jù)點(diǎn)做曲線擬合，然后對(duì)其他數(shù)據(jù)點(diǎn)分別做擬合時(shí)，在已擬合好的曲線上取若干點(diǎn)作為邊界約束。用第一種方法得到自上而下的約束擬合圖4(a)和自下而上的約束擬合圖4(b)，相應(yīng)的曲線擬合優(yōu)先順序分別是f20(t)＞f18(t)＞f16(t)＞f15(t)和f15(t)＞f16(t)＞f18(t)＞f20(t)。從圖中可以看出，由于曲線擬合順序的不同，圖4(a)和圖4(b)有不同的擬合效果，這里較明顯的不同表現(xiàn)在曲線f18(t)和f16(t)上。與圖3沒有約束情況下的擬合結(jié)果相比，顯然在圖4(a)中，f18(t)和圖3中相應(yīng)的曲線形狀大致相同，而f16(t)則因?yàn)閒18(t)的限制而有明顯的變化，在圖4(b)中則相反。這一結(jié)果驗(yàn)證了先擬合的曲線對(duì)整體擬合結(jié)果影響較大，而后擬合的曲線對(duì)先擬合的曲線沒有約束效果。

圖4 帶邊界約束的爆炸數(shù)據(jù)擬合

用第二種方法，式(14)，四條擬合曲線相互約束，互為上/下邊界，擬合結(jié)果如圖5。與圖4(a)和圖4(b)的按順序分別擬合的效果相比，圖5的相互擬合在直觀上有折中的效果，尤其在f18(t)和f16(t)上表現(xiàn)得比較明顯。從原理上來說，實(shí)驗(yàn)測(cè)得的數(shù)據(jù)由于一些未知原因產(chǎn)生誤差，因此不能確定哪一組數(shù)據(jù)比較準(zhǔn)確或者靠近理論值。如果較先擬合的數(shù)據(jù)不準(zhǔn)確，那么整體的擬合效果就會(huì)受到較大的影響。因此，采用第二種方法將曲線彼此作為約束條件能夠較好地平衡這種關(guān)系，減小由于測(cè)量誤差導(dǎo)致的擬合誤差的增大。

圖5 互為邊界約束的爆炸數(shù)據(jù)擬合

從理論上看，固定氧氣的濃度，實(shí)驗(yàn)爆炸時(shí)的壓強(qiáng)對(duì)于油氣的濃度作圖會(huì)呈現(xiàn)上凸的形狀。前面的擬合曲線都呈現(xiàn)出波動(dòng)的狀態(tài)，要避免這種波動(dòng)可以加上上凸曲線的約束。如果只用式(18)的約束來做擬合，得到圖6(a)，與圖3類似，曲線之間有交點(diǎn)。圖6(b)在式(18)的基礎(chǔ)上加上式(14)，更加符合理論上對(duì)爆炸壓強(qiáng)擬合曲線走向的預(yù)期。

為了得到爆炸強(qiáng)度的安全曲線，我們對(duì)四組數(shù)據(jù)點(diǎn)做相互約束的外邊界擬合，得到圖7(a)，直觀上是將爆炸點(diǎn)都包在曲線內(nèi)部。但如前討論，曲線呈現(xiàn)出不合理波動(dòng)狀態(tài)。加上了上凸約束之后的圖7(b)與圖6(b)相比，圖7(b)的爆炸強(qiáng)度曲線形狀相似度更高一些。綜合考慮到擬合曲線的準(zhǔn)確性、曲線形狀與真實(shí)情況的比較、以及危險(xiǎn)預(yù)警的安全性，圖7(b)是本文推薦的爆炸強(qiáng)度曲線的擬合結(jié)果，當(dāng)氧氣濃度為15%、16%、18%、20%時(shí)，其對(duì)應(yīng)的油氣燃爆界限分別為[0.98,1.81]、[0.95,1.91]、[0.75,2.53]、[0.72,3.02]。值得注意的是，具有凸性的安全曲線之間沒有交點(diǎn)，因此這組實(shí)驗(yàn)數(shù)據(jù)擬合中，不必添加互為邊界的約束。

圖6 上凸約束的爆炸數(shù)據(jù)擬合

圖7 爆炸強(qiáng)度的安全曲線擬合

4 結(jié)束語

本文針對(duì)某燃?xì)馊急瑢?shí)驗(yàn)數(shù)據(jù)，給出了多約束條件樣條擬合的一個(gè)應(yīng)用實(shí)例。該實(shí)例本質(zhì)上是給定數(shù)據(jù)的多條件曲線擬合問題，本文主要在于約束條件數(shù)學(xué)化、條件合理簡(jiǎn)化以及擬合算法的實(shí)現(xiàn)。經(jīng)過分析和比較，我們選擇具有上凸性質(zhì)的安全曲線擬合為最終分析方案。擬合得到安全曲線和爆炸界限都比較合理，并且得到合作消防公司的認(rèn)可。更為重要的是，該方法可以推廣到多因素燃爆模型。實(shí)際工程生產(chǎn)中，燃爆影響因素除了本文涉及的各種混合氣體濃度和壓強(qiáng)，還可以有溫度、濕度等。添加新的影響因素，在我們的分析方法中轉(zhuǎn)化為增加樣條擬合維數(shù)，即多約束曲線擬合轉(zhuǎn)化為(超)曲面擬合。

[1]Pottmann H,Leopoldseder S,Hofer M.Approximation with active B-spline curves and surfaces [C]//10th Pacific Conference on Computer Graphics and Applications,Proceedings,2002: 8-25.

[2]Shen Liyong.Error Bounded Conic Spline Approximation for NC Code [C]//Proc.of International Conference on Machine Vision (ICMV2011),SPIE 8349,2011: 1-7.

[3]Wang Wenping,Pottman H,Liu Yang.Fitting B-spline curves to point clouds by curvature-based squared distance minimization [J].ACM Transactions on Graphic,2006,25(2): 214-238.

[4]Liu Yang,Wang Wenping.A revisit to least squares orthogonal distance fitting of parametric curves and surfaces [J].Lecture Notes Computer Science,2008,4975: 384-397.

[5]Yang Zhouwang,Deng Jiansong,Chen Falai.Fitting unorganized point clouds with active implicit B-spline curves [J].The Visual Computer,2005,21(8-10):831-839.

[6]Fl?ry S.Fitting curves and surfaces to point clouds in the presence of obstacles [J].Computer Aided Geometric Design,2009,26(2): 192-202.

[7]Rasbash D.Review of explosion and fire hazard of liquefied petroleum gas [J].Fire Safety Journal,1980,2:223-236.

[8]Bakke J R,Van Wingerden K,Hoorelbeke P,Brewerton B.A study on the effect of trees on gas explosions [J].Journal of Loss Prevention in the Process Industries,2010,23(6): 878-884.

[9]Cheng Jianwei,Yang Shengqiang.Improved coward explosive triangle for determining explosibility of mixture gas [J].Process Safety and Environmental Protection,2011,89(2): 89-94.

[10]Chang J I,Lin Chengchung.A study of storage tank accidents [J].Journal of Loss Prevention in the Process Industries,2006,19: 51-59.

[11]CCPS.Guidelines for Chemical Process Quantitative Risk Analysis [M].2nd ed.Wiley-AIChE,1999: 57-283.

[12]Vapnik V.The Nature of Statistical Learning Theory [M].NY,Springer,1999: 137-170.

[13]王濟(jì)川,郭志剛.Logistic回歸模型——方法與應(yīng)用[M].北京: 高等教育出版社,2001: 146-181.

[14]Brunnet G,Hagen H,Santarelli P.Variational design of curves and surfaces [J].Surveys on Mathematics for Indestry,1993,3: 1-27.