專利名稱：用于水下電纜的雙傳感器降噪系統(tǒng)的制作方法
用于水下電纜的雙傳感器降噪系統(tǒng)
背景技術(shù)：
本發(fā)明總體上涉及海洋地震勘探，并且更具體地涉及降低拖在勘探船后的傳感器中或位于海底的傳感器中不期望的地震波反射的影響的設(shè)備和方法。在拖曳式海洋地震勘探中，在近海面的船舶后拖拽著水聽器陣列。位于多個傳感器電纜中的水聽器通常被稱為拖纜。同樣拖在近海面的震源周期性地發(fā)射聲能。該聲能通過海洋向下傳播，從底層結(jié)構(gòu)反射，然后通過海洋向上返回到水聽器陣列。該水聽器陣列記錄來自海床向上傳播的地震聲波。隨后，將水聽器的記錄處理成為底層結(jié)構(gòu)的地震圖像。聲阻抗是媒介密度P和媒介中聲速c的乘積P Co任何時候產(chǎn)生的反射都會改變聲波遇到的聲阻抗。聲阻抗變化越大，反射能量越多。由于空氣和水的聲阻抗差別很大，因此海面是一個近似理想的聲能反射器。從海底或感興趣目標(biāo)返回之后，該能量再次由海面 反射回拖纜。由于水聽器具有全向響應(yīng)，因此水聽器陣列記錄重影響應(yīng)，該重影響應(yīng)是從海面反射并最后到達(dá)延時到達(dá)的地震聲波，且由于直接反射而顛倒極性。該重影是向下傳播的地震聲波，當(dāng)添加到希望的波形上時該地震聲波會有損記錄的地震圖像。該重影在水聽器響應(yīng)的頻譜中產(chǎn)生一個陷波，fnotch=c/2d，其中c是聲速，并且d是拖纜深度。按照慣例，地震拖纜拖曳在10米或更淺的深度。在10米深度，陷波頻率(fnotch)是75Hz。對于高地震圖像分辨率，頻率響應(yīng)必須擴(kuò)展到IOOHz之外。因此，拖纜有時會拖曳在較淺的深度以提高地震圖像的分辨率。該重影引起的反射還可延續(xù)到海底或其他強(qiáng)反射物，并向上返回再次干擾期望的反射波并降低圖像質(zhì)量。這些反射通常稱為多次波。由于來自海面的噪聲干擾了期望的地震信號，因此在淺水中拖曳是有問題的。此夕卜，近海面的循環(huán)水流會在拖纜外殼導(dǎo)致流動噪聲。這些影響隨天氣變壞而惡化，有時會導(dǎo)致工作人員中斷操作，直到天氣變好。拖曳越深，海面噪聲和天氣影響越小。如果可以消除重影陷波影響，則令人希望的是在更深處拖曳。海底或海床、將地震電纜或傳感器放置在海床的系統(tǒng)通過普遍共知的技術(shù)(如p-z總和)來抑制重影。在聲波中，壓力P是個標(biāo)量，而粒子速度U是個矢量。水聽器采用正全向響應(yīng)來記錄地震聲波壓力P。垂直定向地震檢波器或加速計采用向上信號的正響應(yīng)和向下信號的負(fù)響應(yīng)來記錄地震聲波粒子速度的垂直分量UZ。在P-Z求和中，速度信號由海水的聲阻抗P C來衡量，并添加到壓力信號。如果使用加速計，可集成其輸出以獲得速度信號，或可以區(qū)分水聽器信號，如此以來，可以更好的與加速計信號光譜匹配。這可生成具有對向上傳播波的全響應(yīng)和對向下傳播波的零響應(yīng)以抑制重影和多次波的組合傳感器。在由Monk等人發(fā)明的美國專利號為6，539，308的專利中描述了一種實現(xiàn)單一去重影痕跡的信號調(diào)節(jié)和信號整合方法。由于不是期望信號導(dǎo)致的因素，當(dāng)粒子速度傳感器或加速計不受多余的動作影響時，該技術(shù)和類似技術(shù)有效。當(dāng)有強(qiáng)大的底部電流時，這類多余的加速度通常位于部署在海浪帶或區(qū)域的海床系統(tǒng)中。目前，使用水聽器和粒子運(yùn)動測量的組合來降低它們在地震拖纜中的影響已經(jīng)引起了大家的興趣。因為該拖纜受到由拖曳或海面效應(yīng)引起的加速度比由期望反射引起的加速度大，因此操作地震拖纜中的粒子運(yùn)動傳感器存在問題。此外，這些多余的加速度與期望的反射響應(yīng)位于相同的光譜帶中。地震拖纜及海床地震電纜經(jīng)歷從O度到360度的所有側(cè)傾角和適度的螺旋角。為實現(xiàn)垂直定向的地震檢波器，海底系統(tǒng)已經(jīng)使用(a) —個萬向動圈式地震檢波器；(b) —個3-組件、全傾斜動圈式地震檢波器，該檢波器具有與該傳感器無關(guān)的姿態(tài)感知和計算以求解與重力有關(guān)的測量；以及(c ) 一個3-組件、微電機(jī)系統(tǒng)(MEMS )加速計，該加速計具有與該傳感器無關(guān)的內(nèi)部姿態(tài)感知和計算以求解與重力有關(guān)的測量。授予Rouquette的美國專利號7，167，413在地震拖纜中使用加速計來抑制重影-陷波效應(yīng)。Rouquette使用質(zhì)量_彈簧系統(tǒng)來降低電纜動力對加速計和負(fù)載傳感器系統(tǒng)的影響以測量和抑制該加速計上的電纜運(yùn)動感應(yīng)噪聲。Rouquette系統(tǒng)依賴眾所周知的復(fù)雜機(jī)械關(guān)系，該關(guān)系不能與制造公差、老化和環(huán)境狀況保持一致。Rouquette使用信號處理自適應(yīng)算法來推導(dǎo)負(fù)載傳感器和質(zhì)量彈簧系統(tǒng)與作用于原地加速計上的加速度的關(guān)系Rouquette描述了一種復(fù)雜的機(jī)械電子系統(tǒng)。
由Tenghamn等人發(fā)明的美國專利號為7，239，577描述了一種使用聲波粒子速度傳感器抑制重影陷波的設(shè)備和方法。Tenghamn等人傳授了流體阻尼、萬向架固式地震檢波器的使用。本領(lǐng)域眾所周知的是選擇封裝該地震檢波器的流體來提供懸掛在其平衡架上的傳感器的阻尼。然而，本領(lǐng)域眾所周知而在Tenghamn等人中沒有描述的是質(zhì)量_彈簧隔振系統(tǒng)可降低電纜機(jī)械運(yùn)動對地震檢波器響應(yīng)的影響。由電纜機(jī)械運(yùn)動引起的地震檢波器的運(yùn)動與地震檢波器響應(yīng)中的聲波粒子運(yùn)動是不可區(qū)分的。期望的地震波粒子運(yùn)動被Tenghamn等人的電纜機(jī)械運(yùn)動所掩蓋。由Vaage等人發(fā)明的美國專利號7，359，283描述了一種結(jié)合負(fù)載傳感器和粒子運(yùn)動傳感器來解決機(jī)械運(yùn)動對粒子運(yùn)動傳感器的影響的方法。在該方法中，在某一頻率fO以下不使用粒子運(yùn)動傳感器的響應(yīng)，而只是從壓力傳感器響應(yīng)和已知的壓力傳感器深度來估計。這些抑制的頻率是拖纜的機(jī)械運(yùn)動所期望的。在更低的感興趣頻率上，估計響應(yīng)具有差的信噪比。某一頻率以下的這種抑制不是最優(yōu)的，因為它還抑制了重要低頻段中的有用信號，在該低頻段內(nèi)可能存在深度目標(biāo)數(shù)據(jù)。盡管以上提到的專利都描述了抑制地震拖纜中的重影陷波的方法，但是所有這些專利都缺乏充分說明拖纜繩和水流導(dǎo)致的運(yùn)動對粒子運(yùn)動傳感器的影響。所有這些專利還缺乏生成高保真、具有低至感興趣的最低頻率的良好信噪比的感應(yīng)聲波成分。概述這些不足可由體現(xiàn)本發(fā)明特點的水下電纜的噪聲抑制系統(tǒng)來克服。此類系統(tǒng)的一個版本包括布置在水下電纜中的粒子運(yùn)動傳感器。這些傳感器生成原始傳感器信號，該信號包括對由地震和電纜運(yùn)動引起的粒子運(yùn)動的響應(yīng)。一種自適應(yīng)處理裝置使用水下電纜的自適應(yīng)水動力模型從原始傳感器信號生成估計傳感器信號。該估計傳感器信號表示沒有地震時水下電纜對電纜運(yùn)動的響應(yīng)的估計。從原始傳感器信號中減去估計傳感器信號的裝置用于確定對由地震引起的粒子運(yùn)動的響應(yīng)。水下電纜噪聲降低系統(tǒng)的另一個版本包括布置在水下電纜中的粒子運(yùn)動傳感器。該傳感器生成原始傳感器信號。自適應(yīng)處理裝置使用水下電纜的自適應(yīng)水動力模型來從原始傳感器信號中生成估計傳感器信號。對地震靈敏的地震探測器在地震期內(nèi)禁用自適應(yīng)處理裝置。這樣，該估計傳感器信號表示沒有地震時粒子運(yùn)動傳感器對噪聲源的響應(yīng)。在本發(fā)明的另一個方面，一種降低水下電纜中的噪聲的方法，包括(a)使用水下電纜的自適應(yīng)水動力模型來從原始傳感器信號計算由電纜運(yùn)動產(chǎn)生的估計傳感器信號，該原始傳感器信號由布置在水下電纜上的粒子運(yùn)動傳感器提供；以及(b)通過從原始傳感器信號中減去估計傳感器信號來確定對由地震引起的粒子運(yùn)動的響應(yīng)。附圖
簡要說明通過參照以下描述、所附權(quán)利要求和附圖可更好的理解本發(fā)明的這些方面和特征，其中
圖I是體現(xiàn)本發(fā)明特征的水下電纜噪聲降低系統(tǒng)的框圖。詳細(xì)說明在圖I的框圖中，水聽器10生成水聽器信號H。該水聽器可實現(xiàn)為安裝在水下傳感器電纜(比如，拖曳的拖纜或連接到勘探船的海底電纜)中的一個單一的信號水聽器或一組水聽器。本例中的水聽器是一個加速度消除水聽器，該水聽器對壓力變化靈敏，但是對由傳感器電纜動態(tài)變化引起的加速度不靈敏。與傳感器電纜中的水聽器大致共同定位的是第
二傳感器，粒子運(yùn)動傳感器14——本例中，是三軸加速計——測量由地震波反射引起的粒子加速度。該粒子運(yùn)動傳感器可以是任何對粒子速度或加速度響應(yīng)的傳感器。加速計、地震波探測器、加速度敏感水聽器、類似此類傳感器組或不同此類傳感器的組合是粒子運(yùn)動傳感器實現(xiàn)的其他例子。粒子運(yùn)動傳感器14生成原始傳感器信號A，該信號包括對由地震弓I起的粒子運(yùn)動和對電纜運(yùn)動的響應(yīng)。由于傳感器電纜動態(tài)旋轉(zhuǎn)的可能性，與粒子運(yùn)動傳感器14相關(guān)聯(lián)的方向傳感器16用于確定與重力矢量有關(guān)的粒子運(yùn)動測量的方向。多軸傳感器(比如，對直流作出響應(yīng)的MEMS加速計)結(jié)合粒子運(yùn)動感應(yīng)和單個部件中的方向感應(yīng)。傳統(tǒng)計算裝置可用于確定與重力有關(guān)的粒子運(yùn)動測量的方向。如果方向的一個軸已知，而沿著該軸的電纜運(yùn)動不是一個影響因素，則該粒子運(yùn)動傳感器可以是一個兩軸單元。因為與拖纜平行的軸可以由其他裝置(比如導(dǎo)航系統(tǒng))來確定，因此兩軸傳感器可用于拖曳的拖纜中。另一個替代方案是使用能夠自我修復(fù)的萬向架固定式粒子運(yùn)動傳感器(比如萬向架固定式檢波器)。在萬向架固定式粒子運(yùn)動傳感器中，方向傳感器不是必需的。另一個替代方案是與該粒子運(yùn)動傳感器聯(lián)合使用的單獨(dú)的方向傳感器。無論如何，用于該框圖和貫穿本說明書的原始傳感器信號A指的是重力參考測量。地震探測器18分析該水聽器信號H以確定是否存在地震，比如反射的地震信號。無論何時水聽器信號的振幅超過了預(yù)定頻譜內(nèi)的預(yù)定門限，該地震探測器輸出地震信號
20。地震的實際探測可由電纜中的電子設(shè)備或探測船上水聽器信號的分析來完成。在數(shù)據(jù)采集或后采集處理期間，可實時完成探測。典型地，地震導(dǎo)航系統(tǒng)28 (比如，美國德克薩斯州的休斯頓離子地球物理公司的Spectra或Orca)還用在海洋數(shù)據(jù)采集期內(nèi)。該地震導(dǎo)航系統(tǒng)提供了一個地震源觸發(fā)信號19。該地震導(dǎo)航系統(tǒng)28可以將該地震源觸發(fā)信號傳遞到地震探測器18以提高探測可用性。自適應(yīng)水動力模型22可用作自適應(yīng)處理裝置以估計粒子運(yùn)動傳感器10對電纜運(yùn)動的響應(yīng)。可使用Kalman濾波器作為自適應(yīng)處理裝置的一部分來改進(jìn)該模型。當(dāng)沒有探測到地震時，原始傳感器信號A與估計傳感器信號A’在減法框24內(nèi)進(jìn)行比較，該減法框構(gòu)成從原始傳感器信號中減去估計傳感器信號的裝置。A-A’的差用于更新自適應(yīng)處理裝置中Kalman濾波器的協(xié)方差矩陣。自適應(yīng)的其他選項包括遞歸最小二乘法、神經(jīng)網(wǎng)絡(luò)、或模糊邏輯。如果確定所期望的地震響應(yīng)使該模型惡化，假設(shè)在水聽器和加速計處同時探測所期望的地震響應(yīng)，在地震期間通過地震信號20禁用模型自適應(yīng)，如交換器26所示。模型應(yīng)用和自適應(yīng)可通過布置在拖纜中、船上或陸地上的計算裝置來實現(xiàn)?？稍跀?shù)據(jù)采集或后處理期間實時完成。視情況需要，可使用可用于該模型的其他信息來提高水動力模型的精確度。這可包括實時船舶動力、GPS數(shù)據(jù)、艦船速度的測量，或與傳感器電纜動力有關(guān)的任何其他測量。船用導(dǎo)航系統(tǒng)28估計探測到地震能量位置處的拖纜位置、速度和加速度。所計算的導(dǎo)航數(shù)據(jù)29可用于提高水動力模型的精度。該水動力模型還可由該導(dǎo)航系統(tǒng)用于增加支持其估計的數(shù)據(jù)，如導(dǎo)航數(shù)據(jù)29的雙向流所示。由于在地震期間通過自適應(yīng)動態(tài)模型使得地震探測器18無法自適應(yīng)，水動力模型(估計傳感器信號A’)的輸出將表示沒有所期望的地震反射時加速計對電纜運(yùn)動和其他噪聲源的響應(yīng)。在用于減法的第二裝置30中產(chǎn)生估計傳感器信號A’與原始傳感器信號A之間的差。該差表示對粒子運(yùn)動的降低噪聲響應(yīng)P，這是所期望的地震響應(yīng)。通過形成結(jié)合 對粒子運(yùn)動P的響應(yīng)和水聽器信號H的裝置，PZ求和32用于抑制重影陷波或由多余的多次波所產(chǎn)生的響應(yīng)以生成去重影地震響應(yīng)信號S，該水聽器信號可被區(qū)分從而光譜地匹配原始粒子運(yùn)動傳感器信號A。盡管參照單一版本的變化詳細(xì)描述了本發(fā)明，但是其他版本是可能的。例如，結(jié)合水聽器和粒子運(yùn)動測量的方法可通過PZ求和來實現(xiàn)，或通過正在使用或研究的降低多次波或重影陷波效應(yīng)的任何其他方式來實現(xiàn)。此外，可實時或在后處理中離線完成單個水聽器和粒子運(yùn)動測量的調(diào)整和縮放及所有其他計算。這可包括由入射角或聲阻抗所引起的粒子運(yùn)動測量的縮放。還可包括信號的區(qū)分或集成，如此以來它們可以被光譜地匹配。因此，如這些少數(shù)例子所建議，詳細(xì)描述的版本旨在幫助舉例說明而非不限制本發(fā)明。
權(quán)利要求
1.一種用于水下電纜的噪聲抑制系統(tǒng)，包括 一個布置在水下電纜中且生成原始傳感器信號的粒子運(yùn)動傳感器，該原始傳感器信號包括對由地震和電纜運(yùn)動所引起的粒子運(yùn)動的響應(yīng)； 使用水下電纜的自適應(yīng)水動力模型來從原始傳感器信號生成估計傳感器信號的自適應(yīng)處理裝置，其中該估計傳感器信號表示在沒有地震時水下電纜對電纜運(yùn)動的響應(yīng)的估計； 從該原始傳感器信號中減去該估計傳感器信號來確定對由地震引起的粒子運(yùn)動的響應(yīng)的裝置。
2.如權(quán)利要求I所述的噪聲抑制系統(tǒng)，進(jìn)一步包括對地震靈敏且在地震期內(nèi)禁用該自適應(yīng)處理裝置的地震探測器。
3.如權(quán)利要求2所述的噪聲抑制系統(tǒng)，進(jìn)一步包括加速度消除水聽器，該加速度消除水聽器向該地震探測器發(fā)送水聽器信號以從該水聽器信號中探測地震。
4.如權(quán)利要求3所述的噪聲抑制系統(tǒng)，其中該加速度消除水聽器和該粒子運(yùn)動傳感器協(xié)同定位于該水下電纜上。
5.如權(quán)利要求I所述的噪聲抑制系統(tǒng)，進(jìn)一步包括生成水聽器信號的水聽器及結(jié)合對粒子運(yùn)動的響應(yīng)和該水聽器信號來生成地震響應(yīng)信號的裝置。
6.如權(quán)利要求I所述的噪聲抑制系統(tǒng)，其中該自適應(yīng)水動力模型從該原始傳感器信號和該估計傳感器信號之間的差值中導(dǎo)出該估計傳感器信號。
7.如權(quán)利要求I所述的噪聲抑制系統(tǒng)，進(jìn)一步包括向該自適應(yīng)水動力模型提供導(dǎo)航數(shù)據(jù)的船用導(dǎo)航系統(tǒng)。
8.如權(quán)利要求I所述的噪聲抑制系統(tǒng)，進(jìn)一步包括接收來自該適應(yīng)水動力模型的數(shù)據(jù)以提高由該船用導(dǎo)航系統(tǒng)所計算的導(dǎo)航數(shù)據(jù)的質(zhì)量的船用導(dǎo)航系統(tǒng)。
9.如權(quán)利要求2所述的噪聲抑制系統(tǒng)，進(jìn)一步包括向該地震探測器提供地震源觸發(fā)信號的船用導(dǎo)航系統(tǒng)。
10.如權(quán)利要求I所述的噪聲抑制系統(tǒng)，其中，該自適應(yīng)處理裝置實時操作。
11.如權(quán)利要求I所述的噪聲抑制系統(tǒng)，其中，該自適應(yīng)處理裝置對所存儲的原始傳感器信號進(jìn)行離線操作。
12.—種降低水下電纜中的噪聲的方法,包括 使用水下電纜的自適應(yīng)水動力模型從由布置在水下電纜上的粒子運(yùn)動傳感器提供的原始傳感器信號計算由電纜運(yùn)動導(dǎo)致的估計傳感器信號； 通過從該原始傳感器信號中減去該估計傳感器信號來確定對由地震引起的粒子運(yùn)動的響應(yīng)。
13.如權(quán)利要求12所述的方法，進(jìn)一步包括 探測地震；以及 在地震期間禁止該估計傳感器信號的計算。
14.如權(quán)利要求13所述的方法，進(jìn)一步包括，從來自與該粒子運(yùn)動傳感器協(xié)同定位的水聽器的水聽器信號中探測地震。
15.如權(quán)利要求12所述的方法，進(jìn)一步包括，結(jié)合對粒子運(yùn)動的響應(yīng)和來自與該粒子運(yùn)動傳感器協(xié)同定位的加速度消除水聽器的水聽器信號以生成地震響應(yīng)信號。
16.如權(quán)利要求12所述的方法，進(jìn)一步包括，從船用導(dǎo)航系統(tǒng)向該自適應(yīng)水動力模型提供導(dǎo)航數(shù)據(jù)。
17.如權(quán)利要求12所述的方法，進(jìn)一步包括，使用該估計傳感器信號以提高由該船用導(dǎo)航系統(tǒng)所計算的導(dǎo)航數(shù)據(jù)的質(zhì)量。
18.如權(quán)利要求12所述的方法，該方法實時操作。
19.如權(quán)利要求12所述的方法，進(jìn)一步包括，離線使用該自適應(yīng)水動力模型來計算該估計傳感器信號。
20.一種用于水下電纜的噪聲抑制系統(tǒng)，包括 一個布置在水下電纜中且生成原始傳感器信號的粒子運(yùn)動傳感器； 使用該水下電纜的自適應(yīng)水動力模型從該原始傳感器信號生成估計傳感器信號的自適應(yīng)處理裝置； 對地震靈敏且在地震期內(nèi)禁用該自適應(yīng)處理裝置的地震探測器，其中該估計傳感器信號表示在沒有地震時該粒子運(yùn)動傳感器對噪聲源的響應(yīng)。
21.如權(quán)利要求20所述的噪聲抑制系統(tǒng)，進(jìn)一步包括，從該原始傳感器信號中減去該估計傳感器信號來確定對由地震引起的粒子運(yùn)動的響應(yīng)的裝置。
22.如權(quán)利要求20所述的噪聲抑制系統(tǒng)，進(jìn)一步包括生成水聽器信號的水聽器及結(jié)合對由地震引起的粒子運(yùn)動的響應(yīng)和該水聽器信號來生成地震響應(yīng)信號的裝置。
23.如權(quán)利要求20所述的噪聲抑制系統(tǒng)，進(jìn)一步包括，向該地震探測器發(fā)送水聽器信號以從該水聽器信號中探測地震的加速度消除水聽器。
24.如權(quán)利要求23所述的噪聲抑制系統(tǒng)，其中該加速度消除水聽器和該粒子運(yùn)動傳感器協(xié)同定位于該水下電纜上。
25.如權(quán)利要求20所述的噪聲抑制系統(tǒng)，其中該自適應(yīng)水動力模型從該原始傳感器信號和該估計傳感器信號的差值中推導(dǎo)出該估計傳感器信號。
26.如權(quán)利要求20所述的噪聲降低系統(tǒng)，進(jìn)一步包括，向該自適應(yīng)水動力模型提供導(dǎo)航數(shù)據(jù)的船用導(dǎo)航系統(tǒng)。
27.如權(quán)利要求20所述的噪聲抑制系統(tǒng)，進(jìn)一步包括，從該自適應(yīng)水動力模型接收數(shù)據(jù)以提高由該船用導(dǎo)航系統(tǒng)所計算的導(dǎo)航數(shù)據(jù)的質(zhì)量的船用導(dǎo)航系統(tǒng)。
28.如權(quán)利要求20所述的噪聲抑制系統(tǒng)，其中該自適應(yīng)處理裝置實時操作。
29.如權(quán)利要求20所述的噪聲抑制系統(tǒng)，其中，該自適應(yīng)處理裝置對所存儲的原始傳感器信號進(jìn)行離線工作。
全文摘要
一種抑制水下傳感器電纜(比如，拖纜或海底電纜)中的噪聲的系統(tǒng)和方法。該電纜的自適應(yīng)水動力模型(22)從來自粒子運(yùn)動傳感器(比如加速計(14))的原始傳感器信號中生成估計傳感器信號。該估計傳感器信號表示在沒有地震時水下電纜對電纜運(yùn)動的響應(yīng)的估計。通過從原始傳感器信號中減去(24)估計傳感器信號以抑制原始傳感器信號中的電纜運(yùn)動和其他噪聲影響來生成單獨(dú)對粒子運(yùn)動的噪聲降低響應(yīng)。地震探測器(18)使用來自加速度抑制水聽器(10)的水聽器信號以在地震期內(nèi)禁止(26)適配水動力模型。通過PZ求和(32)將該水聽器信號與對粒子運(yùn)動的響應(yīng)相結(jié)合，以生成去重影的地震響應(yīng)信號。
文檔編號G01V1/36GK102792188SQ201180012675
公開日2012年11月21日 申請日期2011年1月19日 優(yōu)先權(quán)日2010年1月19日
發(fā)明者D·J·蘭伯特 申請人:離子地球物理學(xué)公司