胡文祥 王 浩 馬 琦 陶愛(ài)華 劉西恩

(1 同濟(jì)大學(xué)聲學(xué)研究所 上海 200092)

(2 中海油田服務(wù)股份有限公司油田技術(shù)事業(yè)部 北京 101149)

0 引言

油氣井巖壁與井內(nèi)鋼套管之間的環(huán)形空間通常灌注不同類型水泥進(jìn)行井壁保護(hù)和不同深度巖石層之間的封隔。該結(jié)構(gòu)中，套管銹蝕等損傷以及管外水泥膠結(jié)質(zhì)量問(wèn)題會(huì)對(duì)油氣井的生產(chǎn)產(chǎn)生明顯影響。比如，環(huán)形空間由于水泥灌注缺陷會(huì)導(dǎo)致油氣井不同深度水層與油氣層之間的串漏，產(chǎn)生嚴(yán)重問(wèn)題。因此，進(jìn)行套管井檢測(cè)，以獲取套管損傷及水泥膠結(jié)缺陷信息是十分重要的工作。

聲波測(cè)井技術(shù)是目前進(jìn)行套管井檢測(cè)的最主要甚至唯一的手段。傳統(tǒng)的檢測(cè)方法是所謂水泥膠結(jié)測(cè)井(Cement bonding logging,CBL)。通過(guò)井軸上液體中20 kHz 左右低頻的對(duì)稱單極子聲源激發(fā)和不同深度位置接收獲得的聲信號(hào)，來(lái)了解和分析鋼套管與井壁巖石之間的水泥膠結(jié)信息。其中主要信息是由鋼套管管壁內(nèi)傳播模式(俗稱套管波)的幅度與衰減。該信息對(duì)鋼管外材料的特性阻抗是敏感的。因此，管外為良好的固體水泥粘結(jié)時(shí)幅度顯著減小，而水泥缺失為液體時(shí)，則有明顯的幅度，這實(shí)際上是管外材料的一種阻尼作用。通常認(rèn)為這種管壁內(nèi)的波模式為最低階的蘭姆波對(duì)稱模式S0。

傳統(tǒng)的CBL 測(cè)井通常用于定性判斷水泥膠結(jié)質(zhì)量。近20年來(lái)由于生產(chǎn)需求，進(jìn)行了一些研究來(lái)探討定量應(yīng)用的可能[1?2]。但很顯然，CBL 測(cè)井技術(shù)的確存在明顯的不足：如無(wú)法獲得水泥膠結(jié)缺陷的方位信息，無(wú)法有效進(jìn)行輕質(zhì)水泥情形評(píng)價(jià)，測(cè)量結(jié)果明顯受井內(nèi)測(cè)量聲探頭的偏心影響，快速地層時(shí)無(wú)法有效評(píng)價(jià)等。

超聲脈沖回波測(cè)井技術(shù)是20 世紀(jì)90年代后期以來(lái)開(kāi)始應(yīng)用的套管井檢測(cè)技術(shù)。該技術(shù)采用沿井周向旋轉(zhuǎn)掃描的寬帶自發(fā)自收超聲探頭垂直于套管壁發(fā)射和接收寬帶超聲脈沖。由于井液-鋼管顯著的阻抗失配，只能通過(guò)激發(fā)鋼管的厚度共振模式來(lái)實(shí)現(xiàn)對(duì)鋼管的透射。通過(guò)采用約200 kHz～700 kHz 范圍的寬帶超聲脈沖來(lái)覆蓋較大的套管厚度范圍。通過(guò)該檢測(cè)技術(shù)可以獲得管外介質(zhì)的聲阻抗圖像，判斷管外介質(zhì)性質(zhì)。該方法對(duì)管外水泥膠結(jié)缺陷有很好的周向分辨率。該檢測(cè)技術(shù)還可以較為有效地檢測(cè)套管的厚度，對(duì)套管壁的腐蝕與銹蝕進(jìn)行成像，是一種有效的套管井超聲成像技術(shù)。近年來(lái)的研究顯示，所激發(fā)的套管的這種波模式實(shí)際上與所謂S1零群速度蘭姆波模式(S1-ZGV共振)相關(guān)[3?5]。

