Petrochem1.0巖石化學軟件及其應用

丁衛(wèi)平

中化地質礦山總局地質研究院，河北涿州，072754

0 引言

隨著地質勘查的大面積展開，屬于基礎地質工作的巖石化學應用到各類巖漿巖和火山巖的研究顯得很常見，這方面雖有類似的軟件來解決，但都各有不足之處，比如某些功能的缺陷，作者利用既熟悉軟件開發(fā)工作，又熟悉地質方面知識的優(yōu)勢，通過緊密結合中國地質調查局地質勘探行業(yè)標準開發(fā)了Petrochem1.0巖石化學軟件，除了常用的功能外，還補充了地球化學中有關的部分，比如鈦鐵銅鎳鉻的成礦可能性方面。經過了近一年的黑白盒測試，已經完全能應用到有關巖石化學方面的數(shù)據處理中，為地質研究提供了很好的基礎。

1 主要功能介紹

巖漿巖尤其是火山巖中的氧化物含量對巖石定名、成礦可能性、構造環(huán)境分析和火山巖系列劃分等是很重要的基本數(shù)據，這些常用的氧化物通常為 SiO2、Al2O3、Fe2O3、FeO、MgO、CaO、Na2O、K2O、H2O、TiO2、P2O5、CO2、MnO、NiO、ZrO2、Cr2O3。這里利用主界面進行功能介紹，以五個大類來進行分別說明。主界面及主要功能如圖1所示，該圖是為了顯示所有主要功能而拼出顯示，右鍵菜單只對選出的數(shù)據（背景變深）起作用，總計有15項大的功能，其中彈出式菜單中有14項功能，第15項是對表中的所有樣品（不必選擇樣品）計算標準礦物和巖石化學指數(shù)（CIPW及指數(shù)），用按鈕來彈出對話窗，這些大的功能根據情況可以有更詳細的次級功能。

2 巖石化學計算方法

巖石化學研究包括全巖特征數(shù)值（扎氏法）、標準礦物和巖石化學指數(shù)計算等三方面，本軟件中的全巖特征數(shù)值法采用的是常用的扎瓦里茨基特征數(shù)值法，而標準礦物計算則采用的是美國的CIPW標準礦物計算法，現(xiàn)分別說明：

2.1 扎瓦里茨基特征數(shù)值法【1】

該方法實際上是求出若干個特征數(shù)據來表示巖漿巖的化學特征，并用圖解表示出來，選擇該項菜單后，彈出扎瓦里茨基圖，將所有選擇的巖石序號顯示在圖2中，至于具體的特征數(shù)據可在“圖信息”中彈出窗口來顯示。

圖 1 巖石氧化物數(shù)據及功能菜單分布情況圖Fig.1 Oxid data and function menu dialog for rock analysis

圖 2 扎瓦里茨基圖解Fig.2 Zavaritskii diagrammatizing

該圖解不能用于超基性巖的投圖，對于該圖如何理解特別設計了一項“讀圖法”按鈕，點擊即可彈出讀圖說明，所有圖形可被保存用于放到需要的文件中，可以將圖形顯示部分硬拷貝輸出成一個bmp格式圖形文件，以供被選用（以后圖形界面有此功能的話，不再介紹）。圖2中的1號樣為內蒙西烏旗的花崗斑巖樣，17號樣為玄武安山巖樣，27號樣為輝長巖樣。

2.2 CIPW 標準礦物計算【1】

在圖1中的下方選擇“CIPW”按扭后，即可彈出圖3所示的對話窗，顯示的標準礦物如下：磷灰石(Ap)、鈦鐵礦(Ilm)、榍石(Tn)、鉀長石(Or)、鈉長石(Ab)、鈣長石(An)、赤鐵礦(Hm)、磁鐵礦(Mt)、透輝石(Di)、硅灰石(Wo)、紫蘇輝石(Hy)、石英(Q)、橄欖石(Ol)、碳酸鈉(Ne)、白榴石(Lc)、剛玉(C)、鉻鐵礦(Cm)、霓石(Ac)、硅酸鈉(Ns)，該窗口也有右鍵菜單來對標準礦物進行進一步的操作。

2.3 巖石化學指數(shù)【2】

