?作者簡介：熱愛科研的Matlab仿真開發(fā)者，修心和技術(shù)同步精進(jìn)，matlab項(xiàng)目合作可私信。

個(gè)人主頁： Matlab科研工作室

個(gè)人信條：格物致知。

更多Matlab仿真內(nèi)容點(diǎn)擊

智能優(yōu)化算法       神經(jīng)網(wǎng)絡(luò)預(yù)測       雷達(dá)通信       無線傳感器         電力系統(tǒng)

信號處理               圖像處理               路徑規(guī)劃       元胞自動機(jī)         無人機(jī)

? 內(nèi)容介紹

近年來,幾兆赫的陣列電磁傳播電阻率測井(EPRL)已成為探測油氣層侵入剖面電阻率分布的重要手段。然而,要實(shí)現(xiàn)侵入剖面電阻率分布的精細(xì)解釋,就必須借助于陣列傳播電阻率測井資料的反演。多參數(shù)陣列電磁傳播率測井資料的反演問題,具有較強(qiáng)的非線性和多值性,在這方面的研究國內(nèi)外已有很多,但都無法兼顧反演的適定性與反演效率,因而需要新的思路,研究新的反演方法。 反演結(jié)果能較好地反映地層模型參數(shù).實(shí)例資料的應(yīng)用結(jié)果表明,遺傳算法反演出的地層電阻率使測井解釋結(jié)論更加接近試油結(jié)論.

? 部分代碼

function main

%%*************************************************************************

%                     main file of forard modeling                        *

%%*************************************************************************

clear; clc; close all;

tic

L = [0.571  0.9144];                            %發(fā)射線圈至接受線圈中點(diǎn)距離

Delta_L = 0.216;                                         %接受線圈間距

Freq = [400 2000]*10^3;                       %發(fā)射頻率

COt = [0.001 0.002 0.005 0.01 0.02 0.05 0.1 0.5 1];      %目的層電導(dǎo)率

DEt = [1];                                               %目的層介電常數(shù)

SNR = 100; %[1,10,20,40,60,80,100];                 %儀器測量信號信噪比

%%

% a = fopen(’LWD_APS_Homogeneous_analytical_model_inverion_results_mewsuredawgn.data’,’w’);

% a = fopen(’LWD_inversion_methods_comparison_20180428.data’,’w’);

a = fopen(’SNR=100.data’,’w’);

fprintf(a,’%-14s%-14s%-14s%-14s%-14s%-14s%-14s%-14s%-14s%-14s%-14s%-14s%-14s%-14s’,’Freq’,’COt’,’DEt’,’COC’,’DEC’,’ATT’,’ATT_aw’,’ATTINV’,’PHS’,’PHS_aw’,’PHSINV’,’L’,’ID’,’SNR’);

fprintf(a,’ ’);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

ID = 1; %[1,2,3];                 %反演約束類型{1-聯(lián)合約束; 2-幅度比;3-相位差

Rc = 0.12; Nc = 10;                                             %線圈系半徑、匝數(shù)

% sum = zeros(100,1);

% IRE = zeros(100,1);

% ITterNum = 20;

for m=1:1:length(SNR)

    for i=1:1:length(Freq)

        for k=1:1:length(L)

            for j=1:1:length(COt)

                for n=1:1:length(DEt)

                    [ATT,PHS]=EMLOG(Freq(i),Rc,Rc,Nc,Nc,L(k)-Delta_L/2,L(k)+Delta_L/2,COt(j),DEt(n));

                    ATT_aw = awgn(ATT,SNR(m),’measured’);

                    PHS_aw = awgn(PHS,SNR(m),’measured’);

                    for ii = 1:1:length(ID)

                        disp(strcat(’  Freq = ’,  num2str(Freq(i)),’  COt = ’,  num2str(COt(j)),’  L = ’,  num2str(L(k)),’  ID = ’,  num2str(ID(ii)),’  SNR = ’,  num2str(SNR(m))))

                        [COC,DEC,ATTINV,PHSINV]=DE(Freq(i),Rc,Rc,Nc,Nc,L(k)-Delta_L/2,L(k)+Delta_L/2,ATT_aw,PHS_aw,ID(ii));

                        fprintf(a,’%-14.0f%-14.6f%-14.6f%-14.6f%-14.6f%-14.6f%-14.6f%-14.6f%-14.6f%-14.6f%-14.6f%-14.4f%-14.0f%-14.0f’,Freq(i),COt(j),DEt(n),COC,DEC,ATT,ATT_aw,ATTINV,PHS,PHS_aw,PHSINV,L(k),ID(ii),SNR(m));

                        fprintf(a,’ ’);

                    end

                end

            end

        end

    end

end

toc

fclose(a)

? 運(yùn)行結(jié)果

? 參考文獻(xiàn)

[1] 朱信鵬. 基于汽車毫米波雷達(dá)的目標(biāo)跟蹤算法研究[D]. 沈陽理工大學(xué).

[2] 賈善坡. 基于遺傳算法的巖土力學(xué)參數(shù)反演及其在ABAQUS中的實(shí)現(xiàn)[J]. 水文地質(zhì)工程地質(zhì), 2012, 39(1):5.

[3] 信曉麗. 混合遺傳算法及其在電磁傳播測井反演中的應(yīng)用[D]. 燕山大學(xué).

[4] 沈濤, 石曉龍, 陳智華. 基于遺傳算法的核磁共振測井反演算法布點(diǎn)方式[J]. 測井技術(shù), 2017, 41(3):6.

[5] 姚磊華, 李競生. 綜合改進(jìn)的遺傳算法反演三維地下水流模型參數(shù)[J]. 巖石力學(xué)與工程學(xué)報(bào), 2004.

?? 代碼獲取關(guān)注我

??部分理論引用網(wǎng)絡(luò)文獻(xiàn)，若有侵權(quán)聯(lián)系博主刪除

?? 關(guān)注我領(lǐng)取海量matlab電子書和數(shù)學(xué)建模資料