但上述超聲脈沖回波測(cè)井技術(shù)亦有明顯不足，如對(duì)與鋼管接觸的管外材料的切向耦合不敏感。因此，對(duì)管外膠結(jié)良好的輕質(zhì)水泥與無(wú)水泥時(shí)的液體二者無(wú)法區(qū)分。此外，該技術(shù)亦無(wú)法檢測(cè)水泥-巖石界面(俗稱二界面)，當(dāng)然也無(wú)法判斷環(huán)形空間內(nèi)的材料性質(zhì)，以及井眼結(jié)構(gòu)等幾何信息。

近年來(lái)的研究顯示，針對(duì)上述檢測(cè)方法的不足，可以進(jìn)一步利用鋼套管蘭姆波的其他模式來(lái)實(shí)施套管井的超聲檢測(cè)。最低階反對(duì)稱A0模式是一種很好的可以補(bǔ)償上述不足的檢測(cè)手段[6?7]。由于管內(nèi)外介質(zhì)的低聲速，套管A0模傳播時(shí)會(huì)產(chǎn)生向管內(nèi)外介質(zhì)的聲能泄漏，所謂泄漏A0模式(或泄漏彎曲波模式)。一方面可利用泄漏A0模式對(duì)管外介質(zhì)的切向耦合敏感來(lái)區(qū)分輕質(zhì)水泥和液體；另一方面通過(guò)檢測(cè)管外泄漏A0模式的轉(zhuǎn)換波在井壁的反射來(lái)檢測(cè)和成像井孔幾何形狀以及套管在井中位置圖像。最重要的性質(zhì)是，由于A0模式傳播時(shí)的泄漏衰減與管外耦合材料特性緊密相關(guān)，因此通過(guò)其衰減的測(cè)量，可實(shí)現(xiàn)對(duì)管外材料性質(zhì)的有效定征。

上述不同檢測(cè)方法都利用了鋼套管內(nèi)的導(dǎo)波模式，它們用于鋼套管銹蝕及管外介質(zhì)的檢測(cè)與成像各有其優(yōu)缺點(diǎn)。因此，如果組合使用這幾種檢測(cè)方式，將可以有效利用它們所有的長(zhǎng)處，克服各自的不足。具體的工程實(shí)施方式，可以脈沖回波井周掃描系統(tǒng)為基礎(chǔ)，增加一組斜入射的泄漏A0模式收發(fā)探頭實(shí)現(xiàn)組合掃描成像；同時(shí)可以集成傳統(tǒng)的CBL測(cè)試方式，實(shí)現(xiàn)綜合檢測(cè)分析。本項(xiàng)研究是在多年前本課題組所開(kāi)發(fā)的套管井多功能超聲成像系統(tǒng)基礎(chǔ)上的一種全新改造與功能升級(jí)。

本文將主要闡述A0模式特點(diǎn)，其衰減的測(cè)量，結(jié)合脈沖回波測(cè)量進(jìn)行導(dǎo)波檢測(cè)與成像系統(tǒng)的開(kāi)發(fā)與綜合應(yīng)用。

1 負(fù)載鋼板A0 模式的基本特性

由于鋼套管內(nèi)徑較大，在較高的超聲頻段(比如大于100 kHz)時(shí)，沿管內(nèi)壁縱向斜入射的超聲脈沖，可近似看成與平面鋼板的作用。由此理論計(jì)算其相速度與衰減，相關(guān)理論分析與計(jì)算方法作者已在其他文章討論[7?9]，不在此一一列出。圖1 分別為板兩側(cè)均為液體時(shí)A0與S0模式的相速度與群速度頻散曲線。可以看到，較低頻厚積時(shí)，相對(duì)于S0模式，在一定頻厚積范圍，A0模式相速度與群速度的變化較小。對(duì)于通常的8～10 mm 厚套管，在100 kHz～400 kHz 頻率范圍，雖然其相速度有一定變化，但群速度較為平直，近似為“常群速度”。因此，對(duì)該模式波包可以實(shí)現(xiàn)較為精確的聲時(shí)與包絡(luò)幅度檢測(cè)，適合于實(shí)際工程超聲應(yīng)用。