巖石化學指數(shù)是判斷巖性特征的重要參數(shù)，包括里特曼指數(shù)(σ)、長英指數(shù)(FL)、固結指數(shù)(SI)、分異指數(shù)(DI)、鎂鐵指數(shù)(MF)、拉森指數(shù)(LI)、結晶指數(shù)(CI)、風化指數(shù)(WI)、氧化度(OXD)、氧化率(OX)、鎂鐵比值(M/F)、鎂鐵比值(m/f)、堿度率(AR)、ALK、K2O/Na2O、A/CNK、戈蒂里指數(shù)、顏色指數(shù)、斜長石號碼等。這些指數(shù)直接與標準礦物一同顯示在同一窗口中。對這些指數(shù)可在圖3所示對話框的“指數(shù)說明”按鈕中彈出圖4所示窗口，了解指數(shù)的具體含義。說明這些指數(shù)一般用在什么情況下，方便隨時查詢。

3 巖石的命名及相關圖解

巖石的命名圖解采用先大類命名，然后鼠標點按不同區(qū)域彈出不同窗口進行更詳細的定名方法。

圖 3 CIPW標準礦物及巖石化學參數(shù)圖表Fig.3 CIPW standardize minerals and index dialog for petrochemistry

圖4 巖石化學指數(shù)說明窗Fig.4 Index-explaining dialog for petrochemistry

3.1 全堿-二氧化硅分類(TAS)圖解【3】

該圖解是中國地質調查局所確定的國家標準圖解，彈出對話窗如圖5，可以依據SiO2、K2O、Na2O投影到這個圖中，看點落在什么區(qū)域就基本確定什么巖性名，再根據每個區(qū)域不同的功能菜單再繼續(xù)分類，該圖分上下兩部分，對鎂含量低者，用上部分定名，對高鎂者，則用下部分命名。圖5中，T區(qū)能右鍵彈出“過堿性粗面巖的分類”菜單，R區(qū)可彈出“過堿性流紋巖的分類”菜單，B區(qū)可彈出“玄武巖分類”菜單，S1+S2+S3區(qū)可彈出“粗面玄武安山巖分類”菜單。此外，B區(qū)中SiO2含量大于48%的區(qū)域和O1、O2、O3及R區(qū)都能彈出“高鉀-中鉀-低鉀分類”菜單。下面一一介紹這些菜單彈出界面。

圖 5 TAS分類圖解Fig.5 TAS classification diagram

3.1.1 過堿性粗面巖的分類【3】T區(qū)菜單——該菜單彈出以下對話窗（圖6）：圖中1、2號樣為內蒙古梅勒圖組鈉閃堿流質粗面巖投影圖。

圖6 粗面巖的分類圖Fig.6 Trachyte Classification diagram

3.1.2 按鉀質含量對玄武巖、安山巖、流紋巖等繼續(xù)分類——高鉀-中鉀-低鉀分類【3】彈出對話窗如下圖7。

圖7 高、中、低鉀熔巖分類圖Fig.7 High-middle-low potassium lava classification diagram

同樣的巖石類型，從陸殼俯沖帶到板塊內部，K2O的含量是逐漸升高的，使用此圖可以做輔助的判斷。

3.1.3 流紋巖區(qū)“過堿性流紋巖”的進一步分類【3】同樣是依據 Al2O3和全鐵的含量進行劃分，圖 8是內蒙古西烏旗的火山巖數(shù)據（3號和6號樣），說明這些巖石都屬于鈉閃堿流質流紋巖。

圖8 Al2O3-FeO(全鐵)圖解Fig.8 Al2O3-（FeO+Fe2O3）(total Fe) diagram

3.1.4 粗面玄武安山巖的分類【3】彈出這菜單項，以便對玄武巖、玄武粗安巖、粗安巖等繼續(xù)按Na2O、K2O含量進一步分類，顯示的窗口如圖9。

圖 9 粗面玄武安山巖的分類命名Fig.9 Classification denomination for trachybasaltic andesite

看投影點落在什么格中就是什么類型的巖石，多個點落入同一格中的話就從左往右，從上往下排列。

3.1.5 玄武巖分類【3】在TAS圖解中的B區(qū)，右鍵選擇該菜單后依據計算的CIPW標準礦物，可以將該區(qū)玄武巖繼續(xù)分類如圖10，可以將玄武巖分出11種詳細類型。B代表Basalt，即玄武巖。“幫助”按鈕彈出對話窗說明各帶括號的數(shù)學表示是什么類型玄武巖，如碧玄巖、橄欖玄武巖、苦橄巖等，圖10中的17號樣表示石英拉斑玄武巖，該樣是內蒙古西烏旗一帶的玄武巖數(shù)據。

