WWW.DISSERTATION.XLIBX.INFO FREE ELECTRONIC LIBRARY - Dissertations, online materials

<< HOME
CONTACTS

Pages:     | 1 |   ...   | 5 | 6 || 8 | 9 |   ...   | 16 |

# «OPTIMIZATION AND MONITORING OF GEOLOGICAL CARBON STORAGE OPERATIONS A DISSERTATION SUBMITTED TO THE PROGRAM IN EARTH, ENERGY AND ENVIRONMENTAL ...»

-- [ Page 7 ] --

3.2 Optimization with brine cycling We now consider cases with varying amounts of brine cycling. These optimizations include additional parameters that deﬁne the timing, duration and pumping fraction of the injection events, as detailed in Section 2.1. The total pore volume (PV) of brine injected over the life of the project is speciﬁed in each case. We consider PVs of 0.01, 0.03, 0.09 and 0.27 as fractions of the core region pore volume. These values correspond to 0.0036, 0.0108, 0.0324, and 0.0972 pore volumes of the storage region (the entire 25 × 25 × 8 block region shown in Figure 2.6), or equivalently, to around 0.14, 0.43, 1.3 and 3.9 times the volume of injected CO2 (at typical aquifer conditions). Three runs are performed for each PV using the same initial guesses for CO2 well placement and control as in Section 3.1. The initial guess for the brine cycling parameters, which is the same in each run, is s0 = [50, 100, 200] years, l0 = [10, 10, 10] years, and p0 = 1/(N K) = [1/12,..., 1/12]. This means that the initial guess entails three cycling events at 50, 100 and 200 years, each lasting for ten years, with equal pumping volume fractions for each well and cycling event.

Results for the HJDS optimizations are presented in Figure 3.8 and Table 3.1.

The ﬁgure displays the progress of the optimization for the best run for each value of PV, while the table presents the ﬁnal optimization results for all three runs. In the table, the value in bold represents the best optimum achieved. It is evident that, for a given value of PV, the optima for the three runs diﬀer slightly.

Table 3.1: Objective function values for the initial guess and optimized solution for each HJDS run.

The best performing run for each PV is shown in bold

–  –  –

Figure 3.8 clearly shows that increasing the PV of brine cycled decreases the optimized mobile CO2 fraction.

It is also evident that the optimizations are quite

3.2. OPTIMIZATION WITH BRINE CYCLING

–  –  –

eﬀective for these systems; for example, for the case of 0.03 PV brine cycling, the initial objective function is 0.252 while the optimized objective function is 0.125. We additionally observe that signiﬁcantly more function evaluations are required for the cases with brine cycling than for the case without brine cycling. This is consistent with the fact that the search space is of a higher dimension for systems with brine cycling (41 versus 24).

We now present more detailed results for the 0.03 PV case. Figures 3.9 and 3.10 show the optimization results for well placement, CO2 injection, and brine cycling.

The well locations are similar to those in Run 1 for the case with no brine cycling (Figure 3.3(b)), though they are not identical. In addition, we again see from Figure 3.10(a) that the well that injects the most CO2 (Well 4 in this case) lies in the east portion of the aquifer.

The optimization results for the brine cycling events (Figure 3.10(b)) are also of interest. There is a relatively small amount of brine cycled around the time injection stops, and then two events toward the upper bound of 500 years. By having such a long time between cycling events, the optimizer allows CO2 to accumulate at the top of

## 46 CHAPTER 3. OPTIMIZATION RESULTS WITH KNOWN GEOLOGY

the aquifer before the ﬁnal brine injection. This would be expected to enhance, if not maximize, the impact of both dissolution and hysteresis in trapping the accumulated CO2. Also of note is that the late-time cycling events in the optimized solution last for only one year, which is the minimum duration allowed.

–  –  –

Figure 3.9: CO2 injector positions and vertically-averaged fraction of mobile CO2 at 1000 years for the best result (Run 3) for the case with 0.

03 PV brine cycling. The inscribed box is the core aquifer, representing the feasible region for CO2 injection wells

–  –  –

0.2 0.2 0.1 0.1

–  –  –

0.4 0.4 0.2 0.2

–  –  –

Figure 3.11: CO2 trapping mechanisms with time for 0.

03 PV brine cycling Figures 3.11(a) and 3.11(b) show the evolution of CO2 trapping for the initial guess and optimized (Run 3) solutions for 0.03 PV brine cycling. The reduction in mobile CO2 fraction is accounted for by an increase of 0.05 in dissolved CO2 fraction and an increase of 0.08 in immobile CO2 fraction. The immediate reduction in mobile CO2 during brine cycling events is clearly visible in both ﬁgures. For example, in Figure 3.11(b), we see sharp reductions in mobile CO2 at around 25 years and again at around 500 years.

