王軍磊 賈愛林 甯波 蔣俊超

中國石油勘探開發(fā)研究院

對于頁巖氣、致密氣等非常規(guī)氣藏，水平井輔以分段壓裂技術(shù)能夠增大地層接觸面積，減小滲流阻力，提高氣井產(chǎn)能，同時(shí)也將引起持續(xù)數(shù)年的不穩(wěn)定線性流動期［1－4］。將日常生產(chǎn)數(shù)據(jù)處理為產(chǎn)量修正的生產(chǎn)擬壓差、擬時(shí)間，以反映出的流態(tài)特征為診斷工具，利用相應(yīng)的數(shù)學(xué)表達(dá)式能夠獲取不同的參數(shù)，以此來評價(jià)壓裂效果與預(yù)測氣井產(chǎn)量［5－6］。

但是，氣體擬時(shí)間計(jì)算中涉及的平均地層壓力往往難以確定。Agrwal［7］和 Mattar［8］分別利用壓力迭代法和地質(zhì)儲量迭代法確定平均地層壓力，但都只適用于擬穩(wěn)態(tài)生產(chǎn)階段；Anderson［9］為提出使用探測邊界內(nèi)平均地層壓力的觀點(diǎn)，隨后 Nobakht等［10－11］給出了利用物質(zhì)平衡方程結(jié)合探測邊界傳播規(guī)律的求解思路，該方法雖然突破了擬穩(wěn)態(tài)流動條件的限制，但難以應(yīng)用到更具實(shí)際意義的變產(chǎn)量生產(chǎn)情況。

這樣，如何計(jì)算變產(chǎn)量下的探測邊界移動規(guī)律成為合理分析生產(chǎn)數(shù)據(jù)的關(guān)鍵。傳統(tǒng)探測邊界公式是基于不穩(wěn)定滲流研究的結(jié)果，多是通過脈沖波的最大響應(yīng) 位 置［12－13］或 聯(lián) 立 不 穩(wěn) 定 與 擬 穩(wěn) 態(tài) 壓 力 導(dǎo) 數(shù)［14－15］確定。這些公式形式基本一致，均認(rèn)為探測邊界移動速度與氣井工作制度無關(guān)，但 Wattenbarger［16］計(jì)算表明，定壓和定產(chǎn)下的探測邊界公式并不相同。鑒于此，筆者借助低滲透儲層中的“動邊界”［17－19］概念，將探測邊界視為壓力擾動的外邊緣，通過積分方法獲得變產(chǎn)量下的探測邊界通用表達(dá)式，利用Nobakht方法計(jì)算平均地層壓力，并給出相應(yīng)迭代算法計(jì)算氣體擬時(shí)間函數(shù)，同時(shí)對比分析使用真實(shí)時(shí)間對線性關(guān)系產(chǎn)生的影響，進(jìn)而利用擬時(shí)間分析氣井不穩(wěn)定生產(chǎn)數(shù)據(jù)，計(jì)算壓裂參數(shù)。

1 探測邊界內(nèi)的平均地層壓力

1.1 探測邊界模型及近似解

致密介質(zhì)滲透率極低，無自然產(chǎn)能，氣體只有在壓裂區(qū)域（SRV）內(nèi)才能夠發(fā)生有效流動。假設(shè)各條裂縫等長、等間距分布，SRV區(qū)域內(nèi)的滲流可以等效為一系列的線性流模型［20－22］。

借助低滲透“動邊界”概念研究探測邊界問題：假定地層均質(zhì)，含氣飽和度100%，氣體全部為CH4，忽略重力和毛細(xì)管力的影響，不考慮地層孔隙的可壓縮性，壓力只在探測邊界內(nèi)傳播。利用無量綱量處理氣體擬壓力控制方程，其中變產(chǎn)量下無量綱定義為：

式中K為滲透率，mD；h為地層厚度，m；φ為地層孔隙度；μg為氣體黏度，mPa·s；Bg為氣體體積系數(shù)；cg為氣體壓縮系數(shù)，MPa－1；Zg為氣體偏差因子；T為地層溫度，K；y為空間坐標(biāo)，m；yf為動邊界空間位置，m；xf為裂縫半長，m；qsc為氣井產(chǎn)量，m3／d；qref為參考產(chǎn)量，m3／d；Gp為累積產(chǎn)量，m3；p為地層壓力，MPa；t為生產(chǎn)時(shí)間，h；下標(biāo)sc為標(biāo)準(zhǔn)狀態(tài)，psc＝0.1MPa，Tsc＝293.15K；下標(biāo)i為地層原始狀態(tài)；下標(biāo)D為無量綱量。

