使用核磁共振測量來測定流體的速度和特性的方法和設備的製作方法

2023-06-24 09:14:51

專利名稱：使用核磁共振測量來測定流體的速度和特性的方法和設備的製作方法
技術領域：
本發明總的涉及使用核磁共振(NMR)儀器來進行流體特性。
背景技術：
石油和天然氣行業已經開發了多種能夠測定和預測地層特性的工具。在不同類型的工具中，已經證明核磁共振(NMR)儀器是寶貴的。NMR儀器可以用於測定地層特性，諸如孔洞的部分體積和運動流體填充孔洞的部分體積。NMR測井的總的背景在美國專利No.6140817中描述。
核磁共振是在一組選定的具有核磁矩，即，非零自旋量子數的原子核中出現的現象。當這些原子核被放置在磁場(Bo，「Zeeman場」)中時，它們每個以特定的頻率圍繞Bo場的軸線進動，該頻率為Larmor頻率(ωo)，其為每種核素的特有性能，且依賴於原子核的旋磁比(γ)和在原子核的位置處作用的磁場強度(Bo)，即，ωo＝γBo。
除了提供儲層壓力的測量以外，諸如Schlumberger的ModularDyanmics Testing(MDT)工具之類的井眼流體採樣和測試工具可以提供有關儲層流體的類型和特性的重要信息。這些工具可以使用工具機載傳感器模塊來執行井下流體特性的測量。或者，這些工具可以從儲層提取流體試樣，其可以在瓶子中收集，且帶到地面上來分析。收集的試樣定期送到流體特性實驗室，以分析物理特性，其中包括石油粘度、氣油比、質量密度或者API度、分子組成、H2S、瀝青質、樹脂和各種其它雜質濃度。然而，因為試樣可能由泥漿濾液汙染，所以實驗數據可能沒有用，或者不相應於儲層流體特性。
例如，收集的流體試樣可能是過濾水和原油的乳濁液，或者在由基於石油的泥漿鑽井的井中，可能是儲層原油和基於石油的泥漿濾液的混合物。在兩種情況下，汙染物都可能致使測量的實驗數據與就地儲層流體的實際特性不相關。在帶到地面上的試樣具有低的或者可忽略的汙染物的情況下，實驗結果仍然可能被破壞(例如，由於由溫度變化引起的固體沉澱)。因此，需要在井下條件下對沒有汙染的試樣進行地層流體分析。
一些美國專利披露了用於在井眼中對從地層提取的流體試樣進行NMR測量的方法和設備。頒發給Kleinberg的美國專利No.6346813B1和Kleinberg等的公布的美國專利申請No.2003/0006768披露了一種在MDT工具的流送管上的NMR模塊，用於由核磁信號測定不同的流體特性。該Kleinberg專利和Kleinberg出版物轉讓給本發明的受讓人。從而，通過引證整個包括該專利和申請。頒發給Blades等的美國專利No.6111408也披露了可以包括在類似於Schlumberger MDTTM工具的流體採樣工具中的NMR模塊。該NMR模塊可以執行不同的NMR測量，包括T1(縱向馳豫時間)、T2(橫向馳豫時間)以及D(擴散常數)。頒發給Edwards等的美國專利No.6111409披露了一種用於測量來自碳-13原子核的自旋迴波信號的設備和方法，以及由化學移譜來分析這些測量值來測定脂族烴與芬香烴的比例的光譜方法。頒發給Schwede的美國專利No.3528000披露了不同的方法，包括那些使用裸眼封隔器的，以將流體試樣提取到井眼工具中，以便對流體試樣進行NMR測量。然而，這些現有技術專利都沒有披露用於測定流速或者流體的T1的方法。
頒發給King和Santos的美國專利6046587披露了用於測量在管線中流動的多相流體的流速的方法和設備。King和Santos專利教示了通過使用單個NMR傳感器來使用由不同的衰減時間獲得的信號的NMR自由感應衰減(FID)幅度的比例來推測流速。然而，因為靜磁場變化會導致信號衰減太快而不能可靠地檢測，因此使用永磁體來測量FID信號是困難的。此外，用於測定流速的King和Santos方法沒有說明在管道中有流速分布的事實。在另一個實施例中，King和Santos設備包括兩個分開的NMR傳感器。多相流體的流速和流體體積從兩個傳感器中的FID測量值來計算。然而，該計算需要預先知道流體T1分布。
用於測量流過流體採樣工具的流送管的流體速度的兩種方法在2003年1月22日由Madio等提交的序列號為No.10/349011的美國專利申請中描述。該申請轉讓給本發明的受讓人，因此通過引證整個包括。兩種方法都需要NMR設備的靜磁場在流的方向上具有梯度。第一種方法利用在可變等待時間的脈衝序列中的不同的等待時間測量值之間有Larmor頻移的事實。該測量的頻移與流速和梯度成比例。第二種方法基於這樣的事實，即，在奇數回波和偶數回波之間的測量的相差與平均流速和提到成比例的事實。
雖然這些現有技術的方法對於測定地層流體特性是有用的，但是還是需要更好的方法來使用NMR儀器來測定在地層測試工具中流動的流體的特性。

發明內容
在一個方面，本發明的實施例涉及用於通過核磁共振來測定流體的特性的方法。一種根據本發明的一個實施例的用於通過核磁共振來測定流體的特性的方法包括將靜磁場施加到流體；使用包括擾亂脈衝、等待時間和採集脈衝序列的脈衝序列來採集一組有關流體的核磁共振測量值，其中，對於該等待時間，該組核磁共振測量值具有不同的值；以及使該組核磁共振測量值擬合成響應於流體的正向模型，以導出選自流體的流速、縱向馳豫時間及其組合的參數。
本發明的另一個方面涉及使用核磁共振來監測在被提取到地層流體測試工具中的流體中的汙染物的方法。一種根據本發明的一個實施例的使用核磁共振來監測在被提取到地層流體測試工具中的流體中的汙染物的方法包括將靜磁場施加到流體；使用包括擾亂脈衝、等待時間和採集脈衝序列的脈衝序列來採集一組流體的核磁共振測量值，其中，對於該等待時間，該組核磁共振測量值具有不同的值；使該組核磁共振測量值擬合成響應於流體的正向模型，以導出流體的特性；以及根據導出的流體特性監測流體中的汙染物的級別。
本發明的另一個方面涉及用於測量流體特性的NMR儀器。根據本發明的一個實施例的核磁共振設備包括流管，該流管包括預極化部分和探查部分，其中，預極化部分在探查部分的上遊；圍繞流管設置的磁體，以產生覆蓋預極化部分和探查部分的靜磁場；在探查部分處圍繞流管設置的天線，用於產生具有大致垂直於靜磁場的磁偶極子的磁偶極子的振蕩磁場，以及用於接收核磁共振信號；以及用於控制通過天線產生振蕩磁場和接收核磁共振信號的電路。
從下面的說明書、附圖和權利要求書可以明白本發明的其它方面。


