劉 剛， 李祎宸， 張家林， 劉 闖， 楊 帆， 穆文軍， 王 鍇

(1. 中國石油大學(xué)(華東) 石油工程學(xué)院，山東 青島 266580；2. 濰柴動力股份有限公司 發(fā)動機技術(shù)研究院，山東 濰坊 261000； 3.中國電波傳播研究所，河南 新鄉(xiāng) 453000)

大型叢式井、大位移水平井以及加密調(diào)整井等技術(shù)的廣泛應(yīng)用推進(jìn)了海上油氣資源的高效開發(fā)[1]，同時也使得井間距離不斷減小，增加了井眼碰撞的風(fēng)險。為了有效監(jiān)測與控制井眼軌跡，保持良好的井身質(zhì)量，需要時刻關(guān)注鉆頭與風(fēng)險鄰井套管間的距離，并通過有效的防碰措施來保證安全的井間距離[2]。鉆井施工中，目前主要通過防碰掃描計算來進(jìn)行井間距離的評估，包括法面距離掃描、平面距離掃描和最近距離掃描三種方式[3]，但防碰掃描存在時間和井深的“滯后性”，且計算精度受諸多因素的干擾[4]。本課題組基于監(jiān)測鉆頭破巖時誘發(fā)的振動波完成了渤海與南海等多個區(qū)塊的防碰預(yù)警，并且結(jié)合彈性波的球面擴散和衰減吸收理論，建立了鉆頭趨近鄰井的距離計算模型，由于傳感器布置在套管頂端，使得信號在長距離傳輸過程中產(chǎn)生衰減，從而限制了距離模型的計算精度[5-6]。

本文以地球物理勘探中的“井間地震”為啟發(fā)[7]，通過在風(fēng)險井段內(nèi)的環(huán)空套管壁處布置四個傳感器接收振動信號能量，分析鉆頭震源與四個傳感器在空間所構(gòu)成的幾何關(guān)系，推導(dǎo)得出了基于多傳感器計算鉆頭與套管間距離的理論模型。根據(jù)振動波在地層中的傳播特征，設(shè)計了室外模擬多傳感器接收振動信號能量的試驗，在連續(xù)模擬套管的接箍處粘貼四個加速度傳感器形成測點，并將模擬套管與傳感器掩埋于地表淺層，在水泥地面合理劃分激勵點，利用鉛球于高處落地敲擊激勵點產(chǎn)生的振動波作為信號震源，對所建模型進(jìn)行了距離測試，采用遺傳算法對信號數(shù)據(jù)進(jìn)行處理，增加模型的容錯率。結(jié)果表明，該模型能夠大大縮短振動信號的傳播距離，提高采集數(shù)據(jù)的頻率和信噪比，有效預(yù)測鉆頭與鄰井套管的距離，為井間距離的預(yù)測提供新思路。

1 鉆頭振動波的理論概述

鉆頭在井底鉆進(jìn)工作時，會對地層產(chǎn)生一個較強的沖擊作用，進(jìn)而激發(fā)振動波產(chǎn)生[8-9]。與隨鉆地震相似，以鉆頭在鉆進(jìn)過程中破巖誘發(fā)的振動波作為井下信號源，對相應(yīng)的振動波場特征進(jìn)行分析研究。

1.1 鉆頭振動波的產(chǎn)生形式

假設(shè)鉆頭所在地層為無限大均勻各向同性介質(zhì)，鉆頭在其中相當(dāng)于點震源，根據(jù)鉆頭的破巖機理，鉆頭的輻射模式包括橫向剪切作用力引發(fā)的振動輻射，以及軸向沖擊作用力引發(fā)振動輻射，如圖1所示。

P-縱波；SH-橫波垂直平面內(nèi)分量；SV-橫波水平面內(nèi)分量；φ-鉆頭撞擊力與井軸夾角圖1 鉆頭輻射模式圖Fig.1 Bit radiation pattern graph

鉆頭旋轉(zhuǎn)鉆進(jìn)會產(chǎn)生強烈的軸向沖擊作用力，引發(fā)P波和SV波的產(chǎn)生，Haldorsen等[10]求得該點震源沿軸向振動的波場解析式

(1)

鉆頭在井底鉆進(jìn)時，其產(chǎn)生的橫向作用力相當(dāng)于一系列的力偶作用，會產(chǎn)生SH波，Haldorsen等求得其質(zhì)點位移表達(dá)式為