磚子乍進(jìn)門，見(jiàn)纖掌揮來(lái)，明知敵不過(guò)刀馬旦的拳腳功夫，也要拼力反抗、擊退氣勢(shì)洶洶的強(qiáng)敵的，然而氣吐如蘭的動(dòng)口不動(dòng)手，將他軟化在石榴裙下了。趙仙童擁著他坐上沙發(fā)，轉(zhuǎn)身沏杯龍井茶，又削了一只大蘋果，剖成片放進(jìn)碟子里，捏出一片撒著嬌要喂他。磚子涌起潮水般的感動(dòng)，晚上遭受詩(shī)友們不知情的污辱，白天挨的莫名其妙耳光，全都在一片蘋果、一杯清茶中和解了。盡管這是他多年夫妻生活中，從未經(jīng)歷過(guò)的大起大落。

圖1 雙側(cè)水荷載時(shí)鋼板導(dǎo)波模式頻散特性Fig.1 The dispersion curves of guided waves for a steel plate loaded with water on both sides

對(duì)A0模式的激勵(lì)，按斯奈爾定律，需要考慮泥漿聲速的變化范圍以及A0模式的相速度等因素。對(duì)于1300～1700 m/s 的泥漿聲速范圍，可以估算超聲探頭的入射角范圍在約28?～38?。

考慮液體荷載下鋼板及鋼板下部的不同粘結(jié)材料，可進(jìn)一步計(jì)算A0模式的相速度頻散及其泄漏衰減。圖2 為板下部粘結(jié)層分別為空氣、水、慢速水泥、快速水泥四種條件下的結(jié)果。

可以看到，雖然管外不同粘結(jié)材料時(shí)A0模式的相速度頻散變化不大，但其衰減對(duì)管外材料非常敏感?？諝鈺r(shí)，僅單面液體泄漏，衰減最??；慢速水泥時(shí)，A0模式泄漏分別向水泥輻射縱波與橫波，因此，衰減最大；鋼管兩側(cè)均為水時(shí)，泄漏衰減介于二者之間。需要注意的是快速水泥情形，由于其縱波聲速大于A0模式相速度，故A0模式的泄漏僅向水泥輻射橫波，因此，其衰減反而小于慢速水泥情形，甚至小于兩側(cè)液體情形。圖3 為頻率200 kHz時(shí)，理論計(jì)算的A0模式泄漏衰減隨鋼板下側(cè)粘結(jié)材料特性阻抗的變化曲線。計(jì)算時(shí)考慮到鋼板下側(cè)粘結(jié)材料密度變化范圍相對(duì)較小，固定其為1800 kg/m3。在保持其縱橫波速度的比值為1.8708、即泊松比為0.3固定不變時(shí)，改變粘結(jié)材料的縱波速度從1000 m/s增加到4500 m/s，步長(zhǎng)10 m/s，數(shù)值計(jì)算A0模式泄漏衰減隨材料特性阻抗變化的連續(xù)曲線。圖3 中分別標(biāo)識(shí)了幾個(gè)典型材料的位置，旨在直觀說(shuō)明氣體、液體(水)及典型的水泥在特性阻抗-A0模衰減坐標(biāo)平面上是明顯分離的。因此通過(guò)分別測(cè)量套管井管外材料特性阻抗與沿管壁傳播的A0模衰減，可以有效分類管外材料。

圖2 水-10 mm 厚鋼板-粘結(jié)材料變化時(shí)A0 模式頻散與衰減Fig.2 Phase velocities of A0 mode and its attenuation for a steel plate loaded with water on one side and different bonded materials on other side

圖3 A0 模式衰減隨鋼板下側(cè)粘結(jié)材料特性阻抗的變化曲線Fig.3 The attenuation of A0 mode versus characteristic impedance of materials bonded on a steel plate

需要說(shuō)明的是，圖3 中存在一個(gè)所謂臨界特性阻抗點(diǎn)。該點(diǎn)前后A0模泄漏衰減隨鋼板粘結(jié)材料特性阻抗變化的轉(zhuǎn)折是容易解釋的，原因上文其實(shí)已經(jīng)說(shuō)明。該點(diǎn)處粘結(jié)材料縱波聲速與A0模相速度相等，為2740 m/s；而該點(diǎn)之前，A0模相速度大于粘結(jié)材料縱波聲速；該點(diǎn)之后，A0模相速度則小于粘結(jié)材料縱波聲速。按斯奈爾定律，則該點(diǎn)之前A0模會(huì)同時(shí)向粘結(jié)材料一側(cè)輻射縱波與橫波，而該點(diǎn)之后則僅會(huì)輻射橫波。因此該點(diǎn)之后，A0模的泄漏衰減明顯減小。

