如何在 MATLAB 中从大型文本文件的重复格式读取数据?

发布于 2024-09-16 07:44:36 字数 2563 浏览 10 评论 0原文

我想读取水的平均饱和度(%)数据,如下所示。该数据是大文件的部分形式,但是平均水饱和度 (%) 仅以给定格式重复。

Average Pressure 
   Total Pore Volume  psia      3884.9                                                                                              
   HC. Pore Volume    psia      3884.9                                                                                              
 Average P/Z 
   Total Pore Volume  psia      4457.8                                                                                              
   HC. Pore Volume    psia      4457.8                                                                                              
 Average Saturation %
   Oil                          84.911                                                                                              
   Gas                          .08873                                                                                              
   Water                        15.000                                                                                              
 Percentage Recovery 
   Stock Tank Oil               .02211                                                                                              
   STO as a % of Mobile Oil     .02891                                                                                              
   Total Gas                    .02034                                                                                              
   Water                        62e-12

我试图通过使用 readline 来做到这一点。 m 函数,但不幸的是平均含水饱和度(%)数据的位置不是由行号固定的。不同型号的类似输出文件的行号会发生变化。

这就是我想要做的:

%# Reading Water Saturation (Sw) data from output (.OUT) file of reservoir model
    Sw_LineNo=[554,968,1120,1272,1424,1576,1728,1880,2032,2184,2336,2488,2640,2792,2944,3096,3248,3400,3552,3704,3856]; % This column vector contains the line numbers of the .out file with Sw values for year 1 till 20

    for i=1:size(Sw_LineNo,2)
    read_value=readline('ReservoirModel_ExplorWell_CMGBuilder.out',Sw_LineNo(i)); % read_value stores values in form of string
    Swav_Data_E_W(i,j)=str2num(read_value(33:38)); % converts the required portion of string (Sw value) to number
    end

现在,如果我的模型 (ReservoirModel_ExplorWell_CMGBuilder.out) 发生变化,那么文本文件中水的平均饱和度 (%) 所在的行号也会发生变化。因此 Sw_LineNo 会根据不同的模型而变化,并且我有大量模型。

请建议读取水数据的所有平均饱和度 (%) 的正确方法。

I want to read data of average saturation (%) for water as shown below. This data is a partial form of a large file, however the average water saturation (%) REPEATS itself in the given format only.

Average Pressure 
   Total Pore Volume  psia      3884.9                                                                                              
   HC. Pore Volume    psia      3884.9                                                                                              
 Average P/Z 
   Total Pore Volume  psia      4457.8                                                                                              
   HC. Pore Volume    psia      4457.8                                                                                              
 Average Saturation %
   Oil                          84.911                                                                                              
   Gas                          .08873                                                                                              
   Water                        15.000                                                                                              
 Percentage Recovery 
   Stock Tank Oil               .02211                                                                                              
   STO as a % of Mobile Oil     .02891                                                                                              
   Total Gas                    .02034                                                                                              
   Water                        62e-12

I was trying to do it by using readline.m function, but unfortunately the position of average water saturation (%) data is not fixed by line number. The line number changes for similar kind of output file for different models.

This is what I was trying to do:

%# Reading Water Saturation (Sw) data from output (.OUT) file of reservoir model
    Sw_LineNo=[554,968,1120,1272,1424,1576,1728,1880,2032,2184,2336,2488,2640,2792,2944,3096,3248,3400,3552,3704,3856]; % This column vector contains the line numbers of the .out file with Sw values for year 1 till 20

    for i=1:size(Sw_LineNo,2)
    read_value=readline('ReservoirModel_ExplorWell_CMGBuilder.out',Sw_LineNo(i)); % read_value stores values in form of string
    Swav_Data_E_W(i,j)=str2num(read_value(33:38)); % converts the required portion of string (Sw value) to number
    end

Now if my model (ReservoirModel_ExplorWell_CMGBuilder.out) changes, then the line numbers where the average saturation (%) for water lies in the text file also changes. Thus Sw_LineNo changes for different models, and I have large number of models.

