f (X，t) = -Φf (X)ΣG[Φf (x+e;)-il^x+ejilei，Φf (X) =exp(_l/ Pf(x,t)),itf(x+ei) =exp(-l/pf(x+ei;t))；
[0045] 依據(jù)公式（7)和（8)分別獲得所述第一類組分g在位置x、時刻t的平衡態(tài)的分布 函數(shù)和所述第二類組分f在位置X、時刻t的平衡態(tài)的分布函數(shù)
[0048] 依據(jù)公式（9)和（10)分別獲得所述第一類組分g在位置（x+cAt)、時刻（t+At) 的分布函數(shù)gi(χ+cΛt,t+Λt)和所述第二類組分f在位置（x+cΛt)、時刻（t+Λt)的 分布函數(shù)fi(χ+cΛt,t+Λt):
[0051]其中，c2=c!/3;
[0052] 依據(jù)公式（11)和（12)分別獲得所述第一類組分g的濃度Pg和所述第二類組分 f的濃度Pf:
[0055]獲得所述磁化矢量的模值，所述磁化矢量的模值為所述第一類組分g的濃度Pg。
[0056] 在本發(fā)明實(shí)施例中，對時間進(jìn)行離散化處理以使流體的時間步長為Λt，同時對三 維數(shù)據(jù)體的X軸、Y軸和Z軸進(jìn)行離散劃分獲得多個晶格，該晶格優(yōu)選為立方體，即孔隙中 的流體具體以晶格的形式進(jìn)行劃分，并且將流體的組分劃分為第一類組分g和第二類組分 f，所述第一類組分g為所述晶格內(nèi)在XY平面內(nèi)對應(yīng)有磁化矢量分量的流體分子即能夠?qū)?CPMG脈沖序列產(chǎn)生反應(yīng)的被標(biāo)記的組分，所述第二類組分f為所述晶格內(nèi)除所述第一類組 分g之外的流體分子即不能對CPMG脈沖序列產(chǎn)生反應(yīng)的未被標(biāo)記的組分，在初始狀態(tài)下， 合理假設(shè)第一類組分g的濃度初始值為1-10，第二類組分f的濃度初始值為10 '且第一 類組分g的濃度初始值和第二類組分f的濃度初始值不限于此，所述第一類組分g在平衡 態(tài)的宏觀速度》;為〇,第二類組分f在平衡態(tài)的宏觀速度為〇,通過公式（1)計算獲得第 一類組分g在平衡態(tài)的分布函數(shù)= %夂，此處pg= 1-10 '通過公式（2)計算獲得第二 類組分f在平衡態(tài)的分布函數(shù)f7q=wlPf，此處Pf= 10 '其中，ie[1，2，....，N]，N-1 表示與所述晶格相鄰的所有晶格的數(shù)目，Wl表示權(quán)重系數(shù)，Pg表示所述第一類組分g的 濃度，Pf表示所述第二類組分f的濃度，（^表示晶格聲速，^表示離散速度單位矢量，且
，_ *^_表不所述第一類組分g在 平衡態(tài)的宏觀速度，表示所述第二類組分f在平衡態(tài)的宏觀速度，第一類組分g的濃度 初始值為pg= 1-10，第二類組分f的濃度的初始值為pf= 10 '
[0057] 公式（3)中第一類組分g的宏觀速度ug的初始值為0,公式⑷中第 二類組分f的宏觀速度uf的初始值為0,則公式（5)中ug(x，t)初始值為0， <(以)=#:>，/.)//7辦 值為 0，，Pf(x，t) =Pf= 10，其中，F(xiàn)g(x，t) =-Φg (χ)ΣG[Φg (χ+e；) -Φf (χ+e；) ]e；,Φg (x) =exp(-1/pg (x,t)),Φg (χ+e；) =exp(-1/ Pg(x+ei，t))，τ和G為常數(shù)，F(xiàn)f(x，t) =-itf(x)EGl^fU+eJ-il^x+eJilei，itf(x)=exp(-1/pf (x,t))，Φf (χ+e;) =exp(-1/pf (x+e;，t));通過公式（5)和（6)分別獲得所述 第一類組分g在位置X、時刻t的平衡態(tài)的宏觀速度<%^〇和所述第二類組分f在位置X、 時刻t的平衡態(tài)的宏觀速度^
[0058] 將公式（5)計算獲得的《丨?.Μ)和pg(x,t) =pg= 1-10 1()帶入公式（7)獲得第 一類組分g在位置X、時刻t的平衡態(tài)的分布函數(shù)，將公式（6)計算獲得的<1夂,0 和Pf(x，t)=Pf=10w帶入公式⑶獲得第二類組分f在位置X、時刻t的平衡態(tài)的分 布函數(shù)fT (X，t);其中，cs表示晶格聲速，wi表示權(quán)重系數(shù)，ei表示離散速度單位矢量。
[0059] 將公式⑴獲得的#和公式⑵獲得的〇,?)帶入公式（9)獲得第一類組分 g在位置（x+cΛt)、時刻（t+Λt)的分布函數(shù)gjx+cΛt,t+Λt)，將公式（2)獲得的 fT和公式⑶獲得的fVq (x，t)帶入公式（10)獲得第二類組分f在位置（x+cΛt)、時刻 (t+Λt)的分布函數(shù)fjx+cΛt，t+Λt);其中，c2 二</3。
[0060] 將公式（9)獲得的gl (x+cΛt，t+Λt)帶入公式（11)獲得第一類組分g的濃度 pg，將公式（10)獲得的仁(Χ+CΛt，t+Λt)帶入公式（12)獲得第二類組分f的濃度pf。 并且第一類組分g的濃度P#卩為磁化矢量的模值。
[0061] 本發(fā)明實(shí)施例提供了單相多組分模型即Shan-Doolen模型計算磁化矢量的模值 的方法。
[0062] 圖2為本發(fā)明實(shí)施例提供的CPMG脈沖序列不意圖。在上述實(shí)施例的基礎(chǔ)上，所述 依據(jù)公式（11)和（12)分別獲得所述第一類組分g的濃度pg和所述第二類組分f的濃度 Pf之后，還包括：
[0063] 依據(jù)公式（13)和（14)分別計算所述流體在表面弛豫作用下，所述第一類組分g 的濃度Pg隨時間的變化率sg和所述第二類組分f的濃度Pf隨時間的變化率sf，并依據(jù) 公式（15)和（16)分別獲得表面弛豫作用后的所述第一類組分g的濃度P' g和所述第二 類組分f的濃度P' f:
[0066] Prg=P g-sgPg (15)
[0067] Prf=Pf-sfPf (16)
[0068] 依據(jù)公式（17)和（18)分別計算所述流體在體弛豫作用下，所述第一類組分g的 濃度P'g隨時間的變化率bg和所述第二類組分f的濃度P'f隨時間的變化率bf，并依 據(jù)公式（19)和（20)分別獲得體弛豫作用后的所述第一類組分g的濃度P" g和所述第二 類組分f的濃度P" f:
[0071] p" g=p'g-bgp，g (19)
[0072]p"f=Prf-bfPrf (20)
[0073]依據(jù)公式（21)和（22)分別計算所述流體的擴(kuò)散系數(shù)DjPDf，并依據(jù)公式（23)和 (24)分別獲得流體擴(kuò)散后的所述第一類組分g的濃度Pg和所述第二類組分f的濃 度P"，f:
[0078] 將P"，斿為公式⑴中的Pg、P"，f作為公式⑵中的P,，將<令，_作 為公式⑴中的<，將作為公式⑵中的<，循環(huán)執(zhí)行公式（1)-(24)。
[0079] 在本發(fā)明實(shí)施例中，第一類組分g的濃度與第二類組分f的濃度之和為1，在流體 中加如圖2所示的CPMG脈沖序列，其中，CPMG脈沖序列的90度脈沖是指以氫原子的拉莫 爾進(jìn)動頻率將磁化矢量從z軸正方向搬轉(zhuǎn)到xy平面內(nèi)所需要的脈沖強(qiáng)度和時間；CPMG脈 沖序列180度脈沖是指以氫原子的拉莫爾進(jìn)動頻率將磁化矢量從xy平面內(nèi)軸對稱反轉(zhuǎn)所 需要的脈沖強(qiáng)度和時間。在流體中加CPMG脈沖序列的90度脈沖之前，設(shè)置流體中第一類 組分g的濃度初始值為1-10 1(\第二類組分f的濃度初始值為10 1(\在流體中加CPMG脈沖 序列的90度脈沖之后，設(shè)置流體中第一類組分g的濃度初始值為10，第二類組分f的濃 度初始值為1-10 '
[0080] 流體在表面弛豫作用下，第一類組分g的濃度Pg會隨時間發(fā)生變化，變化率為Sg， 第二類組分f的濃度Pf會隨時間發(fā)生變化，變化率為Sf，將公式（11)獲得的Pg帶入公式 (13)中可獲得變化率\，將9,和、帶入公式（15)獲得表面弛豫作用后，第一類組分g的 濃度P'g;同理，將公式（12)獲得的Pf帶入公式（14)中可獲得變化率sf，將pf和sf帶 入公式（16)獲得表面弛豫作用后，第二類組分f的濃度p'f。
[0081] 流體在體弛豫作用下，第一類組分g的濃度Pg會隨時間發(fā)生變化，變化率為bg， 第二類組分f的濃度Pf會隨時間發(fā)生變化，變化率為bf，將公式（11)獲得的pg帶入公式 (17)獲得bg，將公式（12)獲得的Pf帶入公式（18)獲得bf，將公式（15)獲得的P'， 公式（17)獲得的\帶入公式（19)獲得弛豫作用后的所述第一類組分g的濃度P",;同 理，獲得所述第二類組分f的濃度P" f。
[0082]通過公式（21)和（22)分別計算所述流體的擴(kuò)散系數(shù)DjPDf，并依據(jù)公式（23)和 (24)分別獲得流體擴(kuò)散后的所述第一類組分g的濃度Pg和所述第二類組分f的濃 度P" ' f。