地層中的無量綱擬壓力控制方程滿足：

壓力及探測邊界初始條件為：

由于探測邊界yfD隨時(shí)間而變化，直接求解（2）～（4）并非易事，這里采用近似求解方法。首先對控制方程（2）兩側(cè)先后關(guān)于空間變量（yD）和時(shí)間變量（tD）進(jìn)行積分，利用邊界條件將偏微分方程式（2）轉(zhuǎn)化為探測邊界內(nèi)物質(zhì)平衡方程：

考慮到初始條件式（3）及邊界條件式（4），令方程式（5）的近似解滿足二階精度近似，即

上式中α0（tD）和α1（tD）為待定系數(shù)，通過邊界條件可以得到探測邊界內(nèi)的壓力分布公式，即

將壓力近似解式（7）代入物質(zhì)平衡方程式（5），可以得到探測邊界在地層中的移動規(guī)律：

式（8）與經(jīng)典探測邊界公式相比，系數(shù)更大、時(shí)間修正為物質(zhì)平衡時(shí)間。利用式（8）可以得到任意生產(chǎn)制度下氣井探測邊界隨時(shí)間的移動規(guī)律。

1.2 探測邊界內(nèi)平均地層壓力

從氣井流量來源的角度分析探測邊界的物理意義。氣井流量完全來自于地層孔隙內(nèi)氣體的彈性壓縮，以探測邊界作為空間劃分點(diǎn)，氣井流量由兩部分地層組成：探測邊界內(nèi)＋探測邊界外。使用物質(zhì)平衡方程的理想條件為外邊界封閉，即探測邊界內(nèi)地層的流量供給占?xì)饩髁康?00%，所以探測邊界以內(nèi)地層對氣井產(chǎn)量的供給比例決定著在探測邊界內(nèi)使用物質(zhì)平衡方程的準(zhǔn)確性。

本文參考文獻(xiàn)［16］給出了定流量探測邊界表達(dá)式y(tǒng)Df＝（2tD）1／2，定壓為yDf＝（4tD）1／2，而利用式（8）得到的探測邊界定流量表達(dá)式為：yDf＝（6tD）1／2。定壓：yDf＝（6tDmb）1／2＝（12tD）1／2。利用經(jīng)典解關(guān)于yD的導(dǎo)數(shù)可獲得不同位置處的流量qD（yD，tD），通過計(jì)算qD（yD，tD）／qD（0，tD）比值獲得探測邊界內(nèi)地層對氣井流量的供給比例。其中本文參考文獻(xiàn)［16］定產(chǎn)條件下的計(jì)算結(jié)果為68.3%，定壓為63.9%，近似解式（8）分別為91.68%、95.02%。使用近似解增加了壓力擾動的波及范圍，減小了探測邊界外地層的流量供給比重，提高了在探測邊界內(nèi)使用物質(zhì)平衡方程的精度。

利用式（8）并結(jié)合物質(zhì)平衡方程式（9）可以得到不同時(shí)刻探測邊界內(nèi)的平均地層壓力（pavg），這是計(jì)算氣體擬時(shí)間函數(shù)的基礎(chǔ)。

其中，利用式（8）可獲得不同生產(chǎn)制度下的探測邊界內(nèi)地質(zhì)儲量（G）為：

2 氣體擬時(shí)間

氣體具有強(qiáng)可壓縮性，其黏度（μg）、偏差因子（Zg）、體積系數(shù)（Bg）等都是關(guān)于壓力的強(qiáng)非線性函數(shù)，如果直接對氣體滲流控制方程中的擴(kuò)散系數(shù)進(jìn)行強(qiáng)行近似必然會引起較大的誤差。使用擬時(shí)間函數(shù)（ta）能夠?qū)怏w滲流問題等效轉(zhuǎn)化為液體滲流問題，同時(shí)改進(jìn)了“強(qiáng)行近似擴(kuò)散系數(shù)（1／φμgcg）為常數(shù)”的假設(shè)條件，結(jié)果必然使得理論描述更接近于礦場實(shí)際情況，也方便液體研究領(lǐng)域內(nèi)成果的參考和借鑒。其中，擬時(shí)間定義為：