The eﬀect of brine cycling is also evident in Figures 3.12(a) and 3.12(b), which show isosurfaces of the mobile CO2 before and after the late-time brine cycling events in the optimized solution. The ‘holes’ in the middle of three of the plumes represent areas where CO2 has either dissolved into the injected brine or been immobilized via residual trapping. It is thus evident that the optimizer ‘steers’ the plumes (by optimally placing wells and adjusting rates) such that brine injection will signiﬁcantly reduce mobile CO2.

All of the optimization results for the best runs at various PV values (presented in Table 3.1) can be combined into a single plot. This plot, shown in Figure 3.13 (the dashed curve simply connects the points), represents the bi-objective optimization Pareto front for the competing objective functions, mobile CO2 fraction and cycling PV. These objectives are competing because, although we would like to minimize both, it is clear that by increasing cycling PV we act to decrease mobile CO2 fraction.

## 48 CHAPTER 3. OPTIMIZATION RESULTS WITH KNOWN GEOLOGY

(a) Before the two late-time brine cycling (b) Ten years after the ﬁnal brine cycling events (486 years) event (511 years) Figure 3.12: Isosurfaces for mobile CO2 saturation of 0.1 in the optimized solution for 0.03 PV brine cycling (Run 3) The Pareto front provides the minimum mobile CO2 fraction that can be achieved for a given PV of brine cycling. It can also be read to provide the minimum PV of brine that must be injected to achieve a given fraction of mobile CO2. We note ﬁnally that, because mobile CO2 fraction can be related to risk and PV of brine cycling to cost, this plot presents the risk-cost tradeoﬀs associated with this CO2 storage operation.

0.25

–  –  –

Figure 3.13: Pareto front for the bi-objective optimization involving mobile CO2 fraction and PV of brine cycling

## 3.3. OPTIMIZATION WITH DIFFERENT OBJECTIVE FUNCTIONS

3.3 Optimization with diﬀerent objective functions In the results presented thus far, the objective function J was the fraction of mobile CO2 after 1000 years. Other optimization time frames and objective functions may also be of interest. We have therefore additionally tested the optimization procedure with several other objective functions: mobile CO2 fraction at 100 years, time-averaged mobile CO2 fraction over the 1000-year project, and the time-average of the total mobility of CO2 in the top layer of the aquifer over 500 years, deﬁned in Equation 2.2 (this objective function will be used extensively in Chapter 4). We denote the fractions of mobile CO2 at 1000 years, 100 years and the time-average as J1000, J100 and Jave. Similarly, the time-average of the total mobility of CO2 in the top layer over 500 years is denoted Jλ. The optimized solutions for the various cases are designated u∗, u∗, u∗ and uλ. Optimizations aimed at minimizing J100, Jave 1000 100 ave and Jλ were performed for a case with no brine cycling (Run 1) and for a case with

0.03 PV brine cycling (Run 3). For the case where we seek to minimize J100 with

0.03 PV brine cycling, some of the constraints and initial guesses for optimization variables were changed (e.g., cycling now must start by 95 years, maximum duration of cycling is 25 years). In each case, the initial guess is denoted u0. The results for these optimizations are presented in Table 3.2.

The table shows results for the objective functions J1000, J100, Jave and Jλ using the four sets of optimized control parameters (u∗, u∗, u∗ and u∗ ). For the values 1000 100 ave λ shown in bold, u∗ and the objective function correspond. Comparing (bold) values for the ﬁrst three diagonal entries for the case with no brine cycling, we see that the minimum mobile CO2 fraction at 100 years (0.330) is larger than the minimum fraction at 1000 years (0.220) and the time-averaged result (0.264). It is also evident that the beneﬁt from optimization is comparable in all three cases, speciﬁcally a reduction in mobile CO2 from the base case (u0 ) of about 0.1. We additionally note that, even when the u∗ determined for one objective function is applied for a diﬀerent time frame, optimization still provides beneﬁt over the base case. For example, if we perform the optimization to minimize J1000 but then assess J100, the mobile CO2 fraction at 100 years, designated J100 (u∗ ), is 0.344, which is only slightly higher than the mobile CO2 fraction in the optimization problem that seeks to minimize J100 (J100 (u∗ ) = 0.330).

## 50 CHAPTER 3. OPTIMIZATION RESULTS WITH KNOWN GEOLOGY

–  –  –

As mentioned previously, the objective Jλ quantiﬁes the time-average total mobility in the top layer of the aquifer. Since mobility is proportional to relative permeability, and relative permeability increases more than linearly with Sg, the optimal solution u∗ will favor conﬁgurations that reduce the overall CO2 in the top layer while λ also avoiding regions of high CO2 saturation in the top layer. Consequently, optimal solutions for Jλ tend to provide improvement to the other objectives, though the converse is not necessarily true. This is because optimal solutions for other objectives may result in regions of high CO2 saturation in the top layer, leading to a value of Jλ that may exceed the initial guess. This is evident for the case without brine cycling.