(2)

式中：ρ為地層密度；r為震源到波前面的距離；α，β為縱、橫波的速度；A1，A2為比例因子；φ為鉆頭撞擊力與軸向的夾角；d為鉆頭直徑；θ為xy平面測量的方位角；ur為P波輻射的徑向分量；uθ為SH波運動分量；uφ為SV波運動分量。

經(jīng)過上述分析可知P波只沿鉆頭軸向撞擊力方向輻射，其他方向則沒有；SV波則只沿與撞擊力垂直的方向輻射。鉆頭鉆遇地層的泊松比決定P波和SV波的相對輻射振幅，因為橫向作用力比垂向作用力小得多，所以SH波對質(zhì)點位移作出的貢獻(xiàn)最小。對于典型的振動頻率和鉆頭尺寸，SV波的振幅比SH波的振幅大約大三個數(shù)量級，在計算時可忽略SH波[11-12]。

1.2 鉆頭振動波的傳播特性

鉆頭鉆遇地層引起的振動信號屬于機械彈性波，波源的能量向周圍介質(zhì)傳播的過程其實就是彈性波的傳播過程。波動從振源向遠(yuǎn)處傳播時，介質(zhì)從波前到達(dá)時開始振動，能量也隨之逐層向外擴散傳播，地面?zhèn)鞲衅魉邮盏男盘柧褪嵌鄠€應(yīng)力波疊加的結(jié)果，因此可利用球形傳播模型，描述鉆頭振動波的產(chǎn)生與傳播過程[13]。

彈性波在連續(xù)介質(zhì)中逐層傳播時，由于實際地層巖石是非均質(zhì)和非完全彈性的，地層對彈性波會產(chǎn)生吸收作用，使得彈性波的能量逐漸減降低。通過求解黏彈性介質(zhì)的波動方程，可以從理論上證明彈性波在地層中傳播時其能量隨傳播的距離成指數(shù)形式衰減，其表達(dá)式為

Ar=A0e-x(f)r

(3)

式中：A0為震源的初始振幅；Ar為距離震源r處接收的信號振幅；x(f)為衰減系數(shù)；該系數(shù)的影響因素有兩方面，一方面為振動波傳經(jīng)地層的孔隙度，密度，潤濕程度等理化性質(zhì)；另一方面為振動波自身頻率，高頻信號衰減系數(shù)較大，低頻信號衰減系數(shù)較小。衰減系數(shù)x(f)反映了地層對振動波起到綜合濾波作用[14-15]。

2 井間距離理論模型的建立

2.1 多傳感器監(jiān)測模型的建立

對于海洋加密調(diào)整井，具有防碰監(jiān)測需求的多是淺層直井段，該段地層比較疏松[16]，且鄰井之間的距離較小，所以振動波在該段的傳播過中主要是地層吸收衰減。于防碰監(jiān)測井段的井筒環(huán)空套管內(nèi)壁處布置四個傳感器，同步采集鉆頭鉆進(jìn)地層時的振動信號能量，建立基于多傳感器監(jiān)測鉆頭與套管間距離的理論模型，如圖2所示。由于該井段所布置的傳感器距離較近，可以假設(shè)各個傳感器是處于同一個地層，因此衰減系數(shù)x(f)相近。

(a) 井筒環(huán)空套管內(nèi)壁處的多傳感器位置

(b) 鉆井作業(yè)時的井下鉆頭震源圖2 多傳感器監(jiān)測模型示意圖Fig.2 Schematic diagram of multi-sensor monitoring model

S為鉆頭震源，A，B，C，D是布置在井筒環(huán)空套管內(nèi)壁處的四個接收傳感器，為方便計算求解，將傳感器等間距放置，設(shè)AB=BC=CD=a(m)，鉆頭震源到風(fēng)險井套管的垂直距離為H，鉆頭震源到各個傳感器的直線距離分別為SA=rA，SB=rB，SC=rC，SD=rD。

2.2 井間距離計算模型的建立

設(shè)A，B，C，D四個傳感器接收的能量為EA，EB，EC，ED。由彈性波球形傳播理論可知，彈性波能量E正比于波振幅A的平方

E∞A2

(4)

則結(jié)合式(3)可得

(5)

在ΔSAC中，設(shè)參數(shù)為

(6)

對式(6)兩邊同時取對數(shù)可以得到

(7)