2.1 定壓生產(chǎn)

在定壓生產(chǎn)條件下，線性流動階段的產(chǎn)量變化規(guī)律［16］滿足1／qsc＝αt0.5。其中：

式（13）表明探測邊界內(nèi)的平均地層壓力為常數(shù)，擬時(shí)間與真實(shí)時(shí)間關(guān)系簡化為：

用擬時(shí)間取代真實(shí)時(shí)間，產(chǎn)量和探測邊界內(nèi)的平均地層壓力重新修正式為：

其中，擬時(shí)間修正因子用來表示真實(shí)時(shí)間與擬時(shí)間的差異，定義為：

根據(jù)式（16）結(jié)合Newton迭代算法可以獲得不同生產(chǎn)壓差下的平均地層壓力。圖1計(jì)算了不同井底壓力（pw）下對應(yīng)的擬時(shí)間修正因子（fcp）。結(jié)果表明，pw越小，fcp越大，（mi－mw）／qsc—t斜率較（mi－mw）／qsc—ta偏小程度越明顯，直接使用真實(shí)時(shí)間分析生產(chǎn)數(shù)據(jù)的可靠性越差。

圖1 不同井底壓力下的擬時(shí)間修正因子圖

2.2 定產(chǎn)生產(chǎn)

定產(chǎn)生產(chǎn)時(shí)，同樣利用式（9）式（10）可以獲得探測邊界內(nèi)的物質(zhì)平衡方程：

為方便研究，這里做如下假設(shè)：cg≈1／p，Zg＝1，μg＝常數(shù)，結(jié)合利用式（18），可以得到擬時(shí)間與真實(shí)時(shí)間的近似關(guān)系式為：

則產(chǎn)量修正下的擬壓差（mi－mw）／qsc與真實(shí)時(shí)間t的線性關(guān)系應(yīng)滿足：

用圖2計(jì)算得到不同產(chǎn)量下的t1／2與［m（pi）－m（pw）］／qsc變 化 規(guī) 律。 結(jié) 果 表 明：t1／2與 ［m（pi）－m（pw）］／qsc在生產(chǎn)早期呈線性關(guān)系，隨著生產(chǎn)進(jìn)行逐漸偏離直線關(guān)系（線性相關(guān)系數(shù)變?。x程度受氣井產(chǎn)量控制，產(chǎn)量越大偏離直線的起始時(shí)間越小，偏離程度越大（線性相關(guān)系數(shù)越小），這將導(dǎo)致無法直接使用真實(shí)時(shí)間分析氣井生產(chǎn)數(shù)據(jù)。

圖2 產(chǎn)量修正擬壓差與生產(chǎn)時(shí)間的線性關(guān)系圖

2.3 變產(chǎn)量生產(chǎn)

物質(zhì)平衡時(shí)間（tmb）能夠較好地處理產(chǎn)量變化引起的時(shí)間疊加影響，但Palacio和Blasingame［23］研究表明，tmb通常只在擬穩(wěn)定階段精確成立。將物質(zhì)平衡時(shí)間（tmb）修正為1.23tmb，物質(zhì)平衡時(shí)間在不穩(wěn)定流動階段也能夠精確成立（表1）。

表1 修正物質(zhì)平衡時(shí)間效果表

將式（20）中的擬時(shí)間（ta）替換為修正物質(zhì)平衡擬時(shí)間1.23tmba，可得到變產(chǎn)量條件下的解析表達(dá)式為：

其中，探測邊界內(nèi)的平均地層壓力通過下式確定，即

式（21）與式（23）中均含有未知量xf，這里采用迭代方法計(jì)算xf：

1）繪制［m（pi）－m（pw）］／qsc與真實(shí)物質(zhì)平衡時(shí)間（tmb1／2）的曲線，利用式（21）的線性關(guān)系式中的斜率（β）計(jì)算xf，作為初值。

2）使用式（23）計(jì)算生產(chǎn)數(shù)據(jù)記錄點(diǎn)處的平均地層壓力，形成t—pavg的數(shù)據(jù)表。