上述A0模式的衰減特性為管外介質(zhì)性質(zhì)的超聲定征提供了十分有用的思路。

2 不同負(fù)載鋼板脈沖回波與泄漏A0 模特性的實(shí)驗(yàn)測(cè)試

2.1 脈沖回波與泄漏A0 模信號(hào)的基本特征

僅以輕質(zhì)水泥為例說(shuō)明。用于測(cè)試的鋼板厚10 mm，寬200 mm，長(zhǎng)1000 mm。在其長(zhǎng)度的一半制作黏附厚度為30 mm 的水泥，長(zhǎng)度的另一半仍為自由的鋼板。水泥密度為1600 kg/m3，聲速為2100 m/s，為低密度、低聲速水泥。根據(jù)理論計(jì)算，此時(shí)粘結(jié)水泥段鋼板A0模衰減大于水。將探頭陣沿試樣中心線(紅色虛線)掃描，長(zhǎng)度為400 mm，如圖4 所示。其中T2 為脈沖回波探頭，T1 與R1、R2為斜探頭。

圖4 平面測(cè)試試樣與探頭陣示意圖Fig.4 Plane testing sample and ultrasonic transducer array

圖5 雙面水負(fù)荷鋼板與雙面分別為水與水泥負(fù)荷鋼板的脈沖回波及其幅度譜Fig.5 Pulse echoes and their amplitude spectrums for a steel plate loaded with water on both sides and water on one side and cement on other side

圖6 分別是鋼板雙面水負(fù)荷，雙面分別為水和水泥時(shí)A0模式測(cè)量信號(hào)。藍(lán)色是近探頭接收信號(hào)，紅色是遠(yuǎn)探頭接收信號(hào)。當(dāng)衰減大時(shí)(如鋼板背后是水泥)，接收幅度較小，而衰減小時(shí)(如鋼板背后是水)，接收幅度大。

圖6 雙面水負(fù)荷鋼板與雙面分別為水與水泥負(fù)荷鋼板的A0 模式測(cè)量信號(hào)Fig.6 Experimental signals of A0 mode for a steel plate loaded with water on both side and water on one side and cement on other side

2.2 鋼板兩側(cè)分別為水與空氣時(shí)的測(cè)試結(jié)果

測(cè)試水槽中，在鋼板背后的一段放置氣球(夾持使其貼合鋼板)，隔離水，模擬空氣情況。需要指出，該模擬條件可能導(dǎo)致所測(cè)量的空氣段衰減比實(shí)際偏大。圖7(a)為脈沖回波測(cè)試結(jié)果。圖7(a)中比值為共振波與表面反射波之比，厚度基于共振頻率計(jì)算。鋼板一側(cè)為空氣時(shí)，共振波幅度明顯增大，但厚度測(cè)量結(jié)果不受鋼板負(fù)荷變化影響。圖7(b)為A0模式泄漏衰減測(cè)量結(jié)果。鋼板雙面為水時(shí)，衰減1.16 dB/cm，背后空氣時(shí)衰減0.67 dB/cm。據(jù)此計(jì)算單面水引起的衰減約為0.49 dB/cm。雙面水引起的衰減為0.98 dB/cm，但實(shí)際上雙面水時(shí)衰減為1.16 dB/cm。其差別0.18 dB/cm 的原因，一方面如上所述，模擬方式導(dǎo)致單面空氣時(shí)的衰減測(cè)量偏大，另一方面與換能器靈敏度差異、幾何擴(kuò)散衰減、測(cè)量誤差等因素有關(guān)。

圖7 鋼板雙面分別為水和空氣時(shí)兩種導(dǎo)波方式測(cè)量結(jié)果Fig.7 Experimental measuring results of two kinds of guided wave for a steel plate loaded with water on both side and water on one side and air on other side

2.3 兩側(cè)分別為水負(fù)荷與輕質(zhì)水泥時(shí)的測(cè)試結(jié)果