將rA，rC用rB來表示可得

(8)

設(shè)∠SAB=α，在ΔSAB和ΔSBC中，根據(jù)余弦定理，可得

(9)

因為cos(180° -α)=-cosα，所以將式(9)上下相加可得

(10)

將式(8)代入式(10)可得

(11)

根據(jù)式(11)可以得到距離rB

(12)

同理，在ΔSBD中，設(shè)參數(shù)為

(13)

與ΔSAC中的推導(dǎo)相似，可在ΔSBD中得到距離rC

(14)

將rB，rC代入式(8)可得

(15)

令

(16)

(17)

則有

Mx2(f)=N

(18)

可以求得衰減系數(shù)x(f)

(19)

根據(jù)傳感器測得的振動幅值A(chǔ)，可由式(6)、式(13)求得各個λ值，式(19)中，衰減系數(shù)x(f)只與λ1，λ2，λ3，λ4以及傳感器間距a有關(guān)，所以可計算得到衰減系數(shù)x(f)值。將x(f)代入式(12)、式(14)可求得rB，rC值。

對于任意三角形，根據(jù)海倫公式[17]，其面積S為

(20)

式中：b，c，d為三角形邊長；l為三角形的半周長。

故在ΔSBC中有

(21)

所以鉆頭震源到風(fēng)險井套管的垂直距離H為

(22)

3 井間距離理論模型的驗證

為檢驗所建立的基于多傳感器計算鉆頭與套管間距離理論模型的合理性，并探究傳感器組相對鉆頭的位置關(guān)系，設(shè)計了室外模擬多傳感器接收振動信號能量的試驗，對試驗所采集的信號數(shù)據(jù)采用遺傳算法進(jìn)行處理，增加模型的容錯率。

3.1 多傳感器信號采集的試驗

室外模擬多傳感器接收振動信號能量的試驗，如圖3所示。選擇周圍機器噪聲較小的露天寬闊水泥面處進(jìn)行試驗，在連續(xù)模擬套管的接箍處采用高強度化學(xué)膠結(jié)劑粘貼四個加速度傳感器形成測點，用來模擬井軸線方向，并將模擬套管與傳感器掩埋于水泥面邊緣處的地表淺層，初步滿足傳感器組近似在同一地層。合理劃分多個激勵點，用來模擬鉆頭位置，并標(biāo)記于水

圖3 多傳感器信號采集試驗Fig.3 Simulations to acquire multi-sensor signals

泥地面，將鉛球于高處落下敲擊地面激勵點產(chǎn)生的振動波作為信號震源，鉛球質(zhì)量為6 kg，每次下落高度均為1 m，以滿足同一個測試點每次產(chǎn)生相近的撞擊能及相近的機械波振動頻率。試驗儀器由加速度振動傳感器、微型電荷放大器、數(shù)據(jù)采集儀、采集工作站組成，加速度傳感器靈敏度均為30 pc/g，電荷放大器的放大功率為1 mv/pc，采樣頻率設(shè)置為102.4 kHz。

為探究當(dāng)鉆頭位置距離傳感器組所確定軸線方向較近(<5 m)和較遠(yuǎn)(<20 m)時的模型效果，分別劃分了激勵點1～15和激勵點16～20，用來模擬鉆頭的不同位置情況，如圖4和圖5所示。激勵點(模擬鉆頭位置)和傳感器組的相對位置關(guān)系，主要可分為兩大類：①激勵點位于傳感器組的正對方；②激勵點位于傳感器組的斜側(cè)方。為保證測試的準(zhǔn)確性與穩(wěn)定性，對每一個激勵點進(jìn)行多次重復(fù)敲擊。

圖4 近距離測試下激勵點分布Fig.4 Impacted points distribution under close distance

圖5 遠(yuǎn)距離測試下激勵點分布Fig.5 Impacted points distribution under long distance

3.2 基于遺傳算法的振動信號處理