3）利用t—pavg數(shù)據(jù)表并結(jié)合式（22），借助數(shù)值積分計(jì)算物質(zhì)平衡擬時(shí)間（tmba），進(jìn)一步形成t—tmba數(shù)據(jù)表。

4）繪制各個(gè)時(shí)間點(diǎn)t對應(yīng)的［m（pi）－m（pw）］／qsc與物質(zhì)平衡擬時(shí)間（tmba1／2）的線性關(guān)系式，確定斜率（β），計(jì)算xf。

5）重復(fù)步驟（2）～（4），直到xf收斂為止。

3 實(shí)例分析

以某壓裂水平氣井為例進(jìn)行計(jì)算驗(yàn)證。其中地層滲透率為 0.187mD，厚度為 11.7m，孔隙度為12.4%，原始地層壓力為28.95MPa，原始地層溫度為313.15K，水平段長度為1045m，壓裂10段，氣井生產(chǎn)歷史見圖3。

圖3 氣井生產(chǎn)數(shù)據(jù)圖

按照上述迭代算法分析氣井生產(chǎn)數(shù)據(jù)。圖4計(jì)算了不同時(shí)刻探測邊界移動規(guī)律以及探測邊界內(nèi)的平均地層壓力，結(jié)果表明：隨著生產(chǎn)的進(jìn)行，探測邊界不斷向外傳播，傳播速度逐漸減慢，同時(shí)井底壓力不斷下降，探測邊界內(nèi)的平均地層壓力隨之下降。

圖4 探測邊界及平均地層壓力變化規(guī)律圖

圖5對比了使用物質(zhì)平衡時(shí)間和物質(zhì)平衡擬時(shí)間的計(jì)算效果：①直接使用物質(zhì)平衡時(shí)間，生產(chǎn)數(shù)據(jù)間的線性關(guān)系并不明顯，直線的線性相關(guān)系數(shù)較低，斜率整體偏低，同時(shí)呈現(xiàn)出定壓和定產(chǎn)情況的特征；②使用物質(zhì)平衡擬時(shí)間，能夠明確生產(chǎn)數(shù)據(jù)間的線性關(guān)系，此時(shí)擬合出的線性關(guān)系式是考慮了氣體高壓物性和產(chǎn)量變化共同作用的結(jié)果，更接近氣體的實(shí)際流動情況。利用生產(chǎn)數(shù)據(jù)間的線性關(guān)系斜率（β）來計(jì)算壓裂長度（xf）：

圖5 利用擬時(shí)間分析生產(chǎn)數(shù)據(jù)結(jié)果圖

由于圖5中物質(zhì)平衡擬時(shí)間對應(yīng)的直線斜率（βtmba）大于物質(zhì)平衡時(shí)間斜率（βtmb），故使用真實(shí)時(shí)間分析氣井生產(chǎn)數(shù)據(jù)導(dǎo)致壓裂參數(shù)計(jì)算結(jié)果偏高，物質(zhì)平衡時(shí)間計(jì)算結(jié)果為xf＝53.40m，物質(zhì)平衡擬時(shí)間計(jì)算結(jié)果為41.73m，為進(jìn)一步驗(yàn)證計(jì)算結(jié)果，基于線性流動模型重新計(jì)算井底壓力，同時(shí)對比實(shí)測壓力（圖6）。圖6表明，基于擬時(shí)間（xf＝41.73m）解釋參數(shù)的預(yù)測結(jié)果更為合理，而真實(shí)時(shí)間（xf＝53.40m）解釋結(jié)果則高估了氣井的實(shí)際生產(chǎn)能力，導(dǎo)致相同生產(chǎn)時(shí)間內(nèi)井底壓力的下降幅度更大。

圖6 生產(chǎn)數(shù)據(jù)擬合效果圖

4 結(jié)論

1）完整給出了考慮氣井產(chǎn)量變化和探測邊界影響的氣體滲流數(shù)學(xué)模型，用積分平均方法完成了近似求解，得到了探測邊界傳播規(guī)律的通用解析表達(dá)式，并通過對比經(jīng)典結(jié)果驗(yàn)證近似解的準(zhǔn)確性。

2）探測邊界公式中時(shí)間修正為物質(zhì)平衡時(shí)間，利用新公式計(jì)算的探測邊界內(nèi)地層儲存氣體對氣井產(chǎn)量貢獻(xiàn)率超過90%，高于 Watterbarger結(jié)果，提高了在探測邊界內(nèi)使用物質(zhì)平衡方程的精度。