圖4 的試樣，1 m 長(zhǎng)鋼板僅一半長(zhǎng)度覆蓋輕質(zhì)水泥。圖8為該試樣水中脈沖回波測(cè)試結(jié)果，同前述方法，計(jì)算歸一化共振波幅度及鋼板厚度。圖8(a)中雙側(cè)水負(fù)荷端(自由)歸一化共振波幅度為0.19，水泥部分歸一化幅度0.11，二者區(qū)分明顯，界面處為明顯的臺(tái)階；而所計(jì)算鋼板厚度曲線基本為10 mm的直線，不受鋼板負(fù)荷變化的影響。圖8(b)中將上述共振波幅度轉(zhuǎn)化為特性阻抗，水泥特性阻抗約3.5 MRayls，接近于上述實(shí)驗(yàn)測(cè)量的3.3 MRayls。

圖9 為兩個(gè)接收探頭的A0模式信號(hào)及包絡(luò)B掃描圖。由于水泥導(dǎo)致的衰減比水大，所以靠近水端的幅度較大，靠近水泥端幅度較小。

圖8 由圖4 試樣脈沖回波測(cè)試結(jié)果計(jì)算的鋼板厚度與負(fù)荷材料特性阻抗Fig.8 Thickness of steel plate and characteristic impedance of material bonded on steel plate calculating by using pulse echo measurement signals for sample in Fig.4

圖9 圖4 試樣A0 泄漏模式信號(hào)測(cè)試結(jié)果Fig.9 The experimental signals of leaky A0 mode for the sample in Fig.4

圖10 為上述測(cè)量方式對(duì)試樣掃描測(cè)量得到的衰減。很明顯，水泥段衰減大，而雙側(cè)水時(shí)衰減小。所得到的結(jié)果在水和水泥部分衰減曲線都相對(duì)平直，在分界面附近衰減呈線性變化，呈現(xiàn)斜的臺(tái)階，其值取決于兩個(gè)接收器所覆蓋的范圍內(nèi)水泥和水所占的比例。與前述的自發(fā)自收回波測(cè)量方式相似，該測(cè)量結(jié)果同樣很好反映了鋼板背面介質(zhì)的性質(zhì)。

顯然，上述脈沖回波測(cè)量方式在掃描方向上有更理想的分辨率，而A0測(cè)量方式則由于衰減計(jì)算涉及兩接收探頭之間的間距，測(cè)量分辨率降低。

在不考慮測(cè)量導(dǎo)致的因素前提下，兩側(cè)均為水造成的A0模衰減約為1.16 dB/cm，所以單側(cè)水的衰減約為0.58 dB/cm，由此可估算單側(cè)水泥衰減約為1.6?0.58=1.02 dB。

圖10 由圖9 信號(hào)計(jì)算的泄漏A0 模式的衰減Fig.10 The calculating result of attenuation for leaky A0 mode signals in Fig.9

組合兩種導(dǎo)波測(cè)試結(jié)果，進(jìn)行鋼板背面材料分類是可能的。考慮套管外環(huán)形空間環(huán)境條件下氣-液-固介質(zhì)特性阻抗-A0模衰減兩個(gè)參數(shù)可能的變化范圍，測(cè)量可能導(dǎo)致的離正確值的偏差，可在該兩項(xiàng)參數(shù)變化構(gòu)成的坐標(biāo)內(nèi)，給出氣-液-固每種狀態(tài)發(fā)生的可能分布范圍。將上述圖7、圖8與圖10的結(jié)果在該二維坐標(biāo)中標(biāo)識(shí)出，可以看到鋼板背面材料氣、液、固的區(qū)分是清楚的，見(jiàn)圖11。

圖11 鋼板背面氣-液-固粘結(jié)材料的區(qū)分Fig.11 The zoning map of gas-liquid-solid materials bonded on back surface of a steel plate

3 套管井導(dǎo)波成像系統(tǒng)開(kāi)發(fā)及模型井與實(shí)際井檢測(cè)與成像結(jié)果

在上述原理研究的基礎(chǔ)上，進(jìn)行了套管井導(dǎo)波成像系統(tǒng)研制開(kāi)發(fā)。整個(gè)系統(tǒng)由超聲發(fā)射、接收及信號(hào)采集、處理等電路裝置，檢測(cè)與成像軟件，垂直收發(fā)脈沖回波測(cè)量探頭，A0模激發(fā)與接收探頭(一發(fā)二收)，以及上位計(jì)算機(jī)組成。同時(shí)基于工程化考慮，超聲探頭均滿足耐溫耐壓和寬帶收發(fā)等特性，電路系統(tǒng)亦滿足耐溫耐壓要求。