Holland 1975年創(chuàng)造性地提出遺傳算法(Genetic Algorithm，GA)[18]，該算法模擬了自然選擇和遺傳中發(fā)生的復(fù)制、交叉和變異等現(xiàn)象，其本質(zhì)是一種高效、并行、全局搜索的方法[19], 近年來以其簡單通用 、穩(wěn)健性強、適于并行處理及高效實用等顯著特點, 在各個領(lǐng)域得到了廣泛應(yīng)用[20-21]。本文建立的鉆頭震源到風(fēng)險井套管的垂直距離計算公式中，認(rèn)為井筒環(huán)空套管內(nèi)壁處的傳感器所測振動是由鉆頭震源引起，但實際作業(yè)環(huán)境下，由于海浪、鉆臺振動等因素的影響，使得所測振動幅值包含噪音干擾項，根據(jù)此特點，為減少采集數(shù)據(jù)偏差對計算結(jié)果的影響，本文采用GA對信號數(shù)據(jù)進(jìn)行處理，數(shù)據(jù)處理的具體流程，如圖6所示。

圖6 基于遺傳算法的數(shù)據(jù)處理流程Fig.6 The data processing flow based on GA

選擇激勵點5的測試數(shù)據(jù)為例進(jìn)行處理其時域圖，如圖7所示。該激勵點模擬鉆頭位置處于傳感器組的斜側(cè)方，由于該位置處振動波傳至各傳感器的路徑相近，所以其等效衰減系數(shù)比位于傳感器組正對方的激勵點更加一致。

圖7 激勵點5時域圖Fig.7 The time domain of the 5th impacted point

應(yīng)用Matlab編程對激勵點5的測試數(shù)據(jù)進(jìn)行處理，通過峰值檢測分別獲得五次重復(fù)撞擊試驗下傳感器接收能量的峰值，如表1所示。為預(yù)判各數(shù)據(jù)的有效性，將五次重復(fù)撞擊試驗的峰值作圖分析，如圖8所示。

表1 傳感器接收能量的峰值統(tǒng)計

圖8 能量峰值擬合和回歸Fig.8 The fitting and regress of energy peak value

但傳感器C接收到的振動波幅值嚴(yán)重偏離其它三個傳感器接收的振動波幅值所確定的趨勢線，為了減少傳感器C采集的數(shù)據(jù)對距離預(yù)測的影響，采用閾值修正法進(jìn)行處理[22]，將傳感器C接收的振動幅值修正為趨勢線預(yù)測值與實際值和的一半，具體算法如下

當(dāng)(Ap-As)/Ap≤M

(23)

則Ad=0.5(Ap+As)

(24)

式中：Ap為四個傳感器振動幅值所確定的趨勢線預(yù)測值；As為要修正傳感器實際接收到的幅值；Ad為修正后的振動幅值；M為衡量偏離程度所設(shè)定的閾值。本次將閾值M設(shè)定20%，如果偏離程度小于閾值，則對傳感器采集到的振動幅值進(jìn)行修正，若是偏差大于該閾值，則使用實際值。

為優(yōu)化井間距離模型的計算結(jié)果，增加模型的容錯率，本文采用的GA適應(yīng)度值為輸入距離rA，rB，rC，rD與模型預(yù)測值rAp，rBp，rCp，rdp之間誤差的和e，當(dāng)誤差之和e滿足輸出要求時，GA優(yōu)化結(jié)束，輸出最終模型預(yù)測結(jié)果。本次GA的參數(shù)設(shè)定為[23]：個體數(shù)目NIND=40， 最大遺傳代數(shù)MAXGEN=500，變量數(shù)NVAR=4，變量的二進(jìn)制位數(shù)PRECI=3，代溝GGAP=0.9，其他參數(shù)設(shè)定為系統(tǒng)默認(rèn)值。結(jié)合基于GA的數(shù)據(jù)處理流程，優(yōu)化得到激勵點5下模型預(yù)測距離在進(jìn)化到200代，300代，400代，500代的誤差情況，最終優(yōu)化結(jié)果如圖9所示。

圖9 激勵點5遺傳算法優(yōu)化結(jié)果Fig.9 GA optimization result of the 5th impacted point

由圖9可知，個體數(shù)為40的種群，其誤差分布主要從第200代的20附近最后到達(dá)第500代的3附近，同時也有九個個體的誤差在優(yōu)化過程中逐漸趨于0，大部分個體誤差在5以下，GA的尋優(yōu)效果較好。

經(jīng)過GA優(yōu)化的數(shù)據(jù)，利用井間距離模型重新進(jìn)行計算，求解得到激勵點5與各個傳感器之間的距離半徑，并求得與傳感器組所在井軸線的垂直距離。將模型計算距離結(jié)果與實際值進(jìn)行對比，結(jié)果顯示，距離半徑計算值和真實值具有一定誤差，但多數(shù)維持在10%以內(nèi)，達(dá)到了預(yù)期效果，最后計算的激勵點與井軸距離的誤差為30%。具體結(jié)果見表2。