3）同修正擬壓差—擬時(shí)間形成的線性關(guān)系式相比，定產(chǎn)條件下直接使用真實(shí)時(shí)間導(dǎo)致線性關(guān)系不成立，定壓條件下則引起直線斜率減小，其偏差程度分別由氣井產(chǎn)量及井底壓力決定。

4）在實(shí)際氣井生產(chǎn)中，氣體物質(zhì)平衡擬時(shí)間函數(shù)能夠明確生產(chǎn)數(shù)據(jù)間的線性關(guān)系，提高壓裂參數(shù)解釋結(jié)果的可靠性，為準(zhǔn)確評價(jià)水力壓裂效果提供理論支持。

［1］曾凡輝，王樹義，郭建春，等.裂縫面非均質(zhì)流入的氣藏壓裂水平井產(chǎn)量計(jì)算［J］.天然氣工業(yè)，2014，34（5）：100－105.ZENG Fanhui，WANG Shuyi，GUO Jianchun，et al.Yield calculation of a fractured horizontal well with a non－uniform gas flow on fracture surface［J］.Natural Gas Industry，2014，34（5）：100－105.

［2］CLARKSON C R.Production data analysis of unconventional gas wells：Review of theory and best practices［J］.International Journal of Coal Geology，2013：101－146.

［3］KABIR C S，ISMADI D，F(xiàn)OUTAIN S.Estimating in－place volume and reservoir connectivity with real－time and periodic surveillance data［J］.Journal of Petroleum Science and Engineering，2011，78：258－266.

［4］張小濤，吳建發(fā)，馮曦，等.頁巖氣藏水平井分段壓裂滲流特征數(shù)值模擬［J］.天然氣工業(yè)，2013，33（3）：47－52.ZHANG Xiaotao，WU Jianfa，F(xiàn)ENG Xi，et al.Numerical simulation of seepage flow characteristics of multi－stage fracturing in horizontal shale gas wells［J］.Natural Gas Industry，2013，33（3）：47－52.

［5］ANDERSON D M，LIANG P，OKOUMA V，et al.Probabilistic forecasting of unconventional resources using rate transient analysis：Case studies［C］∥paper 155737－MS presented at the Americas Unconventional Resources Conference，5－7June 2012，Pittsburgh，Pennsylvania，USA.New York：SPE，2012.

［6］NOBAKHT M，MATTAR L，MOGHADAM S，et al.Simplified forecasting of tight／shale－gas production in linear flow［J］.Journal of Canadian Petroleum Technology，2012，11：476－486.

［7］AGRWAL R G，GARDNER D C，KLEINSTEIBER S W，et al.Analyzing well production data using combined type curve and decline curve concepts［J］.SPE Reservoir Evaluation ＆ Engineering，1999，2（5）：478－486.

［8］MATTAR L，ANDERSON D M.A systematic and comprehensive methodology for advanced analysis of production data［C］∥paper 84472－MS presented at the SPE Annual Technical Conference and Exhibition，5－8October 2003，Denver，Colorado，USA.New York：SPE，2003.

［9］ANDERSON D M，MATTAR L.An improved pseudotime for gas reservoirs with significant transient flow［J］.Journal of Canadian Petroleum Technology，2007，46（7）：49－54.

［10］NOBAKHT M，CLARKSON C R.A new analytical method for analyzing linear flow in tight／shale gas reservoirs：Constant－flowing－pressure boundary condition［J］.SPE Reservoir Evaluation ＆ Engineering，2012，7：370－384.

［11］NOBAKHT M，CLARKSON C R.A new analytical method for analyzing linear flow in tight／shale gas reservoirs：Constant－rate boundary condition［J］.SPE Reservoir Evaluation ＆ Engineering，2012，2：51－59.

［12］崔迪生，賀子倫，郭寶健，等.徑向復(fù)合油藏探測半徑計(jì)算方法［J］.油氣井測試，2005，14（1）：15－18.CUI Disheng，HE Zilun，GUO Baojian，et al.Calculating method for investigation radius in radial compound reservoirs［J］.Well Testing，2005，14（1）：15－18.