圖 10 玄武巖類巖石的CIPW標準礦物命名圖(OQNH)圖Fig.10 Standardize mineral denomination diagram for basalt(OQNH)

3.2 標準礦物的巖石分類定名

這是在計算出CIPW標準礦物的基礎上對一些巖石進行詳細定名。

3.2.1玄武巖類型【3】對選擇的數(shù)據，首先進行過濾，根據SiO2含量去掉不是玄武巖的數(shù)據，然后利用CIPW標準礦物含量以及Al2O3的百分含量來確定進一步的玄武巖分類名字（圖11），依據標準是中國地調局的巖石分類和命名方案——火成巖巖石分類和命名方案 GB/T 17412.1-1998。

圖11 由CIPW標準礦物確定的玄武巖分類Fig.11 Basalt Classification from the CIPW standardize minerals

3.2.2 石英-鉀長石-斜長石-副長石(QAPF)圖解【3】依據計算的標準礦物，投影到下面的圖中，如果是侵入巖見圖12。

圖12 侵入巖分類圖解Fig.12 Classification diagram for Invasive rock

點擊其中的“幫助”，可以得到各分區(qū)的侵入巖名稱說明，注意該名稱與實際可以不一樣，只要差別不是太大都可以，大多是用在見不到具體礦物的火山巖相互之間的對比上。圖所示野外定名為花崗斑巖，實際所在區(qū)域（4）為花崗閃長巖區(qū)。

如果選擇的是火山巖數(shù)據則讀下面的圖13。

圖13 火山巖分類圖解Fig.13 Classification diagram for Volcanic rock

點按其中的“幫助”按鈕，可以得到各分區(qū)的火山巖名稱說明。在圖13的分區(qū)9、10中，玄武巖與安山巖都落在該區(qū)，如何進一步劃分，依據該區(qū)右鍵菜單“玄武安山巖分類”【5】，彈出以下對話窗（圖 14），從而把二者區(qū)分開來，通過下圖說明實際上是安山巖，而實際樣品定名為玄武巖。

圖14 玄武巖和安山巖的進一步分類圖Fig.14 Classification diagram for Volcanic rock

3.2.3 輝長巖分類圖解【3】如果選擇的是輝長巖樣，可以利用圖15將輝長巖進行繼續(xù)詳細定名。點按其中“幫助”按鈕，可以了解各分區(qū)輝長巖的詳細定名以及三角頂點名（輝石、斜長石和橄欖石）。

圖15 輝長巖分類圖解Fig. 15 Classification diagram for Gabbro

圖中的18、19號樣鏡下鑒定為輝長巖，實際投影區(qū)域（2）也為輝長巖

4 巖石形成環(huán)境分析

形成環(huán)境包括大地構造環(huán)境，也包括巖漿形成深度是地殼還是地幔，巖漿是地殼重熔的還是地幔融化的等等。

4.1 里特曼-戈蒂里圖【4】

主要是對火山熔巖進行分析，確定熔巖的產生環(huán)境特點。利用 SiO2、K2O、Na2O的百分含量，計算里特曼指數(shù)σ=(K2O+Na2O)2/(SiO2-43)；利用 Al2O3、Na2O、TiO2計算戈蒂里指數(shù) τ =（Al2O3-Na2O）/TiO2；將計算結果都取常用對數(shù)，將取對數(shù)后的值投影到圖16中，圖中的數(shù)據是內蒙古西烏旗一帶的火山熔巖數(shù)據，它們投影點落在造山帶區(qū)域，與實際大地構造環(huán)境是吻合的。

圖16 里特曼-戈蒂里圖Fig.16 Rittmann- Cottine diagram

對其中的 A、B、C的分區(qū)說明等在幫助按鈕彈出的窗口中有說明，如其中的B區(qū)表造山帶，A區(qū)表非造山帶，C區(qū)表A、B區(qū)派生的堿性巖

4.2 F1-F2-F3 判別圖【4】