圖1示出了包括適用於實現本發明的實施例的NMR模塊的地層流體測試工具。
圖2A示出了根據本發明的一個實施例的NMR傳感器的示意圖。
圖2B示出了沿著根據本發明的一個實施例的NMR傳感器的流管的自旋極化曲線的示意圖。
圖3A示出了根據本發明的一個實施例的用於採集可變等待時間的NMR測量值的脈衝序列的圖。
圖3B示出了根據本發明的一個實施例的包括擾動脈衝的脈衝序列的圖。
圖4A和4B分別示出了層流和非層流的示意圖。
圖5示出了根據本發明的一個實施例在仿真中相對於輸入流速的計算的流速。
圖6示出了根據本發明的一個實施例將信號幅度擬合成一組具有用於不同流速的可變等待時間的數據。
圖7示出了根據本發明的一個實施例將信號幅度擬合成一組具有用於不同流速的可變等待時間的數據。
圖8示出了根據本發明的一個實施例在仿真中對於具有不同的T1流體的相對於輸入流速的計算的流速。
圖9示出了根據本發明的一個實施例在仿真中對於不同的流速的相對於輸入的T1的計算的T1。
圖10示出了根據本發明的一個實施例的用於測定流體的流速和特性的方法的流程圖。
圖11示出了根據本發明的一個實施例的用於監測在地層流體測試工具中的流體中的汙染物的方法的流程圖。
具體實施例方式
本發明的實施例涉及使用NMR測量值來測定流體的速度和特性的方法和設備。這裡披露的用於測定流速的方法不依賴於靜磁場梯度。
本發明的一些實施例還涉及用於測量流體的縱向馳豫時間(T1)的方法。T1的測量值提供了一種用於監測在從儲層提取到諸如由Schlumberger Technology Corporation(Houston，TX)以MDTTM商品名銷售的Modular Dynamics Tester工具或者任何類似的流體採樣工具之類的井眼流體測試/採樣工具中的儲層流體中的汙染物的級別的方法。監測T1隨著時間的變化對於確定何時轉移或者停止流是有用的，以便(1)就地進行縱向馳豫時間(T1)或者橫向馳豫時間(T2)、擴散係數(D)或者化學移譜的NMR測量，或者(2)收集可以帶到地面且送到實驗室來分析的流體試樣。
本發明的一些實施例涉及用於計算流體的粘度的方法。根據本發明的一個實施例，倘若流體的氣油比從其它測量值已知，那麼流體的T1也可以用於使用已知的經驗關係來估計儲層石油粘度。根據本發明的另一種方法可以使用跨過流送管的壓降(ΔP)的測量值連同例如由NMR測定的流速來計算粘度。
本發明的方法可以由從本領域中已知的那些儀器改進的NMR儀器來實踐。該NMR儀器可以在實驗室中，或者是測井電纜工具的一部分，諸如用於分析從地層提取的流體的地層測試器(流體採樣工具)。
圖1示出了地層流體測試/採樣工具，諸如Schlumberger的Modular Dynamics Testing Tool(MDT)，其包括幾個部分，這些部分使得可以從可滲透的地層提取流體。如所示的，地層流體測試工具10包括下面的模塊電子模塊11，其可以包括處理器和存儲器；液壓動力模塊12；探測模塊13，其可以展開來與地層液壓密封；泵出模塊17；光學流體分析器(OFA)14；以及多採樣模塊16。此外，該地層流體測試工具10包括NMR傳感器15。
根據本發明的實施例的NMR傳感器15包括(1)磁體(例如，永磁體)，其設計為在地層流體測試工具10的流送管(流管)中產生靜磁場，以及(2)RF天線(其可以起發射器和接收器的作用)，其設計為發射使其磁偶極子大致垂直於(正交)靜磁場的磁偶極子的振蕩磁場。振蕩磁場的頻率選擇為等於在探查下的對NMR靈敏的原子核(例如，1H或者13C)的Larmor頻率。由於信噪比的考慮，最好測量在快速流動的流體中的1H原子核。對於固定的測量(即，當流體不流動時)，可以測量來自包括13C的其它原子核的信號。本領域中的普通技術人員可以理解，相同的RF天線可以起發射器的作用，以發射振蕩磁場，且可以起接收器的作用，以接收信號，如在頒發給Kleinberg的美國專利No.6346813中披露的。或者，可以使用分開的發射器和接收天線。磁體可以具有預極化區域，用於在流體進入發射器/接收器天線區域以前極化流體。預極化對於從快速運動的原子核獲得可測量的信號是必須的。
圖2A示出了根據本發明的一個實施例的用於測定在流送管22中流動的流體的流速和其它特性的NMR傳感器15的視圖。如所示的，傳感器15包括不導電和非磁性的流送管22、永磁體24和用於發射和接收射頻(RF)信號的天線26。
如所示的，流送管(流管)22包括預極化部分28，其在探查部分29的上遊。磁體24圍繞預極化部分28和探查部分29設置。磁體24可以是永磁體或者電磁體。如所示的，流送管22由不導電和非磁性材料製成，諸如複合或者聚合材料。如果流送管22由導電或者磁性材料(例如，鋼)製成，那麼天線26最好位於流送管內。
對於在圖2A中顯示的傳感器，天線(電磁線圈)26起發射器和接收器的作用。或者，NMR傳感器15可以包括兩個天線一個起發射器的作用，另一個起接收器的作用。此外，天線可以包括電磁線圈、迴路或者鞍形線圈。天線26大致覆蓋流送管22的探查部分29。
根據本發明的方法的NMR測量值包括一組可變等待時間(W)測量值。在每個等待時間(W)以前，磁化首先由設計為」消除」磁化的脈衝來擾亂，使得Mx＝My＝Mz＝0。在擾亂脈衝和等待時間以後，將90度脈衝，隨後是180度脈衝(例如，自旋迴波脈衝)施加到發射器，以產生自旋迴波。
接收天線測量自旋迴波信號的同相和正交分量。測量重複多個不同的等待時間。在本發明的其它實施例中，可以施加多個180度脈衝，以產生多個自旋迴波信號。此外，本領域中的普通技術人員可以理解，自旋迴波還可以使用除了90和180度脈衝以外的脈衝獲得，例如，使用混合脈衝。
對於不同的等待時間的檢測的自旋迴波信號的幅度依賴於流速、等待時間、接收器和發射器長度、磁體預極化長度和流體T1分布。除了流體的流速和T1分布以外，所有這些參數通過傳感器設計或者通過脈衝序列確定。如果可以利用足夠組測量值，那麼這些參數可以通過將測量的信號擬合成合適的模型來推導，該模型仿真響應於流體的NMR。也就是，數據可以通過正向模型來解釋。
下面披露的方法使用了理論的正向模型，其允許對於前述的脈衝序列、傳感器參數、流速和T1分布的任何組合來計算檢測的回波信號。正向模型可以通過轉化來擬合到測量的回波信號，以測定表面流速和T1分布。或者，如果流速從其它測量值已知，那么正向模型可以通過轉化來擬合，以確定T1分布。
傳感器和可變等待時間脈衝序列圖2B示出了與速度和位置相關的極化曲線(f(v，z)與接收器線圈(在圖2A中顯示為26)的長度(la)的曲線由兩部分構成——一部分來自新的自旋，當它們在等待時間(W)期間在靜磁場中行進到接收器中時，它們被」預極化」，以及第二部分，其來自接在收器線圈中且在等待時間(W)期間」預極化」的自旋。在圖2A和2B中顯示的預極化長度(lp)相應於預極化區域28的長度。
與位置和速度相關的極化函數由下面的等式給出，f(v,z)=1-exp(-T(z,v)T1).---(1)]]>其中，T(z，v)是對於具有位置z和速度v的自旋的極化時間，T1是自旋的縱向馳豫時間。極化時間T(z，v)定義為T(z,v)=lp+zvfor,0zvW---(2a)]]>以及，T(z，v)＝W for v·W＜z≤la.(2b)從上述等式可以明白，與速度(v)和位置(z)相關的極化函數f(v，z)也依賴於lp、W和T1。然而，因為極化時間T(z，v)用於簡化表示，所以這些依賴性隱含在等式1中。在討論在圖3A和3B中顯示的可變等待時間(VWT)脈衝序列的討論以後，可以更容易理解這些等式。
圖3A示出了一組包括N個測量值的可變等待時間脈衝。N最好是在10的量級。在該組中的每次測量以後，施加一個或者多個擾亂脈衝(由脈衝S共同地表示)，以在開始下一個測量的等待時間以前破壞任何剩餘的磁化。擾亂脈衝的持續時間和頻率是與儀器有關的，且可以根據經驗來確定。
如在圖3B中由」蝌蚪」示意性顯示的，每個等待時間以後施加90度激勵脈衝，以將縱向磁化旋轉到橫向平面上。由於靜磁場中的不均勻性(或者其它因素)，來自橫向磁化的信號迅速地去相，但是通過180度脈衝被重新聚焦，以產生自旋迴波。在獲得自旋迴波信號以後，施加擾動脈衝來去除接收器線圈內(如圖2A中顯示為26)的自旋的磁化。本領域中的普通技術人員可以認識，像在圖3中顯示的一樣的可變等待時間序列也可以用於觀察自由感應衰減信號，替代自旋迴波。也就是，可以省略在圖3B中顯示的180度脈衝，且可以在90度脈衝以後觀察自由感應衰減信號。只需要對正向模型進行較小的和明顯的修改。本領域中的普通技術人員也可以認識，通過對正向模型的明顯的修改，也可以使用諸如倒轉恢復或者飽和恢復脈衝序列之類的其它T1靈敏的序列。倒轉恢復和飽和恢復脈衝序列在本說明書中通常都稱為」T1馳豫探查脈衝序列」。
用於採集NMR信號(FID或者回波)的脈衝序列在本說明書中通常稱為」採集脈衝序列」。採集脈衝序列可以包括單個90度脈衝、自旋迴波脈衝(即，90度脈衝後面是180度脈衝)，以及包括跟隨90度脈衝的多個180度(重新聚焦)脈衝的自旋迴波脈衝的變體(例如，CPMG)。這樣，如在圖3A和3B中所示的」可變等待時間脈衝序列」包括擾動脈衝、等待時間和採集脈衝序列。自旋迴波脈衝序列及其變體(例如，CPMG)在本說明書中通常稱為」自旋迴波脈衝序列」。也就是，」自旋迴波脈衝序列」不但包括單個90度脈衝和單個180度(重新聚焦)脈衝，而且包括90度脈衝以後的多個180度(重新聚焦)脈衝。
在每個等待時間期間，新極化或者預極化的自旋進入天線區域。等式2a示出了，在區域0≤z≤vW中，新的自旋在等待時間期間進入天線。該區域的長度依賴於流速(v)和等待時間(W)，即，其依賴於流速(v)和等待時間(W)的乘積。這些新的自旋的極化時間與等待時間無關。替代的，由於自旋進入了永磁體，所以其依賴於自旋暴露到靜磁場的持續時間。這是因為，當施加擾動脈衝時，流體的該部分在發射器/接收器天線(在圖2A中顯示為26)的外部，且在流體的該部分中的自旋的磁化沒有被擾動脈衝」消除」。另一方面，如等式2b可見，因為當施加擾動脈衝時，這些自旋在發射器/接收器的區域中，所以在長度為la-vW的鄰接區域中的自旋的預極化由等待時間(W)控制。這些自旋的磁化由擾動脈衝去除，以及由接收器檢測的任何極化是由於在等待時間(W)期間的再極化。如果等待時間足夠長，或者流速足夠快，使得乘積vW超過天線長度，那麼只有新的自旋(那些在擾動脈衝以後進入接收器天線的自旋)被測量，且極化函數與W無關。
用於流體的NMR正向模型的推導這裡披露的方法最好使用正向模型來預測NMR傳感器對於具有給定的流速和馳豫時間T1的分布的流體的響應。傳感器的響應也依賴於等待時間和諸如天線長度、磁體預極化長度和流送管的半徑之類的傳感器的物理參數。對於給定的等待時間和傳感器設計，正向模型中僅有的變量為流速和流體的T1分布。流速和T1分布最好由正向模型的轉化來確定。正向模型在下面的段落中推導。
為了精確地對流體的NMR響應建模，應該考慮在管道中流動的粘性流體的層流的速度曲線是拋物線形的事實。例如，如Victor L.Streeter，「Fluid Mechanics，」McGraw-Hill Book Co.，5thEdition.p.244可見。在層流中，最大流速Vm出現在管子的軸線上，而在管子壁處的速度為零。在圓形管子中的層流方式的特徵在於雷諾數的值R≤2000～3000，其中，層流的確切上限依賴於管子的表面粗糙度。雷諾數由下式限定，R=2rov,---(3)]]>其中，ro是流送管(管子)的半徑，v是平均流速，ρ是流體質量密度，以及η是流體的粘度。與層流相反，湍流具有無序分量，且更難建模。湍流的一個特徵是速度拋物線曲線的變平。圖4示出了在圓形管子中的層流(A)和非層流(B)的速度曲線。
為了開發正向模型，需要表示由在圖3A和3B中顯示的可變等待時間脈衝序列測量的回波幅度。如果Ai表示由等待時間Wi測量的相位校正的回波信號的幅度，那麼Ai可以寫成接收器天線的靈敏體積(Ω)的體積積分的形式，即，Ai=Sa(r){P(T1)f(v(r),z,T1,lp,Wi)dT1}d3r.---(4)]]>在等式4中的內部積分是與速度和位置相關的極化函數的T1權重的積分。權重函數是流體的T1分布P(T1)。體積積分由空間靈敏度函數Sa(r)來加權，該空間靈敏度函數可以由測量或者計算來確定。為了簡化下面的討論，考慮具有均勻的空間靈敏度的接收器是有用的，以便，
Sa=1=1ro2la---(5)]]>其中，ro是流管的半徑。使用等式5和柱面坐標(r，θ，z)，等式4中的積分可以寫成下面的形式Ai=2lar020P(T1){0rO[0laf(v(r),z,T1,lp,Wi)dz]rdr}dT1.---(6)]]>在得到等式6中，假設靈敏度函數具有方位角對稱性來進行角積分。注意，等式6的被積函數中的極化函數由等式1和2a-2b給出。為了進行Ai的評價，與速度相關的極化函數可以關於z坐標積分來獲得，即，f(v(r),T1,lp,la,Wi)=0Wiv(1-exp(-lp+zvT1))dzla+Wivla(1-exp(-WiT1))dzla.---(7)]]>上面的積分可以解析地進行來產生f(v(r),T1,lp,la,Wi)=Wivla-vT1la(1-exp(-WiT1))exp(-lpvT1)+(1-exp(-WiTi))(la-Wivla)--(8a)]]>對於Wi·v＜la，以及對於Wi·v≥la，積分由下式給出f(v(r),T1,lp,la,Wi)=1-vT1la(1-exp(-lavT1))exp(-lpvT1).---(8b)]]>在進一步處理以前，進行一些觀察是有用的，以便提供對等式8a-8b的了解。等式8a中的前兩項說明自旋的預極化，而最後項是由於在等待時間(W)期間沒有通過輸出流從接收器天線區域去除的自旋的再極化。注意，等式8b可以通過設置Wi·v＝la來從等式8a獲得。一旦滿足條件Wi·v≥la，靈敏區域完全充滿新的自旋，且極化函數不依賴於Wi。
在等式8a和8b中的函數表示在z坐標積分以後的徑向靈敏函數。接下來，必須關於徑向坐標來積分這些函數。徑向靈敏函數的徑向依賴性是由在流送管中的流速的徑向變化產生的(如圖4可見)。對於層流，速度曲線是拋物線的，且可以寫成這樣的形式，v(r)=(r2-ro2)ro2vm---(9)]]>
其中，ro是流送管的半徑，vm是如圖4所示的流送管的軸線上(即，在r＝0處)的最大流速。根據等式9，層流的平均流速是vm/2。雖然等式9描述了對於層流普遍使用的模型，但是替代的層流模型可以如下描述v(r)=vm(1-rro)1/n---(9a)]]>其中，n通常在5和10之間。本發明的實施例可以使用這兩個層流模型或者其任何變體中的任何一個。為了簡明起見，下面的描述是基於在等式9中描述的模型。然而，本領域中的普通技術人員可以理解，在下面的描述中，在等式9a中描述的模型可以替代等式9的模型。
考慮下面的關於徑向坐標的積分，F(vm,T1,lp,la,Wi)=2ro20rorf(v(r),T1,lp,la,Wi)dr---(10)]]>因為等式9描述了半徑r和流速v之間的關係，在等式10中顯示的關於半徑r的積分可以轉換為通過使用等式9來關於流速v的積分，以改變變量，例如dv=-2rro2vmdr.---(11)]]>等式9和11可以用於改變等式10中的積分變量和界限，以給出F(vm,T1,lp,la,Wi)=0vmf(v,T1,lp,la,Wi)dvvm.---(12)]]>對等式12中的積分的變量進行另一種改變是有用的。改變等式12中的積分中的變量v＝vm·ξ導致積分F(vm,T1,lp,la,Wi)=01f(vm,T1,lp,la,Wi)d,---(13)]]>其採用可以使用高斯積分方法來精確地數值進行的形式。合併項，回波幅度的等式6可以再寫成下面的形式
Ai=0P(T1)F(vm,T1,lp,la,Wi)dT1.---(14)]]>注意，通過lp，la，和Wi的已知值，模型回波幅度只依賴於P(T1)、分布T1和最大流速vm。如在後面部分中描述的，P(T1)和vm的估計可以通過使用多組可變等待時間數據來轉化模型來獲得。對於層流，vm和流送管半徑ro完全確定等式9中的速度曲線。
在等式13中的表示是基於在等式9中描述的拋物線速度流體曲線來推導的。這對於層流是合適的。然而，對於在圖4中顯示的非層流曲線，在等式13中推導的等式是無效的，應該修改。
在圖4中描述的對於非層流的速度曲線可以通過下面的等式來描述v(r)＝vmfor 0≤r≤α，(15a)以及v(r)={(r-a)2-(ro-a)2}vm(ro-a)2,for,arro.---(15b)]]>使用等式15a和15b，等式10中的積分可以表示為兩個積分的和，即，F(vm,a,T1,lp,la,Wi)=2ro20arf(vm,T1,lp,la,Wi)dr+2ro2arorf(v(r),T1,lp,la,Wi)dr.---(16)]]>右邊的第一個積分是明了的。通過使用等式15b將關於r的積分轉換為關於v的積分來執行第二個積分。最後的步驟是引入變量v＝vm·ξ的變化，以將等式16轉化為下面的形式F(vm,a,T1,lp,la,Wi)=a2ro2f(vm,T1,lp,la,Wi)+(ro-a)2ro201f(vm,T1,lp,la,Wi)d.---(17)]]>注意，如果a＝0，等式17減少到用於層流的等式13。此外，在等式14中的模型回波幅度除了依賴於vm和P(T1)以外還依賴於參數a。通過轉化正向模型，估計參數a，從而推測湍流或者非層流的波端或者存在是可行的。
數值積分等式17中的積分可以表示為被積函數關於一組固定的橫坐標的加權和，即，F(vm,a,T1,lp,la,Wi)=a2ro2f(vm,T1,lp,la,Wi)+(ro-a)22ro2k=1nwkf(vmk,T1,lp,la,Wi).---(18)]]>對於n階高斯積分的權重wk和橫坐標xk可以在1964年6月由M.Abramowitz和L.A.Stegun編輯的National Bureau of StandardsApplied Mathematics Series 55的」Handbook of MathematicalFunctions With Formulas，Graphs，and Mathematical Tables，」中給定的表格中找到。例如，在下面的表格中，顯示6階高斯積分的橫坐標和權重。

