測井地層評價與巖石物理建模研究PowerLog中文培訓（試用版）輝固地球科技（）市朝陽區(qū)東四環(huán)中路56號A1207-1209室100025：01059081801傳真：01059081802測井地層評價與巖石物理建模研究PowerLog中文培訓（試用版）輝固地球科技（）市朝陽區(qū)東四環(huán)中路56號A1207-1209室100025：01059081801傳真：01059081802：@Page1?Fugro-Jason2012..目錄Chapter1. 測井地層評價和巖石物理建模工作流程 3Chapter2. 工區(qū)的建立及數(shù)據(jù)管理 4（一）數(shù)據(jù)庫建立 4工區(qū)的建立 5數(shù)據(jù)加載和管理 7測井繪圖 10（二）（三）（四）Chapter目錄Chapter1. 測井地層評價和巖石物理建模工作流程 3Chapter2. 工區(qū)的建立及數(shù)據(jù)管理 4（一）數(shù)據(jù)庫建立 4工區(qū)的建立 5數(shù)據(jù)加載和管理 7測井繪圖 10（二）（三）（四）Chapter3. 測井資料的預處理 12（一）單井測井資料質量12（二）多井一致性處理 22Chapter4. 測井資料地層評價 29（一）泥質含量的計算 29孔隙度的計算 34飽和度的計算 40ShalySandCP-9測井解釋 42Quicklook快速測井解釋 44曲線統(tǒng)計報告和總結報告 45（二）（三）（四）（五）（六）Chapter5. StatMin最優(yōu)化測井解釋 49Chapter6.巖石物理建模與分析 53（一）巖石物理基礎 54儲層流體置換 60彈性曲線的正演 68巖石物理模板 71（二）（三）（四）Page2Chapter1. 測井地層評價和巖石物理建模工作流程藏特征的數(shù)據(jù)體進行標定。之間缺乏一致性或者數(shù)據(jù)缺失。因此，在和資料進行聯(lián)合應用之前，必須對測井資料進行質量石物理建模作為建立油藏特征和巖石彈性參數(shù)之間 的使得可以根據(jù) 反演得到的彈性參數(shù)，對油藏特征進行解釋，其相互如圖1所示。在項目研究過，首先對測井資料開展質量檢查和質量工作，主要對的自然、中子、聲波和密度曲線進行校正處理，并開展多井一致性檢查和標準化處理工作。在獲得了縱向質量可靠、橫向具有一致性的測井資料后，根據(jù)后續(xù)地層評價和巖石物理參數(shù)建模工作的要求，對地層的巖石物性參數(shù)，如礦物含量、孔隙度等進行精細評價，并對參數(shù)結果的一致性進行檢查和成果質量。1Page3Chapter1. 測井地層評價和巖石物理建模工作流程藏特征的數(shù)據(jù)體進行標定。之間缺乏一致性或者數(shù)據(jù)缺失。因此，在和資料進行聯(lián)合應用之前，必須對測井資料進行質量石物理建模作為建立油藏特征和巖石彈性參數(shù)之間 的使得可以根據(jù) 反演得到的彈性參數(shù)，對油藏特征進行解釋，其相互如圖1所示。在項目研究過，首先對測井資料開展質量檢查和質量工作，主要對的自然、中子、聲波和密度曲線進行校正處理，并開展多井一致性檢查和標準化處理工作。在獲得了縱向質量可靠、橫向具有一致性的測井資料后，根據(jù)后續(xù)地層評價和巖石物理參數(shù)建模工作的要求，對地層的巖石物性參數(shù)，如礦物含量、孔隙度等進行精細評價，并對參數(shù)結果的一致性進行檢查和成果質量。1Page3Chapter2. 工區(qū)的建立及數(shù)據(jù)管理（一）數(shù)據(jù)庫建立如果您的機器上已經(jīng)，那么Start→All2所示：2Chapter2. 工區(qū)的建立及數(shù)據(jù)管理（一）數(shù)據(jù)庫建立如果您的機器上已經(jīng)，那么Start→All2所示：2DBfile為目前處于激活狀態(tài)的數(shù)據(jù)庫，DBhost為計算機名，數(shù)據(jù)庫管理界面分別有CreateandSelect、ManageUsersandProjects、BackupandRestoreDBTools四個選項，在CreateandSelect在ManageUsersandProjects選項中，可以對處于激活狀態(tài)下的數(shù)據(jù)庫和連接到此數(shù)據(jù)庫的用戶進行，在BackupandRestore選項中，可以對數(shù)據(jù)庫進行備份和恢復等，在DBTools選項里，可以對數(shù)據(jù)庫加密、壓縮、輸出、輸入、刪除等。如果您選擇已經(jīng)建立了的數(shù)據(jù)庫，選擇ConnecttoDBonthecomputer...，鏈接到目標路徑即可，如果新建數(shù)據(jù)庫，選擇“CreatenewDB…”33Page4數(shù)據(jù)庫管理計算機名字激活狀態(tài)的數(shù)據(jù)庫按鈕，指定數(shù)據(jù)庫的圖4，數(shù)據(jù)庫建立完成之后，會提示您是否將此數(shù)據(jù)庫激活，點擊“Yes”。4這樣目前處于激活狀態(tài)的數(shù)據(jù)庫按鈕，指定數(shù)據(jù)庫的圖4，數(shù)據(jù)庫建立完成之后，會提示您是否將此數(shù)據(jù)庫激活，點擊“Yes”。4這樣目前處于激活狀態(tài)的數(shù)據(jù)庫DBfile即為daqing_2011.fjdb，如圖5所示。5（二）工區(qū)的建立.PowerLog有兩種方式，可以選擇r→owrog2（圖6，或者雙擊DesktopPowerLog圖標。Page562.工區(qū)建立LoadPoject的框，注意數(shù)據(jù)庫里面是空的（圖7，點擊62.工區(qū)建立LoadPoject的框，注意數(shù)據(jù)庫里面是空的（圖7，點擊rjctks下面的CrtaNewrjct建立一個新的工區(qū)。7直接選擇File→Project→New…，建立新的工區(qū)。也可以8（圖9rjctNme命名為qgNxtPage6910，10交匯處，東經(jīng)135度2分30秒；最西端為帕米我國最東端為是黑龍江和爾（烏恰縣3度0910，10交匯處，東經(jīng)135度2分30秒；最西端為帕米我國最東端為是黑龍江和爾（烏恰縣3度03度52分；最北端為漠河以北黑龍（漠河縣3度33經(jīng)73oE-135oE，選擇研究區(qū)合適的投影坐標系，除了UTM，還有多種投影格式可供選擇，OK，工區(qū)建立完成。（三）數(shù)據(jù)加載和管理支持多種數(shù)據(jù)的加載，PowerLog的批量加載。同時可以對多種類型井相關數(shù)據(jù)進行管理，如測井數(shù)據(jù)、巖心數(shù)據(jù)、生產(chǎn)數(shù)據(jù)、地質數(shù)據(jù)、地理數(shù)據(jù)、圖像數(shù)據(jù)等，數(shù)據(jù)庫對工區(qū)井數(shù)與曲線銷”和“重做”的讀寫操作，箱工具還可以恢復已經(jīng)刪除的文件。主要有三個主要窗口：PowerBenchmainwindow、ProjectExplorer和WellPowerLog。ProjectExplorer對研究區(qū)的井WellInformation則對目標井的井曲線、井分層、曲線別名等進行管理。Page71.測井曲線數(shù)據(jù)加載LASPowerBenchMainWindow界面File→Import→LASFile…LASImporter框（圖1。.11點擊Browse1.測井曲線數(shù)據(jù)加載LASPowerBenchMainWindow界面File→Import→LASFile…LASImporter框（圖1。.11點擊Browse按鈕選擇要加載的數(shù)據(jù)，您可以選擇“Create&Load”建立一口新井并加載曲線一步完成或者先“CreateNewWell”創(chuàng)建一口新井，然后“LoadCurves”加載曲線。對于一口井可能有多次測量的測井曲線需要加載，對于同一種曲線要分不同的名字加以區(qū)分，可以通過名字前面加前綴或者后面加后綴實現(xiàn)（如圖11紅框。如RUN1的井徑曲線定義為CLRUN2的井徑曲線定義為CLGrp加以管理，但要注意的是使用前后綴的時候兩次測量的DEPTH曲線要一致，即DEPTH_1改為DEPTH。12Explorer出現(xiàn)工區(qū)的所有井的，WellInformation顯示目標井中的所有的井曲線12所示。2.地質分層數(shù)據(jù)加載自PowerLog3.3版本以后，地質分層的加載有三種方式，一種是在Wellfrmnxcl表格里面的s13wrg與Excel間的交互操作非常方便，支持Ctrl+C和Ctrl+V的粘貼操作。第二種是框，輸入加載分層數(shù)據(jù)。Page813Files→EditTopsTopsEditor（圖14143.井位坐標加載13Files→EditTopsTopsEditor（圖14143.井位坐標加載Excel表格里面的井坐標數(shù)據(jù)過來粘貼即可。154.巖心數(shù)據(jù)加載數(shù)據(jù)格式進行加載，F(xiàn)ile→Import→ASCIIFile…，打框（圖16mrtUs，加載時默認為數(shù)據(jù)的第一行，需要進行重新指正。16Page9（四）測井繪圖支持多種形式曲線成圖，常規(guī)測井曲線的回放、波形圖像、FMIPowerLog電成像、核磁共振測井T2譜、單井和多井直方圖、單井和多井交會圖、數(shù)據(jù)的互顯示等。交17選擇Demo1井，View→Logplot，或者點擊快捷鍵 ，即打開了繪圖模板，在模板上點擊鼠標右鍵，選擇“Format”，即跳出Format（四）測井繪圖支持多種形式曲線成圖，常規(guī)測井曲線的回放、波形圖像、FMIPowerLog電成像、核磁共振測井T2譜、單井和多井直方圖、單井和多井交會圖、數(shù)據(jù)的互顯示等。交17選擇Demo1井，View→Logplot，或者點擊快捷鍵 ，即打開了繪圖模板，在模板上點擊鼠標右鍵，選擇“Format”，即跳出Format：Logplot可以對Depth、Tops、Track及Curves和圖層等進行設置。（圖18，框中我們18在GeneralDrawOrder可以對圖層進行設置。Page10分層編輯網(wǎng)格線深度19在Track選項里，可以對道進行設置， 刻度、線型、顏色、寬度等。20在Curves選項里，可以對曲線進行設置，19在Track選項里，可以對道進行設置， 刻度、線型、顏色、寬度等。20在Curves選項里，可以對曲線進行設置，所屬道、刻度、線型、顏色等。設置編輯完成之后，點擊O，模板就做好了（如圖21。21以后我們可以把這些模板保存起來。選擇File→SaveScreen，保存為splice_rawcurves。點Page11擊鼠標右鍵，可以通過選擇Zoom→ZoomAllWell和PickDepth對其進行整體和局部放大。Chapter3. 測井資料的預處理在儲層研究項目當中，測井數(shù)據(jù)主要被用來計算縱波和橫波阻抗，同反演作的可靠性，是建立在測井資料的質量基礎之上。因此單井測井資料的質量，以及多井測井資料一致性的質量立可靠的巖石物理模型和確保高質量，在測井資料分析過 反演的重要因素。（一）單井測井資料質量現(xiàn)場測量到的測井資料由于受到儀器本身、測量環(huán)境、地層因素等影響，或多或少的都會如各測井項目間的深度匹配良好、測井曲線盡量完整連續(xù)、各測井項目必要的井眼和泥漿影響校正、聲波等測井曲線的異常尖峰有較一致性響應等。單井測井資料的質量深度上的響應特征是否匹配。對密度質量的影響、聲波周波跳躍對時差的影響、軟地層對橫波時差的影響、擴徑對的影響以及泥漿侵入影響的問題等。測井曲線拼接面向的儲層項目研究需要我們提供完整連續(xù)且質量較淺至深進行，盡量減少和避免套管和水泥環(huán)的影響。在實際測井資料的拼接工作中經(jīng)常會遇到這樣的情況（1（擊鼠標右鍵，可以通過選擇Zoom→ZoomAllWell和PickDepth對其進行整體和局部放大。Chapter3. 測井資料的預處理在儲層研究項目當中，測井數(shù)據(jù)主要被用來計算縱波和橫波阻抗，同反演作的可靠性，是建立在測井資料的質量基礎之上。因此單井測井資料的質量，以及多井測井資料一致性的質量立可靠的巖石物理模型和確保高質量，在測井資料分析過 反演的重要因素。（一）單井測井資料質量現(xiàn)場測量到的測井資料由于受到儀器本身、測量環(huán)境、地層因素等影響，或多或少的都會如各測井項目間的深度匹配良好、測井曲線盡量完整連續(xù)、各測井項目必要的井眼和泥漿影響校正、聲波等測井曲線的異常尖峰有較一致性響應等。單井測井資料的質量深度上的響應特征是否匹配。對密度質量的影響、聲波周波跳躍對時差的影響、軟地層對橫波時差的影響、擴徑對的影響以及泥漿侵入影響的問題等。測井曲線拼接面向的儲層項目研究需要我們提供完整連續(xù)且質量較淺至深進行，盡量減少和避免套管和水泥環(huán)的影響。在實際測井資料的拼接工作中經(jīng)常會遇到這樣的情況（1（2）（3）（4）兩次測量沒有重疊。重疊處曲線值有誤差的這種情況，考慮是不是兩次測量的儀器的曲線，將空缺處補齊。下面以GR_1和GR_2為例，首先粗略估計拼接點的位置，拼接點的深度大概在2752m左右。選擇Edit→splice，鼠標將變?yōu)閹в蠸plice的字樣，點擊圖頭GR_1曲線，這樣GR_1曲線就顯示在SP_DISPLAY道中，GR_2曲線，GR_2曲線也顯示在SP_DISPLAY道中，上下移動拼接點的位置，到合適位置為止（如圖22ceCreGR，同時點擊Table，TableSplice_Table。，Page1222點擊按鈕執(zhí)行cGR曲線完成。22點擊按鈕執(zhí)行cGR曲線完成。如果其它曲線拼接點的位置跟GR相近的話，可以對其它曲線同時拼接，PowerBenchMainWindowEdit→spliceothercurvesspliceothercurves完成。選框，curvesubstitutetable對應位（如圖23R。23如果每條曲線的拼接位置不盡相同的話，那就需要單條拼接，步驟和GR曲線的拼24所示。Page1324項目研究的基本需要，GR1對應著其起始深度，GR2對應著其起始深度，在Assign列輸入CONNCTDeph列輸入GR1ss24項目研究的基本需要，GR1對應著其起始深度，GR2對應著其起始深度，在Assign列輸入CONNCTDeph列輸入GR1ssn列輸入“102”對應Depth列輸入GR2曲線終止深度，GR3對應著其起始深度，拼接情況25所示。251.測井曲線深度校正測井曲線的深度誤差可能因為一系列的引起的，多次下井測量、電纜張力的變化、儀器作非勻速及數(shù)字化的時候人為的錯誤等，深度校正以自然或者深探測電阻率等反應巖性較曲線作為標準，每次測井都會測量，在實際測井的過器由井底向井口進的是真實的深度，越往上測量，電纜的張力拉伸會影響到深度系統(tǒng)的誤差，所以校深順序應該從下往上進行。選擇Edit→depthshift，鼠標即變?yōu)閹в小癝hift”的字樣，點擊圖頭選擇要校深的曲線，26擊按鈕將校深得到的新曲線輸入到數(shù)據(jù)庫。Page1426Edit→ShiftOtherCurvesDepth26Edit→ShiftOtherCurvesDepthShiftTable，對其它曲線進行校深?？?，選擇上面生成的OtherCurves27若一起校深的效果欠佳的話，需要單條校深，校深完后曲線如圖28所示。Page15282.測量環(huán)境對密度、聲波測井質量的影響由于反演的特殊要求，密度和聲波測井曲線會參與到整個反演流度和聲波的質量要格外關注。（1）和低頻模型的質量，是高質量反演的基礎，所以密的可靠性。通?？梢圆捎枚冉粫D-井徑和鉆頭的差值來識別壞井眼數(shù)據(jù)，確定問題Format282.測量環(huán)境對密度、聲波測井質量的影響由于反演的特殊要求，密度和聲波測井曲線會參與到整個反演流度和聲波的質量要格外關注。（1）和低頻模型的質量，是高質量反演的基礎，所以密的可靠性。通?？梢圆捎枚冉粫D-井徑和鉆頭的差值來識別壞井眼數(shù)據(jù)，確定問題Format29所示，“Format”Page1629()Fmtn29()Fmtnetpe可以選擇散點圖ctr)和頻率圖(Cell)，IntervalSelection可以選擇。如果顯示Z軸，則要選擇Z-SegmentZSingleColor即可。(2)Discriminators選項中，可以定義判別式條件來選擇樣本點。(3)Grid選項中，可以自定義網(wǎng)格，默認即可。Matrix選項中，可以定義骨架點和流體的測井響應值、中子測量儀器類型等。Histograms選項中，可以定義直方圖的概率積分曲線和是否顯示直方圖。Print選項支持交會圖的直接打印輸出。Normalization主要是用作曲線的標準化處理，后面會詳細講述首先PowerBenchMainWindow界面選擇Compute→Mathpack，利用Mathpack曲線運算器計算鉆頭曲線S和井徑曲線和鉆頭之間的差值C（如圖30。Page17定義顯示Z軸屬性X、Y、Z軸。30CNL-DEN-CALI交會圖為：31從圖3130CNL-DEN-CALI交會圖為：31從圖31上可以看出，井況不是很好，密度測井響應失真，點擊右鍵選擇startareaselect，利用多邊形將壞數(shù)據(jù)點圈出，然后右鍵選擇Exportvalues/flags，輸出曲線名字定E”，將多邊形內(nèi)的數(shù)據(jù)點命名為“1”，多邊形之外名為“0”，見圖32。Page1832E顯示在曲線模板上，檢查問題把輸出的ND的位置，如圖33。33當井眼垮塌非常普遍的時候，采用多元線性擬合的質類似且相距較近的層段，經(jīng)過對質量較測井曲線進行多元線性擬合，建立基準曲程形式如下：，而后在井眼垮塌層 Logref=f(Rt,GR,…)若井壁坍塌普遍時，則要應用巖石物理模型正演32E顯示在曲線模板上，檢查問題把輸出的ND的位置，如圖33。33當井眼垮塌非常普遍的時候，采用多元線性擬合的質類似且相距較近的層段，經(jīng)過對質量較測井曲線進行多元線性擬合，建立基準曲程形式如下：，而后在井眼垮塌層 Logref=f(Rt,GR,…)若井壁坍塌普遍時，則要應用巖石物理模型正演Com→SnecCreGrtrSttcPage19框中，選擇樣本層與密度曲線相關性比較曲線，樣本層段CurveGenerator選取的原則為（1井眼條件良好（2與待校正層段（以免引入聲波和密度曲線的壓實趨勢（3）與待校正層段處于同一地質層段范圍內(nèi)（確保沉積環(huán)境沒有突變；（4框中，選擇樣本層與密度曲線相關性比較曲線，樣本層段CurveGenerator選取的原則為（1井眼條件良好（2與待校正層段（以免引入聲波和密度曲線的壓實趨勢（3）與待校正層段處于同一地質層段范圍內(nèi)（確保沉積環(huán)境沒有突變；（4與待校正層段的巖性相同或相似（5與待校正層段具有相同或相似的流體類型。34Demo1圖35（如。35一條密度曲線，命名為DENED，然后利用Mathpack模塊對首先利用MathpackPage20曲線井和樣本層段基準曲線異常尖峰逐次的，完的密度曲線如圖36。36（2）對于老點的聲波測井儀器，單發(fā)單收或者單發(fā)雙收，由于線路的現(xiàn)聲波時差曲線的跳變，以前我們把這種縱波跳躍視為氣層的標志，現(xiàn)在聲波儀器相比最初的單發(fā)單收儀器，不僅僅作到井眼補償，而且線路的設計也趨于完善，在遇到含氣地層的時候，首波沒激發(fā)，第二個或者三個波至來的時候拾取時差，這從某種程度上削弱了聲波曲線的跳變。對于常用的雙發(fā)雙收補償聲波測井曲線來說，儀器的是首波信號，測量過由于儀器遇卡、打轉、與地層巖石磕碰等會出現(xiàn)尖峰、錯亂異常響應，這些與地質因素無關的可根據(jù)巖石物理極限是否（其中S代表砂巖，S代表石灰?guī)r，DL代表白云巖異常尖峰逐次的，完的密度曲線如圖36。36（2）對于老點的聲波測井儀器，單發(fā)單收或者單發(fā)雙收，由于線路的現(xiàn)聲波時差曲線的跳變，以前我們把這種縱波跳躍視為氣層的標志，現(xiàn)在聲波儀器相比最初的單發(fā)單收儀器，不僅僅作到井眼補償，而且線路的設計也趨于完善，在遇到含氣地層的時候，首波沒激發(fā)，第二個或者三個波至來的時候拾取時差，這從某種程度上削弱了聲波曲線的跳變。對于常用的雙發(fā)雙收補償聲波測井曲線來說，儀器的是首波信號，測量過由于儀器遇卡、打轉、與地層巖石磕碰等會出現(xiàn)尖峰、錯亂異常響應，這些與地質因素無關的可根據(jù)巖石物理極限是否（其中S代表砂巖，S代表石灰?guī)r，DL代表白云巖，應用巖石物理建模方法予以校正。Demo1井縱波總體質量較好，在可37所示。的巖石物理極限范圍之內(nèi)，如圖Page2137（3）了偶極子聲波測井儀器測量撓曲波來近似的得到橫波時差波，可以采用交叉方式測量，在低頻部分(<2kHz)是以地層橫波速度傳播；高頻時低于橫波速度，必須進行頻散校正。（二）多井一致性處理環(huán)節(jié)。由于不同系列的測井儀器可能的系統(tǒng)誤差、各井使用的泥漿性能的差異、以及井眼等因素的影響，使得不同井之間相同測井曲線在標準層的測井響應 很大的差異。這種差異的和受批量井約束的巖性反演，可能存井、震資料結合的重要質控環(huán)節(jié)。的基礎上，多井資料標準化校正，是多井資料標準化的關鍵是選擇標準層，標準層的選取要滿足如下的條件：（1）（2）（3）全區(qū)都有分布，厚度適中，米左右為宜巖性、物性基本，最好選擇泥巖、膏巖等區(qū)域蓋層37（3）了偶極子聲波測井儀器測量撓曲波來近似的得到橫波時差波，可以采用交叉方式測量，在低頻部分(<2kHz)是以地層橫波速度傳播；高頻時低于橫波速度，必須進行頻散校正。（二）多井一致性處理環(huán)節(jié)。由于不同系列的測井儀器可能的系統(tǒng)誤差、各井使用的泥漿性能的差異、以及井眼等因素的影響，使得不同井之間相同測井曲線在標準層的測井響應 很大的差異。這種差異的和受批量井約束的巖性反演，可能存井、震資料結合的重要質控環(huán)節(jié)。的基礎上，多井資料標準化校正，是多井資料標準化的關鍵是選擇標準層，標準層的選取要滿足如下的條件：（1）（2）（3）全區(qū)都有分布，厚度適中，米左右為宜巖性、物性基本，最好選擇泥巖、膏巖等區(qū)域蓋層井眼條件好，測井資料真實可靠標準化還要選取至少一口井作為標準井，標準井要滿足如下的條件：（1）（2）（3）位于構造的有利部位，鉆井較深，代表性好，直井為宜測井系列齊全，測井條件優(yōu)越，油基泥漿最好有必要的成像測井項目的錄井、取心、測試及分析化驗資料測井曲線標準化主要有以下三種：Page22（1）統(tǒng)計標準化量有相同的平均值和標準偏差。直方圖累計概率曲線匹配（2）在一個大綜合直方圖中。根據(jù)統(tǒng)計的“正態(tài)分布”分析該綜合直方圖，產(chǎn)生一符。（3）趨勢面分析標準化標準層的均值響應作為井點處的測井響應標準，再通過插值算法建立“區(qū)域趨勢面”，各井的標準層均值響應應與趨勢面相同。三種都有各自的適用性范圍，統(tǒng)計標準化要求標準地層地質相態(tài)沒有變只要求標準地層地質相態(tài)不變，而考慮到標準層地質相態(tài)有可能發(fā)生的砂地比變化趨勢和砂地比有可能變化的特點，適用性最強，但操作上也最復雜。選取Demo2井作為標準橫波資料套（1）統(tǒng)計標準化量有相同的平均值和標準偏差。直方圖累計概率曲線匹配（2）在一個大綜合直方圖中。根據(jù)統(tǒng)計的“正態(tài)分布”分析該綜合直方圖，產(chǎn)生一符。（3）趨勢面分析標準化標準層的均值響應作為井點處的測井響應標準，再通過插值算法建立“區(qū)域趨勢面”，各井的標準層均值響應應與趨勢面相同。三種都有各自的適用性范圍，統(tǒng)計標準化要求標準地層地質相態(tài)沒有變只要求標準地層地質相態(tài)不變，而考慮到標準層地質相態(tài)有可能發(fā)生的砂地比變化趨勢和砂地比有可能變化的特點，適用性最強，但操作上也最復雜。選取Demo2井作為標準橫波資料套井標準層為（Normal_T）3083m—（Normal_B）3275m。Page23PowerBenchMainWindow界面選擇View→MultiwellPowerBenchMainWindow界面選擇View→MultiwellHistogram，分別顯示兩口井標準層段的、密度、中子和聲波時差的概型分布，如圖40。GRDENPage24DTCNL40從上面概型分布的直方圖可以看出兩口井的自然不同，但中值相同所以均不用做標準化，而兩口井聲波時差曲線著明顯的偏標準化對聲波曲線的意義，對于中子曲線可以通過對比標準化前后的氣層的位置來曲線進行標準化，在中子曲線直方圖上點擊右鍵，選擇“Format”，彈出圖41。41DTCNL40從上面概型分布的直方圖可以看出兩口井的自然不同，但中值相同所以均不用做標準化，而兩口井聲波時差曲線著明顯的偏標準化對聲波曲線的意義，對于中子曲線可以通過對比標準化前后的氣層的位置來曲線進行標準化，在中子曲線直方圖上點擊右鍵，選擇“Format”，彈出圖41。41（Rf（NmiOK。FormatPage2542框（如圖43。432-Point，1-Point：用一個點來標準化有兩種概率積分曲線的左右移動，即X’=X+a；2-Point：用兩個點來44所示。概率積分曲線的左右移42框（如圖43。432-Point，1-Point：用一個點來標準化有兩種概率積分曲線的左右移動，即X’=X+a；2-Point：用兩個點來44所示。概率積分曲線的左右移XXX’aa10thpercentilebcd50thpercentileX’=X+a90thpercentileX’=(X-c)*(b-a)/(d-c)+a44概型分布基本相同。標準化以后的Demo1井和Demo2Page2645sExportvalues/flags，輸出曲線名字定義為“GAS“0”。45sExportvalues/flags，輸出曲線名字定義為“GAS“0”。名為46E曲線與sw_corGAS一致性，這也說明中子標準化是非常有意義的。Page2747同樣對聲波時差曲線進行標準化，并將標準化前后的壓實趨勢進行對比。（1）標準化前（2）標準化后4847同樣對聲波時差曲線進行標準化，并將標準化前后的壓實趨勢進行對比。（1）標準化前（2）標準化后48由圖48可以看出，標準化以前，壓實趨勢是一致的，標準化以后，壓實趨勢發(fā)生Page28了明顯的偏移，結合測井曲線分析認為Demo1聲波時差較Demo2井較小可能是由于Demo1Demo1Demo2井泥質偏輕所導致。交會圖技術進行多數(shù)據(jù)的標準化質量檢,布規(guī)律是否協(xié)調一致。Chapter4. 測井資料地層評價（一）泥質含量的計算泥質含量的計算有很多種，自然、自然電位、電阻率、中子及度交會等，但是每一種都有其應用的有利和不利條件，自然測井計算泥質含量的不利度地層和水層。常用的為交會圖分析（度，中子-聲波，聲波-密度交會圖，度交會法的粘土含量參數(shù)，用來刻度應用自然了明顯的偏移，結合測井曲線分析認為Demo1聲波時差較Demo2井較小可能是由于Demo1Demo1Demo2井泥質偏輕所導致。交會圖技術進行多數(shù)據(jù)的標準化質量檢,布規(guī)律是否協(xié)調一致。Chapter4. 測井資料地層評價（一）泥質含量的計算泥質含量的計算有很多種，自然、自然電位、電阻率、中子及度交會等，但是每一種都有其應用的有利和不利條件，自然測井計算泥質含量的不利度地層和水層。常用的為交會圖分析（度，中子-聲波，聲波-密度交會圖，度交會法的粘土含量參數(shù)，用來刻度應用自然計算的粘土含量。1、度交會確定干、濕粘土點的測井響應參數(shù)View→Crossplots→Neutron-Density，生成度交會圖，右鍵Format，彈出Format字、字母或者“*”字符代表數(shù)據(jù)點出現(xiàn)的次數(shù)，1~9數(shù)字表示，9~35用A~Z表示，大于35的用*表示。利用黃色三角形確定濕粘土點的測井響應，紅色三角形用來確定干粘土點49各點代表的含義為：BCADE49Page29砂巖骨架點（A點：砂巖的骨架礦物主要是石英，其次是長石，此外還有%的重礦微斜長石，而正長微斜長石英的2.57g/cc，中子測井值也相近，故不再區(qū)分長石，常以石英代表砂巖的骨架礦物，當按砂巖刻度的時候，石英骨架點的密度為2.65g/cc，當儀器以純石灰?guī)r為標準進行刻度時，砂巖骨架含氫指數(shù)為負值，約-0.01~-0.05。（B點SS、LS、DOL：依次為飽含孔隙度和密度孔隙度都為100%的點。砂巖線、灰?guī)r線及白云巖線，橫坐標中子為石灰?guī)r刻度的視石灰?guī)r孔隙度，所以只有石灰?guī)r線是線性變化的，其它兩條巖性線都略有彎曲。束縛（C點水，水與干粘土點的連線即為含粘土線。干粘土點（E點砂巖骨架點（A點：砂巖的骨架礦物主要是石英，其次是長石，此外還有%的重礦微斜長石，而正長微斜長石英的2.57g/cc，中子測井值也相近，故不再區(qū)分長石，常以石英代表砂巖的骨架礦物，當按砂巖刻度的時候，石英骨架點的密度為2.65g/cc，當儀器以純石灰?guī)r為標準進行刻度時，砂巖骨架含氫指數(shù)為負值，約-0.01~-0.05。（B點SS、LS、DOL：依次為飽含孔隙度和密度孔隙度都為100%的點。砂巖線、灰?guī)r線及白云巖線，橫坐標中子為石灰?guī)r刻度的視石灰?guī)r孔隙度，所以只有石灰?guī)r線是線性變化的，其它兩條巖性線都略有彎曲。束縛（C點水，水與干粘土點的連線即為含粘土線。干粘土點（E點：干粘土點就是含水濕粘土點的骨架點，只含結晶水，不含束縛水，。（D點粘土線上，粘土的含水量=（干粘土點到濕粘土點的距離）/（干粘土點到的距離。從2.65g/cc；度交會圖上可以看出砂巖骨架點度測井響應分別為：-0.02和度測井響應分別為：1.0和1.0g/cc。泥巖由50~70%的粘土4%70%作為濕粘土點。5%由長 碳酸鹽組成的其它礦物，紅色三角形束縛度測井響應分別為：1.1和1.2g/cc。依據(jù)地區(qū)經(jīng)驗可以適當調整，砂巖骨架點不變，在含 粘土線上，根據(jù)地層束縛水飽和度的大小來確定干粘土點，我們分析確定干粘土點的0.362.76g/cc。PowerBenchMainWindow界面選擇Interpret→CalyVolume→Crossplot，選用度交VCLND。50Page302GR曲線計算泥質含量是由泥質吸附的放射性元素決定的，因此常用自然測確定巖石的泥質含量。利用GRGR（GRH和純巖石的GR2GR曲線計算泥質含量是由泥質吸附的放射性元素決定的，因此常用自然測確定巖石的泥質含量。利用GRGR（GRH和純巖石的GR（GRS。LogplotEdit→Trend/SquareGRSHGRSSGRSSGRSH的GRSH線的偏移。51GRSSGRSHInterpret→ClayVolumeClayVolume框，CurvesGRSSGRSHCleanClay欄。GR計算泥質含量有不同的算法，常用的5種，具體選用哪種，要看計算結果那個與度GRSH。，需要重新回去修訂GRSH和GRSS，尤其52Page3153CurveGenerator計算與法來求取。若是有幾種粘土含量曲線吻合最54的異常尖峰可以通過Mathpack53CurveGenerator計算與法來求取。若是有幾種粘土含量曲線吻合最54的異常尖峰可以通過MathpackVCLAY=0.7+（VCLGR-0.7）*0.3，局部如果效果欠佳的話，可以公式的相應系數(shù)，直至達到理想效果。圖55和圖56分Demo2井求取的泥質含量。Page32PagePage33（二）孔隙度的計算孔隙度的計算也有很多種，可以選擇聲波、密度孔隙度曲線中的一條來計算孔隙度，也可以選擇任意兩條孔隙度曲線的交會分析來計算孔隙度。1的總孔隙度。地方，反映準確中子孔隙度反映的是巖石的總孔隙度，它是巖石含氫指數(shù)的反映，除了孔隙空間中流體的含氫指數(shù)，還粘土結晶水、束縛水的含氫指數(shù)及巖石骨架的等效含氫指可以得到比較準確的總孔隙度。的影響。認為反映原生孔隙度，不反映次生孔隙度，不受 和高角度裂縫2．孔隙度的計算及泥質校正M計算密度視孔隙（二）孔隙度的計算孔隙度的計算也有很多種，可以選擇聲波、密度孔隙度曲線中的一條來計算孔隙度，也可以選擇任意兩條孔隙度曲線的交會分析來計算孔隙度。1的總孔隙度。地方，反映準確中子孔隙度反映的是巖石的總孔隙度，它是巖石含氫指數(shù)的反映，除了孔隙空間中流體的含氫指數(shù)，還粘土結晶水、束縛水的含氫指數(shù)及巖石骨架的等效含氫指可以得到比較準確的總孔隙度。的影響。認為反映原生孔隙度，不反映次生孔隙度，不受 和高角度裂縫2．孔隙度的計算及泥質校正M計算密度視孔隙PHIDPHIDT，計算的沒有經(jīng)過泥質校正和油氣校正的聲波孔隙度和密度孔隙度相差會很大。（1）密度孔隙度泥質校正572.76是干粘土對應的密度響應值。（2）泥質的分布形式對聲波孔隙度影響較大，其中結構泥質和層狀泥質可以直接根據(jù)加權計算，同度一樣，而分散泥質則占據(jù)了部分孔隙，直接減掉就可以。泥質分布泥質很少。Page3458DT_VCL-56)/(189-56)，由于我們評價干粘土點的測井響應參數(shù)的時候，只計算了密度測井響應，這樣我們通過調節(jié)DT_VCL58DT_VCL-56)/(189-56)，由于我們評價干粘土點的測井響應參數(shù)的時候，只計算了密度測井響應，這樣我們通過調節(jié)DT_VCL，直到達到理想的效果。當DT_VCL=90us/ft的時候，兩者吻合較好。所以干粘90us/ft是干粘土對應的聲波時差響應值。（3）選擇Compute→NeutronMatrixConversion，將石灰?guī)r骨架下的中子值轉化到相應骨架下的中子測井值，要選擇儀器的類型和廠商，這是我們測井地層評價必須的資料。59Page35中子孔隙度泥質校正公式：PHINC=PHIN-VCLAY*(-0.02-0.36)/(-0.02-1)，注明：0.36是干粘土對應的中子響應值。（4）度孔隙度泥質校正lotPorosity，定義輸出曲線的名字為PHIND，輸入標準化后的中子和密度以及儀器類型和廠商。60中子孔隙度泥質校正公式：PHINC=PHIN-VCLAY*(-0.02-0.36)/(-0.02-1)，注明：0.36是干粘土對應的中子響應值。（4）度孔隙度泥質校正lotPorosity，定義輸出曲線的名字為PHIND，輸入標準化后的中子和密度以及儀器類型和廠商。60度交會圖板上，由度的測井響應值可以計算得到干粘土對應的18個孔隙度，如圖61所示。61Page36（5）將por_cor投放到繪圖模板，觀察其與上述計算的四種形式的孔隙度是否一致。62Demo1（5）將por_cor投放到繪圖模板，觀察其與上述計算的四種形式的孔隙度是否一致。62Demo1Generator模塊將計算的孔隙度與巖心孔隙度進行相關63。633Compute→MathpackPHIT=PHIDC*0.35+PHINDC*0.35+0.2*PHINC+PHIDTC*0.164所示。Page37Demo1Demo2644在上面孔隙度和泥質含量計算的基礎上，利用MathPackDemo1Demo2644在上面孔隙度和泥質含量計算的基礎上，利用MathPack計算得到砂巖含量VQUA=1-PHIT-VCLAY，兩口井的體積模型見圖65。Page38Demo1Demo265Page39Demo1Demo265Page39（三）飽和度的計算含水飽和度的計算力計算初始含水飽和度等。也有很多種，如公式是一種純巖石的解釋理過。但是由于泥質巖石導電性的復雜性和地域差異性，使得響應復雜多樣，難說哪個響應選用和調整。西門杜表明，巖石的導電性可看成是泥質巖石的適用于地層水礦化度較高的地層。泥質部分和純巖石部分并聯(lián)導電，這種亞表明泥質巖石的導電性是三部分組成：地層水網(wǎng)絡的導電性、粘土網(wǎng)絡的導電性及兩個網(wǎng)絡交叉連接引起的導電性，適用于地層水礦化度較低的地層。W-S模型結為粘土的陽離子交換作用的結果。這兩種模型除了電導率響應不同外，其它測井W-S模型認為泥質砂巖的導電性是砂巖骨架礦物和孔隙空間電解液和粘土陽離子交換并聯(lián)導電的結果，考慮到粘土礦物的附加導電性，但是利用測井資料確定Qv有相當大的難度。雙水模型是含泥質儲集層的一種解釋模型，它把巖石固體干粘土）都看成是不導電的，而把地層水分為兩部分：束縛水。這兩部分水的導電性是不同的，它不需要用巖心測量確定Qv，可廣泛地應用于實際測量解釋。根據(jù)律服從單元后主要的聚集動力就理論。通常認為地層中水界面（WL，在該界面以下地層中充滿理論計算含氣飽和度，引入了描述巖石流體物100%。為了采用理性質的無量綱J函數(shù)，通過J函數(shù)可以把所有的數(shù)據(jù)統(tǒng)一到一條線上，J函數(shù)在對數(shù)坐標下具有線性，J（Sw）=A/SwB，這樣我們就能計算出含水飽和度。1Pickett（三）飽和度的計算含水飽和度的計算力計算初始含水飽和度等。也有很多種，如公式是一種純巖石的解釋理過。但是由于泥質巖石導電性的復雜性和地域差異性，使得響應復雜多樣，難說哪個響應選用和調整。西門杜表明，巖石的導電性可看成是泥質巖石的適用于地層水礦化度較高的地層。泥質部分和純巖石部分并聯(lián)導電，這種亞表明泥質巖石的導電性是三部分組成：地層水網(wǎng)絡的導電性、粘土網(wǎng)絡的導電性及兩個網(wǎng)絡交叉連接引起的導電性，適用于地層水礦化度較低的地層。W-S模型結為粘土的陽離子交換作用的結果。這兩種模型除了電導率響應不同外，其它測井W-S模型認為泥質砂巖的導電性是砂巖骨架礦物和孔隙空間電解液和粘土陽離子交換并聯(lián)導電的結果，考慮到粘土礦物的附加導電性，但是利用測井資料確定Qv有相當大的難度。雙水模型是含泥質儲集層的一種解釋模型，它把巖石固體干粘土）都看成是不導電的，而把地層水分為兩部分：束縛水。這兩部分水的導電性是不同的，它不需要用巖心測量確定Qv，可廣泛地應用于實際測量解釋。根據(jù)律服從單元后主要的聚集動力就理論。通常認為地層中水界面（WL，在該界面以下地層中充滿理論計算含氣飽和度，引入了描述巖石流體物100%。為了采用理性質的無量綱J函數(shù)，通過J函數(shù)可以把所有的數(shù)據(jù)統(tǒng)一到一條線上，J函數(shù)在對數(shù)坐標下具有線性，J（Sw）=A/SwB，這樣我們就能計算出含水飽和度。1Pickett交會確定地層水電阻率Pickett交會圖可以認為是公式F=Ro/Rw=a/?mlgRt=-mlg?+lg(abRw/Swn)所以在雙對數(shù)坐標下，Rt和?之間的是一組斜率為-m，截距為lg(abRw/Swn)的，直線的截距僅隨Sw而變，這樣即為直線。當巖性Rw不變的Pickett交會圖中一組隨Sw變化的平行直線。由于Demo2水井，所以選擇Demo2PowerBenchGRCutoff，也可以通過VCLAY來取Cutoff，將泥巖點剔出。Page4066點擊右鍵，選擇SwLines-1PtFixedM，計算得到Rw=0.073，a=1，66點擊右鍵，選擇SwLines-1PtFixedM，計算得到Rw=0.073，a=1，m=2，n=2。2利用上步計算得到的結果進行下面含水飽和度的計算，PowerBenchMainWindow界面選擇Interpret→WaterSaturation→Archie，Simandoux則相反。亞比較適合地層水比較咸的地層，而67Page41（四）ShalySandCP-9測井解釋（四）ShalySandCP-9測井解釋ShalySand模塊主要已知相同巖性，泥質含量不同的儲層，輸入骨架類型（砂C-9圖表計算視油poupon-gaymard（rPage42經(jīng)驗公式來計算剩余油飽和度。Sand，ShalySandCP-969所示。經(jīng)驗公式來計算剩余油飽和度。Sand，ShalySandCP-969所示。69PhiNeutron：相應骨架的中子孔隙度NeutronType：儀器類型PhiDensity：相應骨架下的密度孔隙度VCL：粘土含量，前面已計算出Temp(degF)：溫度曲線或者定值PhiMax：孔隙度的最大值RxoClay：泥巖層沖洗帶電阻率RtClay：泥巖層真電阻率PhiNeutronClay：泥巖層中子孔隙度PhiDensityClay：泥巖層密度孔隙度CrushPorosityWhenVCL>VCLCutoffCutoff的時候，孔隙度指數(shù)減小，BoostmWhenVCL>：適用于泥質偏重的地層，用來增加含水飽和度BoostSwWhenPorosity<：適用于低孔地層，用來增加含水飽和度Sw/SxoExponent：適用于利用Rxo計算沖洗帶含水飽和度的情況，用來約束Sxo，默認0.2InvasionFactor：適用于利用經(jīng)驗公式計算沖洗帶含水飽和度的情況Rt：深電阻率Rw：地層水電阻率Rxo：微電阻率曲線Rmf：泥漿濾液電阻率a：公式巖性系數(shù)Page43m：n：含水飽和度指數(shù)（五）Quicklook快速測井解釋Quicklook計算孔隙度和含水飽和度，模型中所涉及參數(shù)含義與Shalym：n：含水飽和度指數(shù)（五）Quicklook快速測井解釋Quicklook計算孔隙度和含水飽和度，模型中所涉及參數(shù)含義與ShalySand基本相同，在此不加以贅述。70Page44（六）曲線統(tǒng)計報告和總結報告1. 曲線統(tǒng)計報告曲線統(tǒng)計報告可以提供多井多層段的曲線的統(tǒng)計，例如曲線的最大值和最小值、均值、等（六）曲線統(tǒng)計報告和總結報告1. 曲線統(tǒng)計報告曲線統(tǒng)計報告可以提供多井多層段的曲線的統(tǒng)計，例如曲線的最大值和最小值、均值、等Page4572AddtoSelections可以人為加載曲線；選擇MergefromSelectedWells72AddtoSelections可以人為加載曲線；選擇MergefromSelectedWells一次加載全部曲線，或者利用CurveFilter有選擇性地加載曲線。在ReportOptions選項中，在CurveFields選擇曲線的輸出，例如我們選擇輸出上面選擇曲線的最大值、最小值、測井曲線名稱和測井曲線的ReportFileName一欄，可以選擇您要文件的地址，設置完成后，點擊RUN，生成的Report如下。732. 總結報告和敏感性報告在設定的CutoffSummationReport模塊可以生成多井多層段的ZonesFile→MakeZones，彈出MakeZonesToolZoneZone的名稱，Top一欄選擇層段的頂界深度，Bottom一欄選擇層段的底界深度，點擊“OK”即可。Page4674根據(jù)分層數(shù)據(jù)我們設定了三個Zones，生成總結性報告。75在C74根據(jù)分層數(shù)據(jù)我們設定了三個Zones，生成總結性報告。75在CfZsmmnRrWQ.PHQ1；Options選擇要輸出的曲線和的路徑。生成的報告如下：Page4776敏感性報告：設定臨界值的范圍（開始值、結束值和步長值，W（0，，76敏感性報告：設定臨界值的范圍（開始值、結束值和步長值，W（0，，.；HI（0.05VC（01.組合條件下的曲線的數(shù)值。77生成的敏感性報告如下：78Page48Chapter5. 最優(yōu)化測井解釋傳統(tǒng)的測井解釋和解釋程序（POR程序、CRA程序等）基本上都是應用有限的幾種測井資料進行地層參數(shù)的評價和油氣分析評價。與常規(guī)測井解釋思路不同，StatMin通過建立和求解物理意義明確的線性線來得到地層礦物組分、孔隙度等巖性和物性組，輸入任意反映“四性”特征的測井曲。物組分的極大似然評價，StatMin可以計算出與曲線匹配最佳的結果。對于欠定Mn提供了下列技術來解決這個問題（1（優(yōu)先矩陣（2輸入外部還將輸出所有基于模型結果的重構曲線和DELTAChapter5. 最優(yōu)化測井解釋傳統(tǒng)的測井解釋和解釋程序（POR程序、CRA程序等）基本上都是應用有限的幾種測井資料進行地層參數(shù)的評價和油氣分析評價。與常規(guī)測井解釋思路不同，StatMin通過建立和求解物理意義明確的線性線來得到地層礦物組分、孔隙度等巖性和物性組，輸入任意反映“四性”特征的測井曲。物組分的極大似然評價，StatMin可以計算出與曲線匹配最佳的結果。對于欠定Mn提供了下列技術來解決這個問題（1（優(yōu)先矩陣（2輸入外部還將輸出所有基于模型結果的重構曲線和DELTA曲線，作為模型質量的。最優(yōu)化測井解釋數(shù)據(jù)準備和數(shù)據(jù)預處理的過程，要注意以下幾個方面問題：（1）StatMin利用線性組求解地層的，要求輸入為線性曲線。 地，常用的數(shù)是非線性的，需要在最優(yōu)化求解之前，進行線性化的預處理。①自然：曲線主要用來求解泥質含量和粘土含量參數(shù)，而自然和泥質含量或粘土含量間的響應StatMin最優(yōu)化測井模型前需要選擇模型計算得到泥質含量Vsh或粘土含量Vcl。②中子孔隙度：灰?guī)r地層來說，中子孔隙度與地層物性參數(shù)是線性的響應黃色三角和藍色三角分別為石英砂巖和云巖的線性化骨架點參數(shù)。79③電阻率：由地層因素公式可知，電導率與孔隙度：CtCw*?ma；其中，Ct：巖石的電導率；Cw：地層水的電導率；m：巖石膠結指數(shù)；a：巖性系數(shù)，當m、a值已知時（可根據(jù)巖電分析測試獲得，可將Ct的m層物性參數(shù)間的響應是線性的，再參與到最優(yōu)化求解當中。另外一個就是根據(jù)原狀地層束縛水體積BVW以及沖洗帶束縛水體積BVWSXO等參數(shù)，這些體積組分參數(shù)Page49與物性參數(shù)間的響應顯然是線性的。等公式計算這些參數(shù)的時候會用到孔隙度參數(shù)，所以這是一個迭代工作流程。U與地層物性參數(shù)間的響應可以近似看做是線性的，U的計算表是：U=Pef*ρe=Pef*(ρb +0.1883)/1.0704有可能是線性的，例如時間平均公式。但是，依據(jù)地質條件的差異，它們之間的也能滿足H-R-G模型，是非線性的。總之，在最優(yōu)化求解之前，有必要考查聲波同物性參數(shù)之間的 式，按照適當?shù)倪M行線性化處理。76514V3TimeaverageWood與物性參數(shù)間的響應顯然是線性的。等公式計算這些參數(shù)的時候會用到孔隙度參數(shù)，所以這是一個迭代工作流程。U與地層物性參數(shù)間的響應可以近似看做是線性的，U的計算表是：U=Pef*ρe=Pef*(ρb +0.1883)/1.0704有可能是線性的，例如時間平均公式。但是，依據(jù)地質條件的差異，它們之間的也能滿足H-R-G模型，是非線性的?？傊?，在最優(yōu)化求解之前，有必要考查聲波同物性參數(shù)之間的 式，按照適當?shù)倪M行線性化處理。76514V3TimeaverageWood22V2Reuss1MReussVReuss000.81Porosity80（2）聲波-密度交會圖上的響應范圍，可以發(fā)現(xiàn)70%的測井數(shù)據(jù)點落在這樣的區(qū)間內(nèi)，膏巖、等都可以應用類似確定。81而某些類型的巖石，礦物組成成分多樣，例如泥質砂巖，可以看成是純泥巖與純砂Page50Vp(km/s)1初步定的測井響應誤差范圍參數(shù)。地質資料檢驗來逐步獲得穩(wěn)1典型最優(yōu)化測井解釋測井儀器測井響應誤差范圍（3）依據(jù)最優(yōu)化測井解釋原理，測井響應 可以表達成以下形式：[EP] [V]=[L]其中，EP=各體積模型組分的測井響應骨架點參數(shù)；L=實際測井曲線響應；按照矩陣的運算法則，待求解量V可表示成以下形式，[EP]–1[L]=[V]各符號意義同上：也就是說，骨架點參數(shù)矩陣的逆矩陣著整個矩陣計算過程的陣中出現(xiàn)量級差別很大的數(shù)值時，矩陣計算變得不（1骨架端點值矩陣的行不能[0001]、[0010]、[0100]和[1000]1初步定的測井響應誤差范圍參數(shù)。地質資料檢驗來逐步獲得穩(wěn)1典型最優(yōu)化測井解釋測井儀器測井響應誤差范圍（3）依據(jù)最優(yōu)化測井解釋原理，測井響應 可以表達成以下形式：[EP] [V]=[L]其中，EP=各體積模型組分的測井響應骨架點參數(shù)；L=實際測井曲線響應；按照矩陣的運算法則，待求解量V可表示成以下形式，[EP]–1[L]=[V]各符號意義同上：也就是說，骨架點參數(shù)矩陣的逆矩陣著整個矩陣計算過程的陣中出現(xiàn)量級差別很大的數(shù)值時，矩陣計算變得不（1骨架端點值矩陣的行不能[0001]、[0010]、[0100]和[1000]，這主要GR和P（2）輸入骨架端點值矩陣的兩避免的。情況。這是應該注意和資料比對，計算的孔隙度與巖心分析孔隙度對比檢驗，以考察解釋結果的準確性。下面Demo1StatMin建立其體積模型首先StatMin中提供了中子和密度線性化處理模塊，Interpret→Statmin→Computes，彈出圖82。Page51密度（g/cc）中子（decimal）縱波時差（us/ft）（g/cc*ba/e）自然（API）地層水0.01350.35天然氣0.01150.35砂巖0.02130.35石灰?guī)r0.010.530.55泥頁巖0.03330.55白云巖0.02230.35硬石膏0.02111.55石膏0.02130.35鹽巖0.03110.350.0310.30.3582選82選1誤差如下：83層段，調整骨架值和標準誤差，使得原曲線和重構曲線盡量接近，新生成的吻合不Page52DELTA盡量小。84Chapter6.巖石物理建模與分析巖石物理模型分析工作主要使用RPMDELTA盡量小。84Chapter6.巖石物理建模與分析巖石物理模型分析工作主要使用RPM巖石物理分析模塊進行。首先把多井地層評價處理得到的各礦物含量、孔隙度、飽和度等參數(shù)輸入RPM模塊，然后選擇合理的巖石物理模型進行波速度的建模，通過微調骨架點參數(shù)，使模型數(shù)據(jù)和實測數(shù)據(jù)達到很高的相關性，然后確定可用于工區(qū)內(nèi)的巖石物理模型和骨架參數(shù)點。巖石物理建模Page53可以進行密度、縱波和橫波的正演工作，為疊前反演提供可靠的橫波資料，同時補齊了一些井缺失的密度曲線，校正了一些井由于井眼垮塌造成的密度曲線失真。（一）巖石物理基礎1.巖石物理理論模型利用理論來礦物組分的體積含量和各礦物的彈性模量特征那么只能 出彈性模量的上下（圖85要精確具體值，還需了解其孔隙的幾何特征。圖85彈性模量上下限值示意圖①Voigt、Reuss極限模型成分假設有相等應變時，應力與平均應變的比值；Reuss界限被稱之為等應力平均，它給出了當成分有相等的應力時平均應力與平均應變的比值。Reuss平均精確地描述了在孔隙流體中懸浮顆粒情況下的等效模量。1Voigt_M Mi可以進行密度、縱波和橫波的正演工作，為疊前反演提供可靠的橫波資料，同時補齊了一些井缺失的密度曲線，校正了一些井由于井眼垮塌造成的密度曲線失真。（一）巖石物理基礎1.巖石物理理論模型利用理論來礦物組分的體積含量和各礦物的彈性模量特征那么只能 出彈性模量的上下（圖85要精確具體值，還需了解其孔隙的幾何特征。圖85彈性模量上下限值示意圖①Voigt、Reuss極限模型成分假設有相等應變時，應力與平均應變的比值；Reuss界限被稱之為等應力平均，它給出了當成分有相等的應力時平均應力與平均應變的比值。Reuss平均精確地描述了在孔隙流體中懸浮顆粒情況下的等效模量。1Voigt_M MiReuss_Mi Miiii86VoigtReuss界限示意圖②Voigt-Reuss-Hill平均模量模型Voigt-Reuss-Hill平均是將Voigt上限和Reuss下限取算數(shù)平均來估算巖石的等效彈性模量:MMi i i1iiiVRH_M2③Hashin-Shtrikman界限模型Page54通過交換巖石兩種成分材料計算的Hs-rkmn（圖87當巖石的孔隙幾何細節(jié)未知情況下最窄的可容許上下界限。fKK KKfK43f212fK21通過交換巖石兩種成分材料計算的Hs-rkmn（圖87當巖石的孔隙幾何細節(jié)未知情況下最窄的可容許上下界限。fKK KKfK43f212fK211 1 12 15K431 1187④微分有效模量模型（DEM）Hashin-Shtrikman界限示意圖微分有效模量模型（DEM）通過向固體礦物相加入橢圓形的包含物相（孔隙）來模擬雙相混合物。構建微分有效模量巖石物理模型（圖88）主要以下關鍵步驟：a)選取合適的巖石物理理論模型(Xu&White模型等)開始構建整體流程。b)根據(jù)骨架礦物和充填礦物（泥質）的骨架點 （例如Willy時間平均公式、VRH）計算混合骨架的彈性屬性。合巖石骨架再次進行混合計算，得到干巖的彈性屬性。d）采用Gassmann特征。圖88微分有效介質理論巖石物理建模技術流程圖在建立了巖石物理模型后，可以非常容易地用于建模井中不的彈性曲線，該巖石物理模型也可以進一步開展流體替換研究和侵入校正工作。2.巖石物理經(jīng)驗模型巖石的彈性特征受到多種因素的影響，如果根據(jù)實驗測量，能綜合這些因素到一個巖石物理模型當中去，理論上是可以素相互之間的影響機理尚不清楚，難以建立一個這些特征的。但是模型過于復雜，各因的經(jīng)驗主要考慮影響巖石彈性特征的主要因素，例如巖石的類型、粘土含量、孔隙度等。Page55泥巖骨架Vp,Vs,ρ骨架正演成果流體屬性計算 指數(shù)水Vp,ρ K,ρ Gassmann流 K,ρ 流體置換模型 Vs體溫度、 、油密度、地層水礦化度、氣油比，氣比重 VpBatzle&Wang等K,ρ氣Vpρ油VpρK,ρRhoKμK,,骨架混合Vp,Vs,ρK,μ,ρK,,微分有效介質(DEM)計算模型K,μ,ρ灰?guī)r骨架Vp,Vs,ρK,μ,ρ①速度-密度模型最著名的，也是應用最廣泛的要數(shù)Gardner速度-密度式，在有標定模型的情況下可以取得不錯的效果。模型，如下式所示形 aVb dVbbV c2pfppρb為巖石體積密度；Vp為巖石縱波速度；其它為擬合系數(shù)。②速度-孔隙度模型應用較廣泛的速度-孔隙度模型要數(shù)Wyllie和R-H-G式，Wyllie將孔隙度和速度近似成下述：1 1VpVfl Vma其中Vp，Vma，Vfl分別是巖石，巖石骨架礦物和孔隙流體的縱波速度。R-H-G式是的繼承和發(fā)展，其形式如下：WyllieV(1)VV37%①速度-密度模型最著名的，也是應用最廣泛的要數(shù)Gardner速度-密度式，在有標定模型的情況下可以取得不錯的效果。模型，如下式所示形 aVb dVbbV c2pfppρb為巖石體積密度；Vp為巖石縱波速度；其它為擬合系數(shù)。②速度-孔隙度模型應用較廣泛的速度-孔隙度模型要數(shù)Wyllie和R-H-G式，Wyllie將孔隙度和速度近似成下述：1 1VpVfl Vma其中Vp，Vma，Vfl分別是巖石，巖石骨架礦物和孔隙流體的縱波速度。R-H-G式是的繼承和發(fā)展，其形式如下：WyllieV(1)VV37%1 147%VV V10.4710.3713747%V0.1 V0.1 V③速度-孔隙度-粘土含量模型常被的速度-孔隙度-粘土含量泥質砂巖，Castagna根據(jù)測量得到的水飽和狀態(tài)下的速度、孔隙度和粘土含量，回歸得到經(jīng)驗公式：-2.21CVs(km.s)=3.89-7.07-2.04C其中，?是孔隙度，C是粘土含量。④縱波-橫波速度模型常用的縱波-橫波速度 有HanCastagna等人提出的泥質砂巖Vp-VsPcktCastagnaVp-Vs式(89)，這些式只適合于實驗研究的那些樣本巖石所在地域，實際項目中可借鑒參考的形式，根據(jù)實測資料優(yōu)化調整 的有關系數(shù)，才能保證取得較精度。89Vp-Vs式Page56復雜模型。具體選擇哪一種模型，取決于研究的目的以及資料品質和模型適用性。下面通過巖石物理經(jīng)驗模型來熟悉RM的簡單運用owrB復雜模型。具體選擇哪一種模型，取決于研究的目的以及資料品質和模型適用性。下面通過巖石物理經(jīng)驗模型來熟悉RM的簡單運用owrBchMinWw界面選會跳出一個提示錯誤的框，直接點擊“OK”即可。90選擇Data→SelectWells，彈出圖91，選中Demo2井，點擊“OK”91“EntireWell”，點擊“OK”，在空白工區(qū)內(nèi)，點擊鼠標右鍵，選擇“AddNode”，Node，雙擊“n0”，彈出圖92。92選擇相應的公Page57式模塊，圖92要實現(xiàn)的是縱波聲波時差到縱波速度的轉換，DTCO（ft/s）=1000000/DTCO（s，在Prm1式模塊，圖92要實現(xiàn)的是縱波聲波時差到縱波速度的轉換，DTCO（ft/s）=1000000/DTCO（s，在Prm1“0”，在rm2Cre→“DTCO”。Output可以類型選擇“P_velocity”，設置完畢后，點擊“OK”，選擇Calculate→CalculateWorkflow，或者點擊快捷鍵作，這樣就實現(xiàn)了縱波聲波時差到縱波速度的轉換。，執(zhí)行計算操下面我們利用計算的縱波速度來計算橫波速度，分別采用MudRock和GreenbergCastagna兩種實現(xiàn)，點擊鼠標右鍵，再添加n1和n2，定義如圖93所示。93兩種求取的是橫波速度，Demo2檢驗經(jīng)驗模型的實用性，所以要講橫波速度轉變?yōu)闄M波時差。對兩種94所示。94Page5895設置完成以后，流程如圖95所示，MudRock計算的橫波時差曲線為DTSM_M，Greenberg計算得到的橫波時差曲線為DTSM_C。將兩種95設置完成以后，流程如圖95所示，MudRock計算的橫波時差曲線為DTSM_M，Greenberg計算得到的橫波時差曲線為DTSM_C。將兩種與實測的橫波曲線進行Castagna對比發(fā)現(xiàn)，吻合度較高。96Page59（二）儲層流體置換1.地層流體彈性屬性計算油藏狀態(tài)下的單相流體屬性。首先選擇PowerBenchMainWindow界面Compute→BasicsLogFunctions→General和。圖97間的轉換，比如溫度攝氏度要轉化為華氏度，計算出來的溫度和都要進行同時輸入礦化度、氣油比、氣比重等和輸出曲線分別為：表2，流程圖見圖98，輸入Page60輸入曲線輸出曲線（二）儲層流體置換1.地層流體彈性屬性計算油藏狀態(tài)下的單相流體屬性。首先選擇PowerBenchMainWindow界面Compute→BasicsLogFunctions→General和。圖97間的轉換，比如溫度攝氏度要轉化為華氏度，計算出來的溫度和都要進行同時輸入礦化度、氣油比、氣比重等和輸出曲線分別為：表2，流程圖見圖98，輸入Page60輸入曲線輸出曲線n0Pressuren13地層水的密度Rho-brinen1溫度Temperaturen14地層水的體積模量K-brinen30礦化度Salin15天然氣的密度Rho-Gasn31油比重APIn16天然氣的體積模量K-gasn32氣比重GasSpecificGravityn17油的密度Rho-Oiln33氣油比GORn18油的體積模量K-Oil98計算油、氣和水的體積模量和密度界面輸入?yún)?shù)如圖99所示：Page6198計算油、氣和水的體積模量和密度界面輸入?yún)?shù)如圖99所示：Page6199RPM里流體的溫度和圍分別是：1）油：溫度上限：100攝氏度；都有一定的的范圍，其有效范上限：62兆帕。2）氣：溫度上限：150攝氏度；*P(.2521)，99RPM里流體的溫度和圍分別是：1）油：溫度上限：100攝氏度；都有一定的的范圍，其有效范上限：62兆帕。2）氣：溫度上限：150攝氏度；*P(.2521)，P上限：110兆帕。3）(°C)=180.12為兆帕，上限：140兆帕。計算就可得到各單相流體的流體屬性如圖10012曲線；第3道紅色實線為天然氣的體積模量，藍色實線為地層水體積模量；第道為4道紅色實線為天然氣的密度，藍色實線為地層水的密度。Page62100曼儲層流體置換2.對于兩多相混合流體彈性屬性的計算多用wood公式或Brie1995年發(fā)明的來計算，形式如下所示，BrieK KK SKemixbrinegaswgas在e指數(shù)較大的情況下相當于100曼儲層流體置換2.對于兩多相混合流體彈性屬性的計算多用wood公式或Brie1995年發(fā)明的來計算，形式如下所示，BrieK KK SKemixbrinegaswgas在e指數(shù)較大的情況下相當于wood公式描述的流體飽和度與體積模量e指數(shù)情況體含水飽和度與體積模量差異很大（圖101行標定。地，e指數(shù)在3附近，但嚴格數(shù)值仍需要根據(jù)實測資料進Brie中流體飽和度與體積模量Page63巖石物理分析的重要任務之一就是從一種流體飽和的巖石速度分析中另的彈性模量特征，速度，等同地，可以說根據(jù)干巖骨架的彈性模量和孔隙流體巖石的速度，這就是孔隙流體置換問題。Biot–Gassmann是大多數(shù)完善的巖石物理模型所采用的流體置換102所示。102輸入和輸出曲線分別為：Page64輸入曲線輸出曲線n1Pressuren36流體置換后的密度Densityn2溫度Temperaturen37流體置換后的縱波速度Vpn3含水飽和度WaterSaturationn38流體置換后的縱波時差Psonicn4有效孔隙度Porosity-effn39流體置換后的縱波阻抗Impedancen5干粘土含量DryClayFractionn6體積密度BulkDensityn7縱波速度Vpcurven8橫波速度Vscurven40礦化度Salin41氣比重GasSpecificGravityn42氣油比GORn43油比重APIn44油氣標示符Gas/OilFlagn45Brie指數(shù)BrieExponentn47新含水飽和度NewWaterSaturation巖石物理分析的重要任務之一就是從一種流體飽和的巖石速度分析中另的彈性模量特征，速度，等同地，可以說根據(jù)干巖骨架的彈性模量和孔隙流體巖石的速度，這就是孔隙流體置換問題。Biot–Gassmann是大多數(shù)完善的巖石物理模型所采用的流體置換102所示。102輸入和輸出曲線分別為：Page64輸入曲線輸出曲線n1Pressuren36流體置換后的密度Densityn2溫度Temperaturen37流體置換后的縱波速度Vpn3含水飽和度WaterSaturationn38流體置換后的縱波時差Psonicn4有效孔隙度Porosity-effn39流體置換后的縱波阻抗Impedancen5干粘土含量DryClayFractionn6體積密度BulkDensityn7縱波速度Vpcurven8橫波速度Vscurven40礦化度Salin41氣比重GasSpecificGravityn42氣油比GORn43油比重APIn44油氣標示符Gas/OilFlagn45Brie指數(shù)BrieExponentn47新含水飽和度NewWaterSaturation首先是根據(jù)Batzle&Wang等算法計算單相（油、氣和水）的彈性屬性（密度和體積模量，如紅框所示，這個我們在前面地層屬性正演計算中已經(jīng)提及，在此就不再贅述。首先是根據(jù)Batzle&Wang等算法計算單相（油、氣和水）的彈性屬性（密度和體積模量，如紅框所示，這個我們在前面地層屬性正演計算中已經(jīng)提及，在此就不再贅述。（n44為103計算混合流體的體積模量和密度，其中要輸入含水飽和度Sw，這部分即流體置換以前的混合流體屬性的計算（如綠框所示。104（n21）和孔隙度（n22）n47n26和n27之和。Page65105利用n30新含水飽和度計算流體置換后的混合流體的體積模量和密度，如藍框所示。骨架體積模量的計算主要是粘土和石英骨架值的加權得到。骨架值是可以調節(jié)的，骨架的體積組分是由粘土含量得到。105利用n30新含水飽和度計算流體置換后的混合流體的體積模量和密度，如藍框所示。骨架體積模量的計算主要是粘土和石英骨架值的加權得到。骨架值是可以調節(jié)的，骨架的體積組分是由粘土含量得到。Gassmman，考察流體性質的變化所引起的彈性響應的改變。106Page66，根據(jù)實際油藏的溫度、