For the case with brine cycling, both Jλ (u∗ ) and Jλ (u∗ ) are very close to the initial guess value.

Table 3.2: Optimization results using diﬀerent objective functions.

Values in bold are for cases where u∗ and the objective function correspond

–  –  –

3.4 Sensitivity to grid reﬁnement Because the optimizations considered here require hundreds to thousands of function evaluations, it is important that the simulation models run eﬃciently. In our optimizations, we therefore used a relatively coarse model to achieve reasonable run times. It is, however, useful to assess the impact of grid resolution on the optimized solutions obtained by HJDS. To address this issue, we compare solutions obtained during the optimization (using the 39 × 39 × 8 grid-block model for the case without brine cycling) to a reﬁned model containing 64 × 64 × 16 blocks. The ﬁner model was obtained by reﬁning the central region of the original model by a factor of two in each direction.

A cross plot comparing mobile CO2 fraction on the coarse (39 × 39 × 8) model and mobile CO2 fraction on the ﬁne model is presented in Figure 3.14. The runs on the coarse grid are simply solutions at diﬀerent iterations of the optimization for Run 1. The ﬁne-grid runs are simulations of the same well conﬁgurations (placed appropriately in the 64 × 64 × 16 model) and injection schedules. As the optimization progresses, solutions move down the diagonal – from upper right to lower left. It is evident that there is reasonable correspondence in mobile CO2 fraction between the models, though there is an oﬀset of 0.01 or more for most points. We note that reﬁning by a factor of two in each direction represents a relatively modest level of grid reﬁnement. These results nonetheless suggest that the original model is able to provide useful optimization results.

We can proceed a step further by viewing the 39 × 39 × 8 grid-block model as, essentially, a surrogate for the ‘true’ 64 × 64 × 16 model. There are many ways in which surrogate models can be used for optimization. One approach is to ﬁrst perform the optimization using the surrogate model (either for some number of iterations or until convergence) and to then switch to the ‘true’ 64 × 64 × 16 model and perform additional optimization iterations. The idea with this treatment is that the bulk of the optimization can be accomplished using the much more eﬃcient coarse model, though the ﬁnal optimum is determined using the ‘true’ model.

This procedure was in fact tested and the results are presented in Figure 3.15.

Following optimization with the original model (solid curve), we switch to the ‘true’

## 52 CHAPTER 3. OPTIMIZATION RESULTS WITH KNOWN GEOLOGY

–  –  –

Figure 3.14: Cross plot of ﬁne-scale and coarse-scale mobile CO2 fractions for various points during the optimization (case with no brine cycling, Run 1) model (dashed curve).

There is an immediate jump in the objective function, followed by a gradual decrease with iteration. The optima achieved by the coarse and ﬁne models are close, 0.22 and 0.23 respectively. It is also worth noting that only two of the rate variables changed during the ﬁne-scale optimization (the well positions did not change), which again indicates the eﬃcacy of optimization on the coarse model.

This procedure could be continued with increasingly reﬁned models if necessary. Additionally, ﬁner-scale permeability and porosity descriptions, if available, could be introduced as the grid is reﬁned.

3.5. SUMMARY

–  –  –

0.3 0.28 0.26 0.24

–  –  –

Pages:     | 1 |   ...   | 5 | 6 || 8 | 9 |   ...   | 16 |

Similar works:

«Why Governments Won’t Invest Sir Vince Cable Special Paper No. 33 March 2016 Centre for Economic Performance Special Paper Sir Vince Cable is Professor in Practice on Global Finance, LSE Institute of Global Affairs. All views are the author’s own. The Centre for Economic Performance and the ESRC has no political affiliation or institutional view and the author writes in a personal capacity. Abstract The UK Government is currently exerting a sharp and tightening squeeze on public investment...»

«The Dismal Performance of Non-financial State-owned Corporations in Bangladesh Tanweer Akram* January 28, 2003 *Academic Affiliation: Department of Economics, Columbia University, 420 West 118th Street, New York, NY 10027, USA. Correspondence Address: 5550 Columbia Pike, Apt 803, Arlington, VA 22204, USA. Email: ta63@columbia.edu Acknowledgments: This essay is based on documents and data made available by the Monitoring Cell of the Ministry of Finance, Government of the People’s Republic of...»

«JILL J. McCLUSKEY School of Economic Sciences Pullman, WA 99164-6210 Distinguished Professor of Sustainability Phone: (509) 335-2835 Washington State University Fax: (509) 335-1173 mccluskey@ wsu.edu https://scholar.google.com/citations?user=Qbp4C-oAAAAJ&hl=en DEGREES Ph.D., 1998, Agricultural & Resource Economics, University of California, Berkeley Fields: Economic Theory, Environmental & Resource Economics, Industrial Organization Dissertation: Environmental Contamination and Compensation,...»