表格中給定的橫坐標用於在(-1，1)的範圍中的積分。在等式18中的橫坐標可以使用等式從後面獲得。
k=xk+12.---(19)]]>在進行等式18的總和中，依賴於ξk＜la/vmWi或者ξk≥la/vmWi，再調用」被積函數」是等式8a或者8b中的與速度相關極化函數是有用的。
正向模型的轉化根據本發明的實施例，一旦通過可變等待時間獲得了流體的NMR數據，它們被擬合成正向模型，以提供流體的流速和T1信息。轉化包括改變正向模型中的參數vm、P(T1)和a，直到測量的回波幅度符合由模型預測的回波幅度。有很多不同的方法來將測量的數據擬合成模型。下面描述了兩種方法來作為例子。本領域中的普通技術人員可以理解，在不偏離本發明的範圍的情況下也可以使用其它方法。
轉化方法NO.1測量的NMR數據包括在NMR儀器的同相和正交通道中測量的信號。正交信號相對於在同相通道中檢測的相應的信號為90度異相。
測量的信號包括噪聲。這樣，為了將模型擬合到測量數據，一些測量噪聲應該加到由模型推導的信號中。下面的等式示出了該方法A~i,r=Aicos+Ii,r,---(20a)]]>A~i,q=Aisin+Ii,q.---(20b)]]>其中， 和 分別為對於通過等待時間Wi測量的在同相和正交通道中的測量的幅度；Ai為模型預測的信號幅度；φ為模型預測的信號相位；以及Ii，r和Ii，q分別為信號和正交通道中的噪聲。假設Ii，r和Ii， q為在兩個通道中具有相等的方差σ2的零均值高斯噪聲。使用等式14，信號幅度Ai可以寫成下面的形式Ai=l=1NcalFi,l(vm,a),---(21)]]>其中，在等式14中關於T1的積分由關於用於近似T1分布的Nc幅度(ai)的總和來替代。下面的量定義為Fi，l(vm，a)≡F(vm，a，T1，l，la，lp，Wi).
(22)且再調用包括參數a的等式18，以描述由於非層流產生的拋物線前端變平的半徑。
確定離散組的馳豫時間T1，l，且可以選擇為在對數標度上相等的隔開。因為同相和正交信號相互成90度異相，所以可以從使用下面的等式測量的同相和正交幅度估計的信號相位φ^=arctan(i=1NmA~i,qi=1NmA~i,r).---(23)]]>其中，φ上的脫字符表示其為估計值；以及Nm是在可變等待時間測量值組中的測量值的數量(不同的等待時間)。
然後，根據在頒發給Freedman的美國專利No.5291137(」Freedman專利」)中披露的方法，從等式23獲得的相位估計值可以用於計算相位校正的信號加上噪聲幅度和均方根(rms)噪聲。該專利通過引證在此整個包括。簡要地，A~i(+)=A~i,rcos^+A~i,qsin^---(24a)]]>以及A~i(-)=A~i,rsin^-A~i,qcos^.---(24b)]]>Freedman專利顯示，A~i(+)>=Ai---(25)]]>其中，角括號表示期望值，以及Ai是如在等式14中描述的模型預測的幅度。可以從該等式計算rms噪聲估計值，^=1Nmi=1Nm(A^i(-))2.---(26)]]>通過最小化擬合誤差的平方和，例如通過相對於參數vm，、a和該組幅度a1來最小化目標函數，測量的幅度數據可以擬合成模型幅度
E(vm,a{al})=i=1Nm(A~i(+)-Ai(vm,a,{al}))2+l=1Nmal2,---(27)]]>在等式27中的最後一項用於通過減小在未確定的逆問題中出現的噪聲幅度的影響來穩定該方法。參數γ是調整參數，且γ≥0。
注意，在永磁體的磁場梯度中的」自旋動力學」效果可以導致測量的回波幅度的衰減。如Hürlimann and Griffin，「Spin-Dynamicsof Carr-Purcell-Meiboom-Gill-like Sequences in GrosslyInhomogeneous B0 and B1 Fields and Application to NMR WellLogging，「J.Magn.Reson.v.143，120-135(2000)可見。因此，對於自旋動力學衰減校正了等式27中的回波幅度。
等式27中的最小化可以使用許多市場上可獲得的優化程序來進行，在服從約束條件的情況下，這些優化程序可以最小化線性或者非線性函數。一種這樣的程序是稱為」DLCONG」的雙精度子程序，其可以從由Visual Numerics Corp.(San Ramon，CA)銷售的商品名為IMSLTM的優化代碼庫的一部分可以得到。如果解析函數可以設置為目標函數的偏倒數，那麼該DLCONG代碼更加穩固。目標函數的倒數可以解析地計算。例如，從等式27可以獲得下面的倒數Eak=-2i=1Nm(A~i(+)-Ai(vm,a,{al}))Fi,k(vm,a)+2ak,---(28)]]>以及Evm=-2i=1Nm(A~i(+)-Ai(vm,a,{al}))k=1NmakFi,k(vm,a)vm,---(29)]]>以及
Ea=-2i=1Nm(A~i(+)-Ai(vm,a,{al}))k=1NmakFi,k(vm,a)a.---(30)]]>目標函數在等式29中的倒數需要Fi，l相對於Vm.的微分。根據等式17，該微分需要與速度相關的極化函數的倒數，即，Fi,l(vm,a)vm=a2ro2fi,l(vm)vm+(ro-a)22ro2k=1nwkfl,l(vmk)vm.---(31)]]>在等式31的求導中，積分和微分的階互換，且與速度相關的極化函數的緊縮記法引入如下fi，l(vm·ξ)≡f(vmξ，T1，l，la，lp，Wi)，(32)其可以通過替代v＝vm·ξ來從等式8a和8b獲得。因為轉化決定性地依賴於這些倒數的精度，所以顯示這些詳細的結果是有用的。等式31的總和的倒數如下給定，當ξ＜la/(vm·Wi)時，fi,l(vm)vm=Wlla-[(lp+vmT1,l)vmlaexp(-lpvmT1,l)+Wila](1-exp(-WiTi,l)),--(33)]]>以及當ξ≥ la/(vm·Wi)時，fi,l(vm)vm=-(1-exp(-lavmT1))(lp+vmT1,l)vmlaexp(-lpvmT1,l)+1vmexp(-la+lpvmT1,l)---(34)]]>注意，等式31的右手側的第一項的偏倒數可以通過設置ξ＝l來從等式33和34獲得。等式30中的Fi，l相對於a的偏倒數是明了的，且通過等式17的微分可以容易地獲得。
轉化方法No.2在前面部分中描述的方法將相位校正的數據擬合成馳豫模型。下面描述的替代方法同時將兩通道數據擬合到正向模型。類似於等式27的誤差或者目標函數可以限定為
E(vm,a,{al,r},{al,q})=i=1Nm(A~l,r-Al,r(vm,a,{al,r}))2+i=1Nm(A~l,q-Ai,q(vm,a,{al,q}))2+l=1Nmal2,---(35)]]>其中，Ai，r＝Aicosφ，(36a)Ai，q＝Aisinφ，(36b)al，r＝alcosφ，(36c)以及al，q＝alsinφ.
(36d)等式35的目標函數可以相對於vm，a以及兩組幅度{al，r}和{al， q}最小化。T1分布中的幅度可以使用下面的等式從這兩組幅度直接計算，al=al,r2+al,q2.---(37)]]>信號相位也可以使用下面的等式從這些幅度估計，^=arctan(i=1Ncal,ql=1Ncal,r),---(38)]]>
其中，Nc是用於近似T1分布的分量數。此外，等式35中的目標函數的最小化可以通過許多市場上可得到的程序(例如，DLCONG)中的任何一種以及前面部分中討論的偏倒數來實現。
數值仿真本發明的實施例可以精確地估計流體中的流速和T1值(流體的T1值的分布或者平均)。接下來描述一些計算機仿真實驗，其證明了vm和T1可以使用本發明的實施例來估計的精度。對於具有la＝0.5英寸、lp＝10.0英寸和ro＝0.1英寸的NMR傳感器(如圖2所示)進行仿真。在從0.001秒到0.72秒的範圍的可變等待時間組中有35個等待時間。用於仿真的平均流速為2.5，5.0，10.0，15.0，20.0，25.0，30.0，35.0和40.0英寸/秒。
圖5示出了平均流速(＝vm/2)相對於估計的平均流速的比較。該仿真是在假設0.1p.u.的rms噪聲幅度的情況下進行的。考慮到當在完全極化的固定流體上測量時信號幅度為100p.u.，這是相對高的信噪比。在圖5中顯示的仿真是對於具有T1值的範圍的三種不同的流體進行的。用於仿真中的T1分布是單個定值的(即，Dirac δ函數)。注意，對於所有三種流體，在用於仿真中的估計的平均流速(y軸線)和輸入值(x軸線)之間有好的一致性。在較高的平均流速下(例如，大約30英寸/秒以上)有一些精度損失。然而，考慮到高的流速，這些估計值相當的好。
圖6示出了對於具有T1＝0.316秒的流體的作為流速的函數的可變等待時間數據組，以及正向模型的轉化後擬合到這些數據。注意，在圖6中，由於不完全的極化，測量的回波幅度隨著增加的流速而較大地減小。這可以說明在圖5中的流速估計值在較高的流速下的精度損失。該效果對於具有長T1的流體加劇。
圖7示出了對於具有T1＝3秒的流體在不同流速下的可變等待時間數據，以及轉化後擬合。擬合的質量通過如下定義的無量綱x平方的」擬合質量因數」參數來評價，2=1Nw2i=1Nw(A~i(+)-Ai(vm,a,{al}))2.---(39)]]>
對於完全擬合x2＝1，在噪聲誤差範圍內。
注意，回波幅度隨著流速增加而較大地減小。這轉化為測量信噪比的較大地減小。圖7中的數據的另一個重要的特徵是當流速增加時的數據壓縮(即，曲線更靠近在一起)。這樣，減小了對在高流速處進行的流速測量的相對靈敏度。這兩種效果，減小的信噪比和減小的靈敏度，意味著在較高的流速處精確地轉化可變等待時間數據組更加困難。
圖8示出了與圖5中相同的仿真的結果，但是由1.0p.u.的rms噪聲來替代0.1p.u.。圖5和8的比較顯示，流速估計的質量由信噪比的十倍降低而較大地下降。這清楚地顯示了具有高信噪比測量的重要性。
對於具有寬T1分布的流體，分布的幾何平均或者分布的其它矩可以用於監測和定性流送管中的流體。
從測量值估計流體粘度如上所述，本發明的一些實施例涉及從流體的NMR測量值估計流體粘度的方法。一旦流體的流速根據上述的方法由具有可變等待時間的一組NMR測量值確定，流體的粘性也可以根據在下面的部分中描述的方法來估計。將描述使用對在流送管中的流體進行的可變等待時間的測量來估計流體粘度(η)的兩種方法。本領域中的普通技術人員可以理解，這兩種方法只是示例性的，在不偏離本發明的範圍的情況下可以使用其它方法。
第一種方法使用由NMR測量值估計的平均流速(v)和跨過流送管長度L的測量的壓降ΔP來計算粘度(η)的估計。壓降ΔP可以從裝備在NMR傳感器上或者地層測試工具的壓力傳感器來獲得。從流速(v)和壓降(ΔP)的流體粘度計算可以根據熟知的流體力舉公式來實現(例如，如Victor L.Streeteretal.，FluidMechanics，5thEdition.p.244，McGraw-Hill，NY可見)=Pro28vL,---(40)]]>其中，流體的平均速度v等於vm/2。因為L和ro由設備限定，所以等式40可以轉化為更加普遍的形式，
=KPv,---(41)]]>其中，比例常數K可以由對具有以已知的壓力梯度在不同的速度下測量的已知粘度的流體進行的測量值來經驗地確定。使用v，η和流體的質量密度的獨立測量值，可以從等式3計算雷諾數R。等式40或者41中的計算對於任何粘性的流體是有效的，且不限制於原油或者鹽水。
用於估計原油粘度的第二種方法是使用T1分布的幾何平均或者算術平均。T1分布的算術平均通過已知的相關性來與石油粘度相關(例如，如Freedman et al.，「A New NMR Method of FluidCharacterization in Reservoir RocksExperimentalConfirmation and Simulation Results，」SPE Journal，December2001，452-464可見)，o=kTT1,LMf(GOR),---(42)]]>其中，ηo是釐泊的石油粘度，k是對於原油的經驗確定的常數，T是Kelvin溫度，T1，LM是石油T1分布的算術平均，以及f(GOR)是溶解氣油比(GOR)經驗確定的函數。等式42中GOR的可以從流送管中的其它測量值來確定，或者可以從儲層流體特性的局部或者先驗知識已知。
由等式40、41或者42計算的流體的計算的T1分布或者粘度可以用於監測流送管中的流體的汙染物的變化。該應用將在下面的部分中描述。
圖10示出了根據本發明的實施例的方法100。根據該方法，通過可變等待時間採集一組NMR數據(步驟102)。用於數據採集的脈衝序列可以包括擾動脈衝、等待時間和採集脈衝序列，如圖3A和3B所示。採集脈衝序列可以包括簡單的90度激勵脈衝、自旋迴波脈衝(90-滯後-180)、Carl-Purcell-Meiboom-Gill(CPMG)脈衝或者自旋迴波脈衝或者CPMG脈衝的任何變體。單個自旋迴波脈衝在圖3B中示例性顯示。在不偏離本發明的範圍的情況下，顯示的90度和180度脈衝可以由其它脈衝(例如，混合脈衝)來替代。
如上所述，除了在圖3A和3B中顯示的可變等待時間脈衝序列以外，本發明的實施例還可以使用其它T1探查脈衝序列。例如，倒轉恢復脈衝序列或者飽和恢復脈衝序列可以替代地使用。如果使用飽和恢復脈衝序列模型擬合過程類似於使用擾亂脈衝序列。然而，如果使用倒轉恢復脈衝序列，那麼將修改等式1的自旋極化函數。
一旦採集了該組可變等待時間的NMR，它們可以使用上述的轉化方法的任何一種來擬合成正向模型(步驟104)。該轉化過程產生流體的流速(v)和T1值。這樣獲得的T1值可以是T1分布或者T1分布的算術平均。
由NMR數據擬合到流體模型獲得的流速(v)和/或T1值可以用於計算流體的粘度(步驟106)。如上所述，根據等式40或者41，估計的流速(v)連同跨過NMR傳感器的長度L的獨立測量的壓降(ΔP)可以用於計算流體粘度(η)。類似的，估計的T1值或者其算術平均可以連同獨立確定的GOR函數f(GOR)一起使用，以計算流體粘度(η)，如等式42所示。然而，可以使用這些方法的任何變體。
如上所述，當從地層提取流體時，本發明的實施例可以用於監測流體汙染物。由於根據在圖10中顯示的方法能夠計算流體的T1和粘度，所以當流體被提取到地層流體測試工具中時可以監測地層流體的純度。
圖11示出了根據本發明的方法的用於這樣的監測的方法110。根據該方法，地層流體被提取到測試工具中(步驟112)，且監測流體的粘度或者T1值(步驟114)。因為提取的流體的早期部分可能被泥漿濾液汙染，所以它們的特性(例如，粘度和T1值)可能不同於未汙染的地層流體。當地層流體的流繼續時，泥漿濾液汙染的級別將逐漸降低。因此，通過監測流體的粘度或者T1值，可以辨別何時泥漿濾液汙染物減小到可忽略的水平，例如通過監測反應純地層流體的特性的近似穩態值的粘度或者T1值。
通過訊問粘度或者T1值是否到達穩定值(步驟116)，可以決定是否繼續提取地層流體。例如，如果粘度或者T1值沒有穩定，那麼提取過程和監測將繼續。如果粘度或者T1值穩定，那麼提取的流體試樣可以保存下來，用於以後的分析，或者可以進行就地井下測量(步驟118)。
各種NMR測量已經在井下進行。例如，頒發給Freedman的美國專利No.6229308和頒發給Venkataramanan等的美國專利No.6462541B1披露了在井下條件下的NMR測量。這些測量可以提供T2(橫向馳豫時間)分布、地層流體的D(擴散)分布、T2-D圖、流體分類和流體粘度(η)。
在圖11中顯示的方法110功能上等於目前用於監測由地層流體測試工具提取的流體中的汙染物的很多監測技術。其它方法例如可以監測流體的顏色、光密度或者GOR。這樣，在圖11中顯示的方法110提供了這些現有方法的替代物。
如圖10和11中顯示的方法可以實現為電腦程式，且存儲在與測井一起使用的計算機或者處理器上，例如，包括在圖1中顯示的電子模塊11中的處理器或者電路。
本發明的優點包括如下。本發明的實施例提供了估計流體的流速、T1值和粘度的方便的方法。本發明的方法可以用在地層流體測試工具上或者任何其它類型的設備和需要監測流體流動的應用中。如果用在地層測試工具上，當從地層取出地層流體時，本發明的方法可以用於監測地層流體的汙染物。這些方法確保了取出基本沒有汙染物的地層流體來用於分析。通過提供當地層流體被提取時用於監測地層流體中的汙染物的方法，本發明的方法可以減輕與地層流體分析中的汙染物相關的問題。此外，在監測表示汙染物可以忽略時，本發明的方法還提供了井下NMR測量。能夠在接近儲層條件下測量流體特性可以提供表示就地流體特性的流體特性。
雖然相對於有限數量的實施例描述了本發明，得益於本披露物的本領域中的普通技術人員可以理解，可以在不偏離如這裡披露的本發明的範圍的情況下設計其它實施例。因此，本發明的範圍只由後附的權利要求書限定。
權利要求
1.一種用於通過核磁共振來測定流體的特性的方法，其包括將靜磁場施加到流體；使用包括擾亂脈衝、等待時間和採集脈衝序列的脈衝序列來採集一組有關流體的核磁共振測量值，其中，對於該等待時間，該組核磁共振測量值具有不同的值；以及使該組核磁共振測量值擬合成響應於流體的正向模型，以導出選自流體的流速、縱向馳豫時間及其組合的參數。
2.如權利要求1所述的方法，其特徵在於，採集脈衝序列包括選自自旋迴波脈衝序列和單個脈衝中的一種。
3.如權利要求1所述的方法，其特徵在於，通過正向模型的轉化來進行擬合。
4.如權利要求1所述的方法，其特徵在於，還包括基於推導的流速和跨過流體行進的管子的選擇的長度的壓降來估計流體的粘度。
5.如權利要求4所述的方法，其特徵在於，根據選自=Pro28vL]]>和=KPv]]>的一個來估計，其中，η是粘度，v是流體的平均速度，L是管子的選擇的長度，ΔP是選擇的管子長度上的壓降，以及r0是管子的半徑，以及K是經驗上確定的常數。
6.如權利要求1所述的方法，其特徵在於，還包括基於推導的流體的縱向馳豫時間和氣油比來估計流體的粘度。
7.如權利要求6所述的方法，其特徵在於，該估計是根據o=kTT1,LMf(GOR)]]>其中，ηo是粘度，k是對於流體的經驗確定的常數，T是Kelvin溫度，T1，LM是流體的縱向馳豫時間的算術平均，以及f(GOR)是氣油比(GOR)經驗確定的函數。
8.一種用於通過核磁共振來測定流體的特性的方法，其包括將靜磁場施加到流體；使用包括縱向馳豫探查脈衝序列和採集脈衝序列的脈衝序列來採集一組有關流體的核磁共振測量值，其中，對於縱向馳豫探查脈衝內的衰減時間，該組核磁共振測量值具有不同的值；以及使該組核磁共振測量值擬合成響應於流體的正向模型，以導出選自流體的流速、縱向馳豫時間及其組合的參數。
9.如權利要求8所述的方法，其特徵在於，該縱向馳豫探查脈衝選自倒轉恢復脈衝序列和飽和恢復脈衝序列中的一種。
10.如權利要求8所述的方法，其特徵在於，該採集脈衝序列包括選自自旋迴波脈衝序列和單個脈衝中的一種。
11.如權利要求8所述的方法，其特徵在於，通過正向模型的轉化來進行擬合。
12.如權利要求8所述的方法，其特徵在於，還包括基於推導的流速和跨過流體行進的管子的選擇的長度的壓降來估計流體的粘度。
13.如權利要求12所述的方法，其特徵在於，根據選自=Pro28vL]]>和=KPv]]>的一個來估計，其中，η是粘度，v是流體的平均速度，L是管子的選擇的長度，ΔP是選擇的管子長度上的壓降，以及r0是管子的半徑，以及K是經驗上確定的常數。
14.如權利要求8所述的方法，其特徵在於，還包括基於推導的流體的縱向馳豫時間和氣油比來估計流體的粘度。
15.如權利要求14所述的方法，其特徵在於，該估計是根據o=kTT1,LMf(GOR)]]>其中，ηo是粘度，k是對於流體的經驗確定的常數，T是Kelvin溫度，T1，LM是流體的縱向馳豫時間的算術平均，以及f(GOR)是氣油比(GOR)經驗確定的函數。
16.一種使用核磁共振來監測在被提取到地層流體測試工具中的流體中的汙染物的方法，其包括將靜磁場施加到流體；使用包括擾亂脈衝、等待時間和採集脈衝序列的脈衝序列來採集一組流體的核磁共振測量值，其中，對於該等待時間，該組核磁共振測量值具有不同的值；使該組核磁共振測量值擬合成響應於流體的正向模型，以導出流體的特性；以及根據導出的流體特性監測流體中的汙染物的級別。
17.如權利要求16所述的方法，其特徵在於，流體的特性包括選自縱向馳豫時間的分布、縱向馳豫時間的算術平均，以及它們的組合中的一種。
18.如權利要求16所述的方法，其特徵在於，流體的特性是粘度
19.一種核磁共振設備，其包括流管，該流管包括預極化部分和探查部分，其中，預極化部分在探查部分的上遊；圍繞流管設置的磁體，以產生覆蓋預極化部分和探查部分的靜磁場；在探查部分處圍繞流管設置的天線，用於產生具有大致垂直於靜磁場的磁偶極子的磁偶極子的振蕩磁場，以及用於接收核磁共振信號；以及用於控制通過天線產生振蕩磁場和接收核磁共振信號的電路。
20.如權利要求19所述的設備，其特徵在於，電路包括具有使用包括擾動脈衝、等待時間和採集脈衝序列的脈衝序列來用於採集一組流體的核磁共振測量值的指令的程序。
21.如權利要求20所述的設備，其特徵在於，採集脈衝序列包括選自自旋迴波脈衝序列和單個脈衝中的一種。
22.如權利要求20所述的設備，其特徵在於，該程序還包括用於使該組核磁共振測量值擬合成響應於流體的正向模型，以導出選自流體的流速、縱向馳豫時間及其組合的參數的指令。
23.如權利要求22所述的設備，其特徵在於，通過正向模型的轉化來進行擬合。
24.如權利要求22所述的設備，其特徵在於，該程序還包括用於基於推導的流速或者推導的縱向馳豫時間來估計流體的粘度的指令。
全文摘要
一種用於通過核磁共振來測定流體的特性的方法包括將靜磁場施加到流體；使用包括擾亂脈衝、等待時間和採集脈衝序列的脈衝序列來採集一組有關流體的核磁共振測量值，其中，對於該等待時間，該組核磁共振測量值具有不同的值；以及使該組核磁共振測量值擬合成響應於流體的正向模型，以導出選自流體的流速、縱向馳豫時間及其組合的參數。
文檔編號E21B47/13GK1591042SQ20041006874
公開日2005年3月9日 申請日期2004年9月6日 優先權日2003年9月5日
發明者R·弗裡曼 申請人:施盧默格海外有限公司