分析SiO2、TiO2、Al2O3、FeO、MgO、CaO、Na2O、K2O的含量，并利用它們計算 F1、F2、F3【4】（具體見邱家驤等《巖石化學》），并投影到圖17，此圖只適合于判斷玄武巖的構造環(huán)境。

4.3 板塊移動速度【4】

由于火山巖主要分布在板塊的擴張邊緣（如洋中脊）及閉合邊緣（如消減帶）；板塊擴張閉合速度與火山巖中K2O含量反相關，利用K2O含量可以確定閉合或者擴張速度(圖18)。

圖 17 F1-F2-F3與玄武巖關系圖Fig.17 The relation diagram of F1-F2-F3 and Basalt

圖 18 火山巖中K2O與板塊移動速度關系圖Fig.18 The relation diagram of K2O and Plate moving velocity in Volcanic rock

圖 18中的點是內蒙西烏旗的某流紋巖的數(shù)據，結合前面分析該區(qū)是造山帶環(huán)境，表明這里的封閉速度是2.44cm/a。

4.4 花崗巖構造環(huán)境分析R1-R2圖【4】

圖 19 中 R1=4Si-11(Na+K)-2(Fe+Ti)；R2=6Ca+2Mg+Al(為原子數(shù))，它們是判斷花崗巖構造環(huán)境的常用方法。對其中的參數(shù)和分區(qū)代號所顯示的構造環(huán)境，可按“幫助”在其信息窗中顯示。

圖19 花崗巖構造環(huán)境分析圖Fig.19 Tectonic environment analysis diagram for Granite

4.5 大洋與大陸巖漿巖判別圖解

這是趙崇賀1989年提出的一種判斷大洋、大陸裂谷和島弧造山帶巖漿巖的圖解，比較通用，既可用于侵入巖也可用于噴出巖，圖20的21、23號樣品分別為內蒙古西烏旗一帶的玄武巖和粗玄巖，由圖可知，玄武巖和粗玄巖樣產在島弧造山帶巖漿巖區(qū)中，符合該地的大地構造條件。

4.6 基中性火山巖構造環(huán)境分析圖【4】

顯然圖21只適合中性到基性的火山巖分析，分區(qū)標志等見“幫助”按鈕。圖中所標21、22 、23號樣為內蒙西烏旗一帶的玄武巖數(shù)據，3區(qū)是造山帶，與以其他參數(shù)確定的構造環(huán)境吻合。

4.7 花崗巖巖漿結晶與交代類型分析【4】

花崗巖有巖漿結晶者，也有交代（花崗巖化）者，如何區(qū)分可以采用圖22，對氧化物進行分析，得出其中的K、Na、Ca原子量相對百分含量投影到圖22中，如果點落在中間的深色區(qū)域就是巖漿結晶花崗巖，否則就是交代花崗巖。