表2 撞擊點5距離半徑計算結(jié)果

激勵點5的分析中，已證實式(3)的有效性，所以本文認(rèn)為，對于四個傳感器接收到的能量數(shù)據(jù)，如果符合隨距離增大而降低則為有效數(shù)據(jù)，否則，判定數(shù)據(jù)無效。根據(jù)激勵點和傳感器組的相對位置關(guān)系，本文對位于傳感器組正對方的激勵點和位于傳感器組斜側(cè)方的激勵點分別進(jìn)行平均峰值統(tǒng)計。六個正對方激勵點與九個斜側(cè)方激勵點的統(tǒng)計結(jié)果，如表3所示。

分析表中數(shù)據(jù)可知，正對方的六個激勵點中，有三個激勵點數(shù)據(jù)滿足本文假設(shè)條件，有效數(shù)據(jù)獲得率為50%，獲取率較低。斜側(cè)方的九個激勵點均符合模型假設(shè)條件，有效數(shù)據(jù)獲取率為100%。對比結(jié)果可發(fā)現(xiàn)，傳感器組和激勵點的相對位置關(guān)系嚴(yán)重影響了有效數(shù)據(jù)的數(shù)量，激勵點位于傳感器組的正對方，增大了振動波傳播至各傳感器路徑的差異性，導(dǎo)致各條路徑的等效衰減系數(shù)差別較大，嚴(yán)重影響了模型假設(shè)條件。

對于有效的激勵點數(shù)據(jù)，采用激勵點5數(shù)據(jù)的處理方法和流程計算，求解得到各個激勵點與傳感器組所在井軸線的垂直距離，如表4所示。計算誤差均在30%以內(nèi)，多數(shù)在20%以內(nèi)，模型計算效果較好。但誤差范圍為4.7%～30%，分析后可知，實際的傳播地層理化性質(zhì)不一致，增大了振動波傳播至各傳感器路徑的差異性，導(dǎo)致了公式計算過程中，各個采集傳感器相應(yīng)衰減系數(shù)x(f)相近的假設(shè)不夠準(zhǔn)確，因此在實際的作業(yè)監(jiān)測中，應(yīng)該根據(jù)井深調(diào)整傳感器位置，保證傳感器能接收到鉆頭振動信號，也要盡可能使得鉆頭與采集傳感器之間的巖性、物性基本一致。對于衰減系數(shù)x(f)的計算，要根據(jù)實測數(shù)據(jù)結(jié)合防碰掃描數(shù)據(jù)進(jìn)行擬合驗證。

表3 激勵點數(shù)據(jù)有效性統(tǒng)計

表4 激勵點距離計算結(jié)果

對于遠(yuǎn)距離測試點數(shù)據(jù)，發(fā)現(xiàn)激勵點19的信號脈沖已經(jīng)淹沒，提取不到有效幅值，所以在試驗激勵條件下，傳感器組可以監(jiān)測到的最遠(yuǎn)距離為激勵點18所確定的12 m。同樣采用上述步驟分析，遠(yuǎn)距離測試下激勵點有效性分析結(jié)果與最終距離計算結(jié)果，見表5和表6。

表5 遠(yuǎn)距離點數(shù)據(jù)有效性統(tǒng)計

表6 遠(yuǎn)距離點距離計算結(jié)果

4 結(jié) 論

(1)在風(fēng)險井段內(nèi)的環(huán)空套管壁處布置四個傳感器接收振動信號能量，分析鉆頭震源與四個傳感器在空間所構(gòu)成的幾何關(guān)系，并結(jié)合振動波在地層中的傳播特征，可以推導(dǎo)得出基于多傳感器計算鉆頭與套管間距離的理論模型。

(2)傳感器組和鉆頭的相對位置關(guān)系嚴(yán)重影響有效數(shù)據(jù)獲取的成功率，鉆頭位于傳感器組的斜側(cè)方時，震源產(chǎn)生的振動波傳播到各傳感器的路徑相似，途徑地層的理化性質(zhì)一致強，使得其等效衰減系數(shù)較一致，有助于有效數(shù)據(jù)的獲取。