圖12 模型井測(cè)量與成像結(jié)果Fig.12 The testing and imaging results in model testing well

圖13 實(shí)際井測(cè)量與成像結(jié)果Fig.13 The measuring and imaging results in a cased borehole in an oilfield

所研制的基于導(dǎo)波方式的套管井超聲檢測(cè)與成像儀器系統(tǒng)在中海油田服務(wù)股份有限公司多個(gè)模型井以及陸上與海上實(shí)際井進(jìn)行了測(cè)試，獲得良好效果。圖12 所給出的是其中一個(gè)典型模型井的成像結(jié)果。其套管外環(huán)形空間材料沿深度方向設(shè)置不同周向角度范圍缺失(填充近于液體的低特性阻抗仿體材料)通過(guò)聲阻抗率成像結(jié)果，A0模泄漏衰減結(jié)果圖像獲得很好體現(xiàn)。其中所得到的水泥聲阻抗約4 MRayls，泄漏A0模平均衰減約0.77 dB/cm，符合輕質(zhì)水泥較大衰減的特征。由二者綜合給出的管外介質(zhì)材料分類也得到了較為滿意的結(jié)果，表明了方法的有效性與系統(tǒng)開(kāi)發(fā)與工程應(yīng)用的實(shí)用性。

圖13 給出了實(shí)際井測(cè)量和成像結(jié)果。除了給出聲阻抗率圖像、彎曲波A0模衰減成像、環(huán)形空間介質(zhì)固-液-氣分類結(jié)果外，進(jìn)一步給出了據(jù)此得到的環(huán)形空間水力聯(lián)通圖像。該結(jié)果可以提供給工程人員形象的“竄槽”指示，有非常實(shí)際的工程價(jià)值。圖13 同時(shí)給出了傳統(tǒng)的CBL-VDL 與扇區(qū)水泥檢測(cè)儀器SBT結(jié)果，以作為對(duì)照。顯然，所開(kāi)發(fā)的導(dǎo)波成像結(jié)果是定量、直觀的，且有明確的周向位置，比傳統(tǒng)的CBL-VDL以及SBT結(jié)果明顯優(yōu)越。實(shí)際工程應(yīng)用時(shí)，可以綜合集成脈沖回波、A0彎曲波與傳統(tǒng)的CBL測(cè)量方式，進(jìn)行綜合成像分析。

利用上述的泄漏導(dǎo)波模式在環(huán)形空間巖石壁(二界面)的反射，通過(guò)周向掃描方式，還可進(jìn)一步進(jìn)行井眼形狀的成像、鋼套管在井內(nèi)相對(duì)位置的成像等，該結(jié)果不在此一一給出。

4 結(jié)論

研究表明，傳統(tǒng)CBL與超聲脈沖回波測(cè)井技術(shù)都可歸于鋼套管壁內(nèi)導(dǎo)波模式的利用。為了解決現(xiàn)有技術(shù)無(wú)法有效檢測(cè)輕質(zhì)水泥條件下的套管井水泥膠結(jié)質(zhì)量問(wèn)題，同時(shí)進(jìn)一步對(duì)管外環(huán)形空間內(nèi)的材料性質(zhì)進(jìn)行定征，成像井眼形狀及其中套管的相對(duì)位置，可基于已有的多功能超聲成像掃描成像測(cè)井儀，進(jìn)一步增加A0彎曲波模式的激勵(lì)和檢測(cè)方式。綜合的多模式導(dǎo)波檢測(cè)與成像方法，可以較為全面解決目前套管井檢測(cè)與評(píng)價(jià)所遇到的問(wèn)題。在泄漏A0模式性質(zhì)分析和實(shí)驗(yàn)?zāi)Ｐ蜏y(cè)試研究的基礎(chǔ)上，開(kāi)發(fā)相應(yīng)的超聲檢測(cè)換能器、電路與軟件系統(tǒng)，實(shí)現(xiàn)了套管井綜合的多模式超聲導(dǎo)波成像系統(tǒng)。通過(guò)中海油田服務(wù)股份有限公司在其模型井與陸地和海上油氣井的實(shí)際測(cè)試使用，均獲得良好效果，說(shuō)明研發(fā)工作是成功的。

致謝參加本項(xiàng)研發(fā)工作的還有同濟(jì)大學(xué)梁軍汀、盧杰、周八妹和中海油服王文梁等同仁，在此表示感謝。