«Volume 1 Issue IV Dec 2013 ISSN 2321 7065 Women and Human Rights: Some Reflections in North East India Dr. Th. Binarani Devi Department of Political Science G.P Wome n’s College, Imphal 795001 Emailthokchom2009@gmail.com Abstract The whole fabric of society depends upon human rights. Human rights are essential components of sustainable development and therefore development of a society is not possible without respect of human rights. Respect of human rights makes freedom. They are part and...»

«60 South Sixth Street, Suite 3900 Minneapolis, MN 55402 612-661-3770 www.marquetteam.com March 27, 2015 This Brochure provides information about the qualifications and business practices of Marquette Asset Management, Inc. [“Marquette”]. If you have any questions about the contents of this Brochure, please contact us at 612-661-3783 or jean.gullikson@marquetteam.com. The information in this Brochure has not been approved or verified by the United States Securities and Exchange Commission or...»

«TERMS OF BUSINESS AGREEMENT Law and Language This Terms of Business is subject to English Law and the jurisdiction of the English Courts. We will use the English language for all communications, the contractual terms and conditions, and any information we are required to supply to you, before and during the duration of the contract. Acceptance The purpose of this agreement is to set out our professional relationship and detail the services we will provide to you. For your own benefit and...»

«1 Finance Committee Minutes Date: March 6, 2010 Time: 8:00 a.m. Present: O’Connor, Kupcinskas, Hamilton, Gauthier, Bartosiewicz Absent: Hebert, Girard Meeting called to order at 8:00 a.m. Bob Spain, Town Manager, was in at 8:00 a.m. Brian Turbitt, Finance Director, was in at 8:00 a.m. Richard Handfield, Police Chief, was in to discuss his budget at 8:00 a.m. Michael O’Connor: We are here today to discuss budgets with several department heads. The Town Manager was kind enough to prepare a...»

«ISSN 22775846 The International Journal Of Management THE INTERNATIONAL JOURNAL OF MANAGEMENT A Strategic Resource Analysis of Sustainable Competitive Advantage in an Australian Service Business Environment Dr. Fouad Ichrakie Assistant Professor of Management and Finance School of Business Administration, Rafic Hariri University, Beirut, Lebanon Abstract: In the last seven decades, many scholars tried to explain the variances in performance amongst various firms and industries. One main...»

«ICABE 2014 • International Conference on Applied Business & Economics EDINGS October 23rd – October 25th 2014 European Cultural Centre of Delphi (ECCD) Greece          1    ICABE 2014 • International Conference on Applied Business & Economics CO NT E NT S PROCEEDINGS LIST OF PAPERS CONFERENCE PROGRAM ABSTRACTS APPLICATION FORM OF INTEREST   2    ICABE 2014 • International Conference on Applied Business & Economics EDINGS PROCEEDINGS October 23rd – October 25th 2014 European...»

«NEW MEXICO TRANSPORTATION BY THE NUMBERS: Meeting the State’s Need for Safe and Efficient Mobility January 2014 202-466-6706 www.tripnet.org Founded in 1971, TRIP ® of Washington, DC, is a nonprofit organization that researches, evaluates and distributes economic and technical data on surface transportation issues. TRIP is sponsored by insurance companies, equipment manufacturers, distributors and suppliers; businesses involved in highway and transit engineering and construction; labor...»

«Pharos Journal of Theology ISSN 1018-9556 Vol. 96 (2015) Copyright: ©2015 Open AccessOnline @ http//: www.pharosjot.com Religiousness: An Inter-Disciplinary, Multi-Demographical and Multi-Regional Study Prof. A. Grobler* Area Head: Leadership and Organisational Behaviour Graduate School of Business Leadership University of South Africa Email: grobla@unisa.ac.za Telephone: (office) +27 11 652 0277 Dr. S. Grobler Senior lecturer: Department of Industrial and Organisational Psychology University...»

«Corporate/Commercial & Corporate Finance Practice October 2016 Contents Contents About Us Our Corporate and Commercial Expertise Christopher Williams Katherine Chuang Brian Ho Eddie Lee Chia Ching Tan Jill Wong Christopher Yu Denise Che Antony Yung Gilles Cardonnel Bertie Mehigan Our Team Recent Representative Matters of our corporate team members 1 Howse Williams Bowers is an independent Hong Kong law firm which combines the in-depth experience of its lawyers with a creative, forward-thinking...»

<<  HOME   |    CONTACTS
2016 www.dissertation.xlibx.info - Dissertations, online materials

Materials of this site are available for review, all rights belong to their respective owners.
If you do not agree with the fact that your material is placed on this site, please, email us, we will within 1-2 business days delete him.