圖20 大洋與大陸巖漿巖的判別圖解(趙崇賀，1989）Fig.20 Distinguished diagram of magma rock for ocean and continent

圖21 基中性火山巖構造環(huán)境分析圖Fig.21 Tectonic Environment Analysis diagram for the basite-mesite

圖22 巖漿花崗巖投影點分布區(qū)（據R.D.Raju等，1972）Fig.22 Projective dot distribution zone for Granite

巖漿結晶花崗巖可以繼續(xù)分類為S型和I型花崗巖，在深色區(qū)域，可以選擇右鍵菜單“I/S類型花崗巖分類”，彈出對話窗，見圖23，確定是I或者S型，也可以將類型定為IS混合型或者SI混合型，IS混合型說明是以I型為主的混合型，SI說明是以S型為主的混合型。

圖23 I與S型花崗巖在化學成分上的不同F(xiàn)ig.23 Chemical ingredient characters of the I and S type granite

在確定混合型時，以參數(shù)滿足S型化學成分特征多的就定義成SI型，否則就定義成IS型。有關I型和S型的定義點擊“幫助”能顯示出作為I與S類型確定的參考。

上述方法確定可以避免其他方法把花崗巖巖石截然分為I和S類型的局限。

4.8 I/S 花崗巖分類【4】

該分類是利用Na、K、Ca原子的含量特點來分類（圖 24），與前面的分類方法不一樣，前面的分類有點人為因素，這里直接分類出來，但這種分類沒法確定混合類型，過于死板。

圖24 花崗巖I/S自動分類圖Fig.24 Automatic classification diagram for granite

4.9 花崗巖的分類圖【4】

此種分類是從產出的構造環(huán)境上來分的，所以它從另一種途徑來分析巖石的形成環(huán)境（圖25）。

圖25 花崗巖分類圖Fig.25 Classification diagram for granite

4.10 玄武巖結晶溫度分析

4.10.1 計算法【4】在選擇數(shù)據后，在右鍵菜單“計算法”中彈出以下對話窗（圖 26），即可利用《巖漿巖巖石學》中提出的方法計算出玄武巖中橄欖石和斜長石的結晶溫度。

圖26 玄武巖結晶溫度對話窗Fig.26 Crystal temperature dialog for Basalt

4.10.2 圖解法【4】選擇數(shù)據后，在右鍵菜單“圖解法”中彈出對話窗（圖27），即可以估計礦物結晶溫度。對于圖中的分區(qū)及線的顏色等意義，可見“幫助”按鈕的彈出窗口。紅實線表橄欖石結晶溫度，藍虛線表斜長石結晶溫度，如圖27中的27號玄武巖樣，表明斜長石結晶溫度為1100℃，橄欖石結晶溫度為1070℃左右。說明橄欖石結晶不絕對會比斜長石要早，在特定的環(huán)境和成分比例情況下，斜長石可以比橄欖石結晶早。

圖27 玄武巖中MgO、Al2O3含量與橄欖石、斜長石結晶溫度及玄武巖類型關系圖Fig.27 Relation diagram among Olivine and plagioclase crystal temperature、 basaly type and MgO、Al2O3 Contents

5 巖漿堿性程度等分析

5.1 中酸性巖漿的堿度分析【4】

利用堿度率和 SiO2可以把巖漿巖劃分為鈣堿性、堿性和過堿性三級，它主要適合于中酸性巖石，過于酸性的巖石或過于基性的巖石不要使用此圖來判斷，要用其他方法，圖28為內蒙古西烏旗的火山巖數(shù)據，1為花英崗斑巖，9為石英閃長巖，說明是鈣堿性的，與本區(qū)大興安嶺發(fā)育的火山巖堿性特征吻合。

圖28 AR-SiO2與堿度關系圖Fig.28 Relation diagram of AR-SiO2 relation diagram

5.2 花崗巖堿度分析【4】

彈出圖29對話窗，分為堿性、亞堿性和鈣堿性三個類型，這一般只適合于中酸性巖石，其他巖石類型不合適。

圖29 花崗巖堿度分析圖Fig.29 Alkalinity analysis for granite

5.3 火山巖堿性和亞堿性系列的劃分

在選擇了火山巖數(shù)據后，上述菜單項有兩個子菜單，說明有兩種分類方法。

5.3.1 硅-堿圖法【4】彈出圖29對話窗即把巖石樣分成了堿性和亞堿性系列，在亞堿性系列（S區(qū)）又可右鍵彈出菜單“亞堿性進一步分類”，而彈出圖30的AFM對話窗，其中的A=K2O+Na2O，F(xiàn)=FeO+0.9Fe2O3，M=MgO，均百分含量。

圖30 硅-堿圖Fig.30 Silicon-Alkalinity diagram

圖31 A-F-M圖Fig.31 A-F-M diagram

利用圖31，可以將亞堿性系列劃分為拉斑玄武巖系列和鈣堿性系列。上述樣品是內蒙西烏旗數(shù)據，可見最終確定為鈣堿性系列，與通過其他方法確定的本區(qū)火山巖特征吻合。

5.3.2 Ol’-Ne’-Q’圖解法【5】圖 32 是比較可靠的堿性和亞堿性系列劃分方法，但一般只適合于中基性巖投影，類似地在亞堿性系列（S區(qū)）又可右鍵彈出菜單“亞堿性系列分類”，而彈出圖33對話窗，把亞堿性系列劃分為拉斑玄武巖系列和鈣堿性系列。

圖32 Ol’-Ne’-Q’圖Fig.32 Ol’-Ne’-Q’ diagram

圖33中的21和23號樣為某一玄武巖樣，說明為鈣堿性系列，符合其他方法確定的西烏旗一帶鈣堿性火山巖總特征。

圖33 Al2O3-“An”圖Fig.33 Al2O3-“An” diagram

6 巖漿成礦可能性分析

6.1 鐵銅礦床分析【4】

選擇巖漿巖進行鐵銅礦床分析，判斷是鉀質類型還是鈉質類型，這能判斷如果產礦的話將會是什么類型的礦（圖 34），因鈉質類型與鐵礦有關，而鉀質類型與銅礦有關，但并不能說明一定會有礦，必須依據其他條件一起分析。

圖34 巖漿巖與鐵銅礦床的關系Fig.34 Relation dialog for Magma rock and Iron-Copper mineral deposit

6.2 基性巖含礦性分析【4】

利用鈣、硅的原子百分數(shù)來確定是銅鎳礦還是鈦鐵礦，但只是可能，圖35中的23號樣為玄武安山巖樣，說明可能產鈦鐵礦，具體是否有礦要綜合其他資料分析。

圖 35 基性巖含礦可能性分析圖Fig.35 Mineral-bearing possibility analysis diagram for basite

6.3 金伯利巖的含礦性分析【5】

金伯利巖是種堿性的超基性巖，它與煌斑巖是兩種可能產金剛石的火成巖，據研究，大約只有10%的金伯利巖才產金剛石，池際尚教授提出T.A.值來判斷金伯利巖是否可能產金剛石，其中T.A=Al2O3+K2O+Na2O+TiO2+P2O5。這里都是百分含量，該值4.0～6.5者富礦，6.5～9.5者次富礦-貧礦，大于9.5者基本不含礦。計算以后以圖36所示的窗口表示如下：

圖36 金伯利巖含金剛石分析對話窗Fig.36 Diamond-bearing analysis dialog for Kimberlite

6.4 超基性-基性巖與礦產的關系【4】

這是利用吳利仁的研究，初步判斷超基性-基性巖的m/f與礦產的關系，圖37說明了可能的成礦類型，但并不是肯定有礦，但若有礦很可能是這種礦產。選擇“幫助”按鈕，可以得到m/f的指數(shù)說明窗口。

圖37 超基性-基性巖與礦產關系圖Fig.37 Relation dialog for the ultrabasite-basite and mineral type

6.5 鈦鐵銅鎳鉻的巖漿成礦分析【5】

經過國內已知的百余個鉻、銅鎳和鐵礦床的含礦巖相研究，發(fā)現(xiàn)m/f對M/S（含義見圖形說明）與其成礦類型有如下的分區(qū)規(guī)律，利用該圖可以分析可能的成礦類型，但具體是否有礦還得其他方法補充，圖38的20號樣表明可能產鈦磁鐵礦。

圖38 含礦巖體的m/f對M/S投影圖Fig.38 Projective diagram of m/f and M/S for mineral-bearing rock

7 結語

本軟件在開發(fā)過程中得到過中化地質礦山總局地質研究院李響在花崗巖形成構造環(huán)境分析方面的資料協(xié)助，另得到過欒俊霞、潘林林、劉勇等在圖版方面的協(xié)助，在這里對他們的幫助表示衷心感謝，由于巖石化學涉及內容極其廣泛，圖版更是多種多樣，這里不可能全部涉及到，有待以后來補充。該軟件是用 Qt平臺開發(fā)的【6】，所以既可用在Window操作系統(tǒng)中，也可用在Unix操作系統(tǒng)中，是跨平臺的。由于學識水平有限，可能還存在一些不足之處，請大家諒解。對于本軟件有需要者請聯(lián)系中化地質礦山總局地質研究院辦公室或者直接聯(lián)系軟件開發(fā)者本人丁衛(wèi)平（電郵見首頁）。

1 邱家驤. 巖漿巖巖石學[M]. 北京：地質出版社，1985

2 Hugh R.Rollison著. 巖石地球化學[M]. 楊學明、楊曉勇、陳雙喜，譯. 合肥：中國科學技術大學出版社，2000

3 中華人民共和國國標標準. 巖石分類和命名方案火成巖巖石分類和命名方案 GB/T 17412.1-1998[S]

4 邱家驤、林景仟. 巖石化學[M]，北京：地質出版社，1991

5 戚長謀，李鶴年，鄒祖榮. 地球化學通論[M]，北京：地質出版社，1984

6 Jasmin Blanchette，Mark Summerfield著..,C++ GUI Qt3 編程[M]. 齊亮，譯. 北京：北京航空航天大學出版社，2006