(3)利用GA對傳感器采集數(shù)據(jù)偏差進(jìn)行優(yōu)化，能夠增加鉆頭與套管間距離預(yù)測模型的容錯率。室內(nèi)試驗數(shù)據(jù)的計算結(jié)果，驗證了本文所建立模型的有效性。

參 考 文 獻(xiàn)

[ 1 ] 姜偉. 海上密集叢式井組再加密調(diào)整井網(wǎng)鉆井技術(shù)探索與實踐[J]. 天然氣工業(yè)，2011, 31(1): 69-72.

JIANG Wei. Research and practices of re-infill drilling technology for offshore close-spaced cluster-well pads [J]. Natural Gas Industry, 2011, 31(1): 69-72.

[ 2 ] 劉永旺，管志川，史玉才，等. 井眼防碰技術(shù)存在的問題及主動防碰方法探討[J]. 石油鉆采工藝，2011, 33(6): 14-18.

LIU Yongwang, GUAN Zhichuan, SHI Yucai, et al. Discussion on problems of wellbore anti-collision technology and active anti-collision methods [J]. Oil Drilling & Production Technology, 2011, 33(6): 14-18.

[ 3 ] 王飛躍，劉峰剛，薛燕，等. 防碰掃描方法在叢式鉆井中的應(yīng)用[J]. 遼寧化工，2010, 39(9): 925-928.

WANG Feiyue, LIU Fenggang, XUE Yan, et al. Application of anti-collision scanning method in cluster well drilling [J]. Liaoning Chemical Industry, 2010, 39(9): 925-928.

[ 4 ] POEDJONO B, ISEVCAN E, LOMBARDO G J, et al. Anti-collision and risk management offshore qatar: a successful collaboration ：IPTC 13142[R]. 2009-12-07.

[ 5 ] 劉剛，劉華亮，何保生，等. 叢式井鉆頭距鄰井井筒距離的計算模型推導(dǎo)[J]. 石油機械，2013, 41(10): 46-50.

LIU Gang, LIU Hualiang, HE Baosheng, et al. Calculation model of cluster well drilling bit-adjacent well distance [J]. China Petroleum Machinery, 2013, 41(10): 46-50.

[ 6 ] 劉剛，劉闖，夏向陽，等. 基于 PCA-SVDD 方法的鉆頭異常鉆進(jìn)識別[J]. 振動與沖擊，2015, 34(13): 158-162.

LIU Gang, LIU Chuang, XIA Xiangyang, et al. Drill bit abnormal drilling condition recognition based on PCA-SVDD [J]. Journal of Vibration and Shock, 2015, 34(13): 158-162.

[ 7 ] 段寶平，何正勤，葉太蘭. 井間地震技術(shù)的研究現(xiàn)狀及其工程應(yīng)用前景[J]. 物探與化探，2010, 34(5): 610-616.

DUAN Baoping, HE Zhengqin, YE Tailan, et al. Research progress and engineering application prospect of crosswell seismic technology [J]. Geophysical and Geochemical Exploration, 2010, 34(5): 610-616.

[ 8 ] 劉剛，劉闖，夏向陽，等. 基于HHT方法的鉆頭振動信號識別[J]. 振動與沖擊，2015, 34(12): 212-218.

LIU Gang, LIU Chuang, XIA Xiangyang, et al. Identification of vibration singals of drill bit based on Hilbert-Huang transform method [J]. Journal of Vibration and Shock, 2015, 34(12): 212-218.

[ 9 ] 劉剛，張家林，劉闖，等.鉆頭鉆進(jìn)不同介質(zhì)時的振動信號特征識別研究[J]. 振動與沖擊，2017, 36(8): 71-78.

LIU Gang, ZHANG Jialin, LIU Chuang, et al. An identification method of vibration signal features when bit drills different mediums [J]. Journal of Vibration and Shock, 2017, 36(8): 71-78.

[10] DENNEY D. Optimizing well construction: full-waveform data from while-drilling seismic measurements[J]. Journal of Petroleum Technology, 2003, 55(8): 36-37.

[11] 何保生，劉華亮，劉剛，等. 基于球面波傳播理論的鉆頭趨近鄰井防碰預(yù)測模型[J]. 石油鉆探技術(shù)，2013, 41(3): 62-66.

HE Baosheng, LIU Hualiang, LIU Gang, et al. Anti-collision prediction model for drill bit approaching adjacent well based on spherical wave propagation theory [J]. Petroleum Drilling Techniques, 2013, 41(3): 62-66.

