% network3D(2060,1580,10,1,1,3  );
%REFINE 此处显示有关此函数的摘要
%   此处显示详细说明
% % % well=cell(4,1);
% % % load('perforation1.mat');welloc1 = [perforation1, 7.5*ones(size(perforation1,1),1)];
% % % load('perforation2.mat');welloc2 = [perforation2, 7.5*ones(size(perforation2,1),1)];
% % % load('perforation3.mat');welloc3 = [perforation3, 7.5*ones(size(perforation3,1),1)];
% % % load('perforation4.mat');welloc4 = [perforation4, 7.5*ones(size(perforation4,1),1)];
% % % well{1,1} = welloc1;
% % % well{2,1} = welloc2;
% % % well{3,1} = welloc3;
% % % well{4,1} = welloc4;
% % %      
% % % for i=1:4
% % % line(well{i,1}(:,1)',well{i,1}(:,2)',well{i,1}(:,3)','color','k','LineWidth',2);
% % % hold on
% % % end

% load('injectionWellsCoord.mat')
% injectionWellsCoord = [45,45;
%     655,45;
%     45,625;
%     655,625;];
% for i=1:size(injectionWellsCoord,1)
% line(injectionWellsCoord(i,1)*ones(2,1),injectionWellsCoord(i,2)*ones(2,1),[0;20],'color','k','LineWidth',2);
% hold on
% text(injectionWellsCoord(i,1),injectionWellsCoord(i,2),20,'注水井','FontSize',12,'FontWeight','bold');
% hold on
% end

% wellbore
% wellbore = [250+100,100+100,46;
%     550+100,100+1700,46];
% % wellbore = [50,352,10;
% %     950,352,10;];
% plot3(wellbore(:,1),wellbore(:,2),wellbore(:,3),'color','k','LineWidth',2.5);

% wellbore = [450,520,10;
%     450,1420,10;];
% % wellbore = [50,352,10;
% %     950,352,10;];
% plot3(wellbore(:,1),wellbore(:,2),wellbore(:,3),'color','k','LineWidth',2.5);
hold on
% wellbore = [1175,370,10;
%     1175,1170,10;];
% % wellbore = [50,352,10;
% %     950,352,10;];
% plot3(wellbore(:,1),wellbore(:,2),wellbore(:,3),'color','k','LineWidth',2.5);
% text(wellbore(1,1),wellbore(1,2),10,'压裂水平井','FontSize',12,'FontWeight','bold');
hold on
% hwell=cell(1,1);
% hwell{1,1}=[385,1300,15;
%              385,400,15;];
% hwell{2,1}=[955,1330,15;
%              955,330,15;];
% hwell{3,1}=[1635,1275,15;
%              1635,580,15;];
% for i=1:3
% line(hwell{i,1}(:,1)',hwell{i,1}(:,2)',hwell{i,1}(:,3)','color','k');
% hold on
% end