Please suggest correct way to read all the average saturation (%) for water data.

如果你对这篇内容有疑问,欢迎到本站社区发帖提问 参与讨论,获取更多帮助,或者扫码二维码加入 Web 技术交流群。

扫码二维码加入Web技术交流群

发布评论

需要 登录 才能够评论, 你可以免费 注册 一个本站的账号。

评论(1

夏日浅笑〃 2024-09-23 07:44:36
%# Reading Average Water Saturation (Savw) data from output (.OUT) file of reservoir model
    fid = fopen('ReservoirModel_CMGBuilder.out'); % open the file

    dotOUT_fileContents = textscan(fid,'%s','Delimiter','\n'); % read it into one big array, row by row
    dotOUT_fileContents = dotOUT_fileContents{1};
    fclose(fid); %# don't forget to close the file again

    %# find rows containing 'Average Saturation %'
    Swav_Starts = strmatch('Average Saturation %',dotOUT_fileContents); % Swav_Starts contains the line numbers wherever 'Average Saturation %' is found
    nSwav = length(Swav_Starts); % total no. of Swav values will be equal to the total no. of 'Average Saturation %' read from the .out file

    %# loop through the file and read the numeric data
    for w = 1:nSwav 
        %# read lines containing numbers
        tmp_str = dotOUT_fileContents(Swav_Starts(w)+3); % stores the content of the 3rd row from the row containing 'Average Saturation %' in form of string
        tmp_str = tmp_str{:}; % store the content of the string which contains Swav, as well, in form of a character
        %# assign output
        Swav_yearly(w,j) = str2num(tmp_str(30:35)); % convert the part of the character containing Swav into number
    end

现在 tmp_str = dotOUT_fileContents(Swav_Starts(w)+3); 生成以下字符串:

Water                        15.000 

如果我尝试使用 str2num 将其转换为数字,则会得到一个空矩阵。因此,我选择包含饱和度值的字符串字符(此处为15.000),然后将此字符更改为数字,如下所示,给出平均水饱和度的值:

str2num(tmp_str(30:35))

如果有人有任何方法,请告知从字符串中提取数字,而不像我那样选择字符。

%# Reading Average Water Saturation (Savw) data from output (.OUT) file of reservoir model
    fid = fopen('ReservoirModel_CMGBuilder.out'); % open the file

    dotOUT_fileContents = textscan(fid,'%s','Delimiter','\n'); % read it into one big array, row by row
    dotOUT_fileContents = dotOUT_fileContents{1};
    fclose(fid); %# don't forget to close the file again

    %# find rows containing 'Average Saturation %'
    Swav_Starts = strmatch('Average Saturation %',dotOUT_fileContents); % Swav_Starts contains the line numbers wherever 'Average Saturation %' is found
    nSwav = length(Swav_Starts); % total no. of Swav values will be equal to the total no. of 'Average Saturation %' read from the .out file

    %# loop through the file and read the numeric data
    for w = 1:nSwav 
        %# read lines containing numbers
        tmp_str = dotOUT_fileContents(Swav_Starts(w)+3); % stores the content of the 3rd row from the row containing 'Average Saturation %' in form of string
        tmp_str = tmp_str{:}; % store the content of the string which contains Swav, as well, in form of a character
        %# assign output
        Swav_yearly(w,j) = str2num(tmp_str(30:35)); % convert the part of the character containing Swav into number
    end

Now the tmp_str = dotOUT_fileContents(Swav_Starts(w)+3); produces the following string:

Water                        15.000 

If I try to convert it to number using str2num, then I get an empty matrix. So I pick the characters of this string containing the saturation value (here 15.000) and then changing this character to number as follow, giving the value of the average water saturation:

str2num(tmp_str(30:35))

Please advice if someone has any method to extract number out of the string without picking the characters as I have done.

~没有更多了~
我们使用 Cookies 和其他技术来定制您的体验包括您的登录状态等。通过阅读我们的 隐私政策 了解更多相关信息。 单击 接受 或继续使用网站,即表示您同意使用 Cookies 和您的相关数据。
原文