[12] 黃光南，鄧居智，李紅星，等. P-SV波和P-SH波非線性旅行時層析成像[J]. 石油地球物理勘探，2015, 50(6): 1127-1133.

HUANG Guangnan, DENG Juzhi, LI Hongxing, et al. Nonlinear traveltime tomography of P-SV and P-SH waves [J]. Oil Geophysical Prospecting, 2015, 50(6): 1127-

1133.

[13] 巴振寧，梁建文，張艷菊. 三維層狀黏彈性半空間中球面SH、P和SV波源自由場[J]. 地球物理學(xué)報，2016, 59(2): 606-623.

BA Zhenning, LIANG Jianwen, ZHANG Yanju. Free-field responses of spherical SH-, P- and SV-wave sources in a layered visco-elastic half space [J]. Chinese Journal of Geophysics, 2016, 59(2): 606-623.

[14] 周光泉，劉孝敏. 黏彈性理論[M]. 合肥：中國科學(xué)技術(shù)大學(xué)出版社，1996.

[15] 祝效華，湯歷平，童華. 高頻扭轉(zhuǎn)沖擊鉆進(jìn)的減振與提速機理研究[J]. 振動與沖擊，2012, 31(20): 75-78.

ZHU Xiaohua, TANG Liping, TONG Hua. Rock breaking mechanism of a high frequency torsional impact drilling [J]. Journal of Vibration and Shock, 2012, 31(20): 75-78.

[16] 李紅星. 番禺30-1氣田叢式井淺層鉆進(jìn)防碰繞障技術(shù)[J]. 石油鉆探技術(shù)，2015, 43(6): 125-129.

LI Hongxing. Anti-collision and obstacle bypassing techniques in cluster wells drilling in shallow layers of the PY30-1 gas field[J]. Petroleum Drilling Techniques, 2015, 43(6): 125-129.

[17] 翟春杰，徐建閩，劉永桂. 基于分區(qū)的能耗均衡路由協(xié)議[J]. 傳感技術(shù)學(xué)報，2016, 29(1): 80-87.

ZHAI Chunjie, XU Jianmin, LIU Yonggui. Energy-consumption balancing routing protocol based on regions[J]. Chinese Journal of Sensors and Actuators, 2016, 29(1):

80-87.

[18] ESAT I I, BAHAI H. Vibratory system synthesis for multi-body systems based on genetic algorithm [J]. Journal of Sound & Vibration, 2000, 230(4): 933-948.

[19] 邊霞，米良. 遺傳算法理論及其應(yīng)用研究進(jìn)展[J]. 計算機應(yīng)用研究，2010, 27(7): 2425-2429.

BIAN Xia , MI Liang. Development on genetic algorithm theory and its applications [J]. Application Research of Computers, 2010, 27(7): 2425-2429.

[20] 付江華，史文庫，沈志宏，等. 基于遺傳算法的汽車動力總成懸置系統(tǒng)優(yōu)化研究[J]. 振動與沖擊，2010, 29(10): 187-190.

FU Jianghua, SHI Wenku, SHEN Zhihong, et al. Factors influencing shock characteristics of a cubic nonlinear packaging system with critical component [J]. Journal of Vibration and Shock, 2010, 29(10): 187-190.

[21] 王波，劉樹林，蔣超，等. 基于量子遺傳算法優(yōu)化RVM的滾動軸承智能故障診斷[J]. 振動與沖擊，2015, 34(17): 207-212.

WANG Bo, LIU Shulin, JIANG Chao, et al. Rolling bearings intelligent fault diagnosis based on RVM optimized with quantum genetic algorithm [J]. Journal of Vibration and Shock, 2015, 34(17): 207-212.

[22] 盧娜，肖志懷，MALIK O P. 信號幅值跳變小波降噪方法研究[J]. 信號處理，2014, 30(4): 443-449.

LU Na, XIAO Zhihuai, MALIK O P. Research of signal amplitude changing wavelet de-noising method [J]. Journal of Signal Processing, 2014, 30(4): 443-449.

[23] 譚建豪，章兢. 遺傳算法在模糊設(shè)計中的應(yīng)用[J]. 控制理論與應(yīng)用，2010, 27(4): 501-504.

TAN Jianhao, ZHANG Jing. Application of genetic algorithms in fuzzy design [J]. Control Theory & Applications, 2010, 27(4): 501-504.