［13］王曉冬.滲流力學(xué)基礎(chǔ)［M］.北京：石油工業(yè)出版社，2006：82.WANG Xiaodong.Fundamental mechanics of fluid flow in porous media［M］.Beijing：Petroleum Industry Press，2006：82.

［14］李傳亮.油井探測半徑的計(jì)算公式研究［J］.大慶石油地質(zhì)與開發(fā)，2002，21（5）：32－34.LI Chuanliang.Study of calculated equation of investigation radius of production well［J］.Petroleum Geology ＆Oilfield Development in Daqing，2002，21（5）：32－34.

［15］齊麗巍，王曉冬.探測半徑計(jì)算方法研究［J］.油氣井測試，2007，16（2）：1－4.QI Liwei，WANG Xiaodong.Study on determination of the radius of investigation［J］.Well Testing，2007，16（2）：1－4.

［16］WATTENBARGER R A，EI－BANBI A H，VILLEGAS M E，et al.Production analysis of linear flow into fractured tight gas wells［C］∥paper 39931presented at the SPE Rocky Mountain Region／Low－Permeability Reservoirs Symposium，5－8April 1998，Denverm，Colorado，USA.New York：1998.

［17］李愛芬，劉艷霞，張華強(qiáng)，等.用逐步穩(wěn)態(tài)替換法確定低滲透油藏合理井距［J］.中國石油大學(xué)學(xué)報(bào)：自然科學(xué)版，2011，35（1）：89－93.LI Aifen，LIU Yanxia，ZHANG Huaqiang，et al.Determination of reasonable spacing in low permeability reservoirs by means of series of steady state analysis［J］.Journal of China University of Petroleum：Natural Science Edition，2011，35（1）：89－93.

［18］孫來喜，張宗輝，王仕莉，等.低滲透油藏井間儲量動用狀況分析新方法［J］.西南石油大學(xué)學(xué)報(bào)：自然科學(xué)版，2014，36（2）：99－104.SUN Laixi，ZHANG Zonghui，WANG Shili，et al.A new method of analysis of low permeability reservoir interwell reserves producing status［J］.Journal of Southwest Petroleum University：Science ＆ Technology Edition，2014，36（2）：99－104.

［19］王曉冬，侯曉春，郝明強(qiáng)，等.低滲透介質(zhì)有啟動壓力梯度的不穩(wěn)態(tài)壓力分析［J］.石油學(xué)報(bào)，2011，32（5）：843－847.WANG Xiaodong，HOU Xiaochun，HAO Mingqiang，et al.Pressure transient analysis in low－permeability media with threshold gradients［J］.Acta Petrolei Sinica，2011，32（5）：843－847.

［20］趙海洋，鄔藍(lán)柯西，劉青山，等.不同完井方式下水平井不穩(wěn)定產(chǎn)能研究［J］.西南石油大學(xué)學(xué)報(bào)：自然科學(xué)版，2012，34（5）：133－136.ZHAO Haiyang，WULAN kexi，LIU Qingshan，et al.A study on difference completion methods of horizontal well productivity［J］.Journal of Southwest Petroleum University：Science ＆ Technology Edition，2012，34（5）：133－136.

［21］曾保全，程林松，齊梅，等.特低滲透油藏壓裂水平井新型布井方式研究［J］.西南石油大學(xué)學(xué)報(bào)：自然科學(xué)版，2013，35（2）：115－120.ZENG Baoquan，CHENG Linsong，QI Mei，et al.Novel well pattern for fractured horizontal well in ultra－low permeability reservoir［J］.Journal of Southwest Petroleum U－niversity：Science ＆ Technology Edition，2013，35（2）：115－120.

［22］劉曉旭，楊學(xué)鋒，陳元林，等.頁巖氣分段壓裂水平井滲流機(jī)理及試井分析［J］.天然氣工業(yè)，2013，33（12）：77－81.LIU Xiaoxu，YANG Xuefeng，CHEN Yuanlin，et al.Seepage behavior and well testing in horizontal shale gas wells under multi－stage fracking［J］.Natural Gas Industry，2013，33（12）：77－81.

［23］PALACIO J C，BLASINGAME T A.Decline－curve analysis using type curves－analysis of gas well production data［C］∥paper 25909－MS presented at the SPE Joint Rocky Mountain and Low Permeability Reservoirs Symposium，26－28April 1993，Denver，Colorado，USA.New York：SPE，1993.