{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "\n", "*This notebook contains course material from [CBE20255](https://jckantor.github.io/CBE20255)\n", "by Jeffrey Kantor (jeff at nd.edu); the content is available [on Github](https://github.com/jckantor/CBE20255.git).\n", "The text is released under the [CC-BY-NC-ND-4.0 license](https://creativecommons.org/licenses/by-nc-nd/4.0/legalcode),\n", "and code is released under the [MIT license](https://opensource.org/licenses/MIT).*" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "\n", "< [West Virginia Chemical Spill](http://nbviewer.jupyter.org/github/jckantor/CBE20255/blob/master/notebooks/B.01-West-Virginia-Chemical-Spill.ipynb) | [Contents](toc.ipynb) |

\"Open" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Ajka Alumina Plant Spill" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Summary\n", "\n", "This notebook demonstrates the analysis of the 'Red Sludge' tragedy that took place in Ajka, Hungary, in 2010." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Background\n", "\n", "On Monday, October 4th, 2010, an earthen dam failed at an alumina ore processing facility near [Ajka, Hungary](http://www.theatlantic.com/infocus/2011/09/a-flood-of-red-sludge-one-year-later/100158/), resulting in the release of red mud over a large area in western Hungary. According to the [New York Times](http://www.nytimes.com/2010/10/06/world/europe/06hungary.html), a flood of 700,000 cubic meters of sludge swept cars off bridges, damaged roads, and flooded businesses and homes throughout the affected region. There were reports of up to eight deaths and hundreds of displaced families. The red mud entered creeks about 45km upstream of the Danube raising the potential for devastating environmental damage. At the time there were fears that other earthen dams at the same processing facility were showing signs of impending failure." ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/jpeg": "/9j/4AAQSkZJRgABAQAAAQABAAD/2wCEAAYEBQkJCQkJCQkJCQkJCQcHCQkJCQgHBwcHDAgMDAwICgoOERkTDg8XEAoKFCAUFxobHR4dDRMtJCEcJBkcHRwBBwcHCwkLDQoKDRwODA4cHBwcHBwcHBwcHBwcHBwcHBwcHBwcHBwcHBwcHBwcHBwcHBwcHBwoHBwcHBwcHCgcHP/AABEIAWgB4AMBIgACEQEDEQH/xAAcAAABBAMBAAAAAAAAAAAAAAAEAgMFBgABBwj/xABSEAABAwIEAgUEDAsGBQQDAQACAAMSBCIBBTJCE1IGB2JykhEjgqIUM0NTk7LBwtHS0/AVFyExVIGRobHh8hZBUWHi8QhjcYOzJCVzozREwxj/xAAbAQADAQEBAQEAAAAAAAAAAAAAAQIDBAUGB//EACwRAQEAAgIBBQACAgIABwAAAAABAhEDEiEEExQxUQVSIkEGkSMyM0JhcdH/2gAMAwEAAhEDEQA/APPmPQyr98p/E59RZj0Nq/fKb4Rz6i6CWCRiuSc+Tf24oP8AY2r98pvG59RRz2SOtnwyNqWHaKP8F0zGSrGZ+2ld2k/eyOcc2iKDojUvlFtxj0ic+QFM0fVfmLxRF6j/AFuVHyNqa6LEIueiukdHzGQ9oS8S0487knl45jNxy0eprNsRlx8v+FqvsVrDqazX9Jy/4Wq+xXe2HSiIxS1tfDluVcC/Ezmv6Tl/wtV9is/Ezmv6Tl3w1V9iu/eVKCMtQqdl3rzdn/VlmGXtC+8/QkJFEeG7UEUv1t4Krlk7uBeSTXlw/wA3PoXo/rOZJyhGOkXZS7K4m6xIii4JeiXyI20wu/tWfwO7zN/tL6Fn4JPna8RfQp9yncErbkyTZd0ku1a9Yhzykx/ObXiL6E0GXmW4P2l9CmTAkjAS5SHtRtU96cwgBvJ3cfySD9pfQpeg6E1j+h2mw7xu4fwBFU1K5aQkXhV16L0xCJSKJeqp9zJdwx0rdD1S5nUEQg/Q4RGWMnagR/c0isepbNv0nLvhaz7Bda6NVrbAmRERWoj+0BHa0Mi3SERJaTKuS27cbx6ms1/Sst+FrPsFmHUzmv6Tl3w1V9iuxY1lTtGKQdXWR2j6Ke6c7VyHDqXzX9Jy74Ws+wWfiYzX9Jy34Wr+xXWCr6uKb/CNSJXEIqeyphnf9OVfiYzX9Jy34Ws+xTdV1PZq01i6VRl+Ij/cLtVL9mLK6zhmT/NKMkNUZk+Q8IpRldbtT7NMeLP/AHHDnOhVWJeTiU3je+om3uh9WGo2PE59RdXqWGhc7XrJt9tiJD82SO1Zz8rl39jKyQDNjzgyEpORj4FmHQ2pnw+PSy7731F1HBhshbEdvekmsMubJ1xy6JWiOkpD2f7kdqrw5tW9CqpgMXCfpCEY6DeIv3treHQis4XHwcpib8glrdl8RdNfo2hCLsiEubUg2aWJcNt3zRN6ZCjdGooR9Aq6LZ8WlLi4SHDBx2XxE010HriKPlY/WTkfiLq9FgLTQtlHiD2pRRJvgTgxIRjq2o3Rpymu6uswYHAnHKXyFym98rabPoBXCTY8SlxxdGYxcext8C6t0yIn2mhEo3alA4kYcIpSJq3mtS7Vcwlc/qOhlU2RCTlMRD+cRN3Eu9o/Mh3ei1SO5rTL8hOfUXVcGxKpbMtLolLs/Qt1NC04+2I3NiMXOb/dReSjpHHmMldcKIk3y3EX0KXb6DVhDIXKYu6b31FaKfL+EZRHzZEXejL6FZKKmINN0hEolyqPeyP23MMOhVVd5XqUI8xvXd3yNoWn6MPuFEXGPE59RdTzBwfOW7YqCySkGV1pCSPep+0rOHQKt99pPG99RZ/YKulHiU2HpvR+IupRISlq7Kbq61uBStLu/Qj36XtVzY+ruvEcT4tJEe299mhy6DVnPTF3Se+ouj/hhuJCJXFzDJOZJjIiIrhK3tCX38i1xzuRZYac3Y6BVzmk6bVGOJuy+Ip+j6mM3dESF7Lxl747VB//ACXWujmSEPnHBtIpDLUrvRlERGOm2K2xlsc2XJq+HnjDqKzn9Jyz4es+wTw9Qmd4/wD7WV/D1n2C9INYyR1PttkrmCPcrzPh/wAPeefpWVfD1v2C3/8A56z39Kyn4et+wXqIcIreOKrpB7leW8f+HvPf0vKfh637BIx/4fs8/Ssq+HrPsF6jLFMuEjpBeWvL34hc6/Ssr+HrPsE3+IrOf0rK/h6z7Bem8be0hDxISuHUn7cKcteW6/qlzNjXUUGPddqvlaQbHVpmDn5n6P8AW5UfZruXTB7b2lAUlQQ2yjtR0gvLXMvxW5j7/Q/C1H2aCLoBWYC4XFpfNCRl5x6RCIyt82u1C9ISGWnSq/mFrT8iG5p//wARJdImc12jcUg0vFMOOd5eVI9UotyrWbYNzIt3hU8VSKgc5ISKQ6opniN6MFJ3vCW5dL6On5xrukK5h0bcHjjHliuiZU4LZSL/AFK8Lo88dzS9BbHtJaro9IBESHtW27UoMzfc9qEi7o2rbuyx9HnknsSjqtHmUfmGcNNbhktU9BmVT7mQjq0iKcDoBVukJOXelaqcvNw3C6V7OM2CsYcaLTIS+/kVZp+jjBFLhiXduXccg6EtMD59tsiKMhIRJTzHR6jGMWGh7o6lUicfDgFL0TG4RYIpaZDcpij6CyGRU2rmFdwxyxsdIiPhTg0I9lPqrtXGm+r6Q+0N+kK0fV2W1hrwrs40w81q2VEN1xau8l0HuVxnDq5dj7QOq3ai2+gtS1dwtQxIZCuwk2SaxYIh1JdNq7uN1PR6spGjccaHhiJSLdFRXRkeOX/phkQlpJdxr6An2ibK4XBISlHSoLoz0OaoScJsRKVw3ff8qcx0zsVhqmqxKTlMMe6iDccjL2Np0+bGKv8A7GKXtfiFbwoeZtvuimcyyjnYVjelymIebzacwqsvL3CQ7iIdKv2NIPvTXpCkfg5jcw16IxS0v3soo7buVF7m2MrbkR7AyriCPm4ldq1K2OZNSOaqZrl0x/gqJ1i9H3W76ESGIja2RW/t/N+RKRXyc01/YjL3/OCIxLtbUhzq4o+Xu3fyVKyDpBmVIIi4LpNj2Zetj/cus9Gs6brGh2uCN0uZPqz7XJVvxcMD4ltzq6YiUR9b+S6BisRobcuzPq3EmpDcQ6ZFd/BVXor0JN06kS9wdjLsyJd6xwSGKdsSMRbbGUSKIxkXa8n50aHZyDNert+ROAI3bZaVFOdXlYO2XpLuxuDplpimnKxsdRfFTmI7uHD0IrOGQut3bbhuFRNT0NrmxIuFaXqiu+PZ1SbjElHVWfZfuIfSEUdFTl088VlK+3ESEhIdyw2HWh4kZStkuy5jV5NUyEuAJEMRjbcuZZ5kbrRlFwnGykTZSt7v7FzcuNxbcNnJVYedIo7Ue3VaSlpGJIqn6PvmIuRG7uqUaydpvWQiXKueV6OPDJFZcORXSKSbaovOCQjFXNjK2CutJSDdA0PufxVNV7UnlUSYdLSPpJNLk5uFdd3tKupA1yx7sUxIR1DbzaU8cd0r1insZZSC642IiRNXOFK4RVq6OZKBENTw4sDaI++9olXM0aEX+KwQ3e2DzfSpIM/cabFtsrR22xXo8eGnl+o5fOnTWGGxEbe6nhHcufUXTE4jxZDcO21WOk6UUzjdzgiXaW81HFasrJI9jHxKsUvSKm0zGSLPpLRtDJx0YkqJYpS3LXlUTlvSGhfIRbdC7tCpdp9gpecHxI2mtEdqYxJFatOn0U3i1cmAuGpC5qQiMh1CpR1rsqs9KHyAezdJGxpzjpTWlIi7SgmamQlK0tQklZ9VE46QjK0u8iXqZsm2xk22XDuIktnMdhgribK60SG1Qeb5sRC4I7m3xLwkjXKQGy86+LktI8TSSkKvK2m6ZwhbEi4D5f8A1FpQc4/Ln+OdOaRld3Vpx6pc0y2q5tdG2uYUe3lVM3aQivM6vV251hTVJc3iSTy58iiRejJdRCmph/pSipRK8Rt9HSjqO0UXIacmiHzZF2hirzR4EQ6CkjsuoDdcEW2xKRdldFyro7TUzfHqbYjIh1EPorTDj2j5HVWej/Q918hMhERjKUl0fK+i9NTCPEbGRbiGXy/kwVWzDpsLUmqQSiOkhiKrlX0kzCpIZPkPpRW84i5fXZWaxdfdqqamEtMdNooQuk1IO7dyrlYE+ftrpF6S05GI2lq5lfRzTO5ecnY6HN2Kn2u7mRr7wAMit5uyqX0AoSHiGWmIxuuSus7N+BTC0255wyjqjsR/9Jyukz/aaju856pJ1vpHSFv9VcRo25CJEWoRlcjm8BHcUu8Svoy7Oxt51TFuLwp38JsbSXIGH4yjKXNKSZOtduiXrI6n2dobzRgvdIlKNwp1uva2kJdqMblwTiujIhMhKUtSf/CdYAl50ox5t3ZR1Pu78ZtxEpS7qG4jY7lWervMXanL5PyJwRIZFq0rmnSLpTWMVhNtukQjq2jq0/lS0qV3XiDzLYOdr0l59/trmG0yLwj/AARTfTasEZERWlzCp0XZ3vEh7y1auJ03WXUjrl6MVNU/WS2UZC5Hup9T7R1HEU1UsC4MSES5lSqHrDpiIW4kMuZWiizhh/SQ3afvipsG23clpCtJofDqS6DKaamKTTYgRd5GgUvSGSXipqoTjgh6ivBu0tq3mDhNtOENpCMhXBs9z+rN10eKWoh9ZOQsq7JmPSNhoZSt7tyrGYdPWhHzbly5gIPuCJOO2lpFHUNM0VpCKvTPsmq/pu+egiHwqFe6QVLvurnrInCmb5Uv2IPKpXijvZlS57o4SYxF0vbDLxKdbY7q09SSEtPiuRsqrjjJFuKXaJTLFWVS17GIouiPm3Cu2aP8kxjaSacEpCY2kJakZyZTyOHO4UIVc+wPDKUhtjqiShqysfIpcTVarnXMtuttuiNxDEuYiFCt0LRe2CNq8/PxX0XDlOTGaVZiufHS4UUT+EKnURF4vWVxboqS0uH6qj87fo2mi4YtkUS0rOY20uX/ABiCZzYhK4iLsyTmY52RN2iQiQ8w2qtuV7ZRiMeaSSJGW0o8q6MePTgvL2E4PERXERJ5yot3IcadwtIopjKn3NI3d6K6OznuFyNcUo3ERJOLpcxKR/s/VltHtXDpSv7NVctviFL3B8XYWnzAhHV6ydxrActdEi9K1GNdFqkrpDy6rlhdFKkd4il7pfEplnM22Ci3IbbbkS10pfHS654pJhzow/uISjquWH0cJOcyb6WrVl3T50RES4heH5VOUPTkSG4SkqK1lkSHujpUs3SNjqa9Za+7Ki+mq+f2pbIZDLTaKrufZ1xRK0u9aop0IjIVFvuFG6VyV5FY+nU3Ps0cYdJ1sdUtSrONVXVzhRfJsR2yVjzFibpSEreZKpsW2JDbcovJpr8dE5Nl74FJ90nB2yLcp7Ns2k04IkUeE6OrsEKjK6vKVpKKfdIm3LtpfEU+8r41dWBhsdQ3ClOONiWkVGv5kPMKCdzdsZXdnSsdvVw9JN+U2dQ2W0RR+U0DtS4Lbeki7qB6LZY7XEJDa3K4o7VbavMmMracCkIeKNpORld2fKtccduX1V4+KJUHqbKmyExEn46RuGX61Ta7M6uscJwiIWiK0ZRG3s4f3IGoMny4jjhOEVxERJ1uoiMRtHSunDB4Vtyy2fabER7RLaDJwkrE1ouXQ0D5fWRVAHHcFsSuJQ+JR1Kz9XNALtXxBu4YlaVyKO7quTU4sMNjpi2MlxfrCzHj1hCNwtkS7Hn9WNNSOu6bbV5zxqHHXTcJy5wiItKWMRl5SVDiQoocRLVtQpYRiIlcm8XC5lonaRdqBEdyjzqZFuUPm2YEMRkXoqIwzMRIiIiVaDoOLdolptUFnFURONtiUYkMi5hUQ30idK2XZHuoTF0n3BkUtXZU3wcei+rLCVJISERIikMeUFxvpU8JVz//AMnqyXZ+r2nFvKhIRja6Uuaxea+ldc4Nc7Fz3Qo+NQ0WfLxEpaUZwuYRiqnkOZHEic5lM09eJXSkhFSRNN8ooqnomiErVGYVbfMjcvqWyG0pERJlKIqssAtNtvKo9yqfpLhMrS2kQkrSyI280bu6oTN6cSac7JKMoqV0bqx6UDVti0fEIxlqL7/3LoGBLzV0EzIqbMGIyuMmiHsxXo6mOTbZcyjJrjTtSMm3B5hJebOkuMayrHTF0l6XPCVq839YDPDzWsbH3ySnGnkVlrkmm+yj6fVJQmVFaIqTwciUluw1pKz8QrROR7I+qh23Lbk/ASEVjm0wrY4iW5IHAeb0f91sxEd3hTThjHclFhnSGUSQ9SMpD85Kfxut5u8hscZlK4Y22qrNoF5fVCJCJSjpGOmSZr5CQkNpCSG8pNlaUiEtRKQdeaEWnSlArTj76uXPDVel6Pm/0jHa90dJXF2lDY5TUuyLiDdtkpQ81Ybd2xt2qwt5uwLYxEfDqU4ZYx183bJTsp6MOuOjKMZXXSV3YySmbGMRlttWqbMBKKJ9liP9Sq5M8OMptlseXsxEVvyDt1eFJwqx3aU42+BEs7k1x4zJCZdr0oph8XNpR9KSksHG+X1UjEWylEUSruKMYacEpESfwFzmtRgCMZcqT5BK5CdaDYtDzFLvLXAiiPIKZddEf6kSA2bS06+KEqa0R5lD1WZjzI3oSS+EhmFeIqCzHMhEVH1+ZiNxSVYzHNgKUdSO9K8UG5nmpEo114SuUa8UubxIjL2ZKbV446EBItKx6kciXLwy+KrVkeVDG7VHl+MpTPqJsKZwmyEiFp3TErYKYMp4QQtDKJSl2RkrB0T6NFXODARIZRKS3n9ONG7wIiRQkXMPd/xVQpOmr9HViDUhEiuKUdX/AEw/wW+HHuufm9brxHY+mPSBjIaH2IwI8UhuIY6iP/Jc/oa8nGhfMpEdxXfSqV0+z4q2pGJEUoyIpf7KfyV7zDTcrR1Lsxx08jktzu7VhbzMtItyHtInjEVxDH0lVqyp4UiGRCox7pE4JW6dwyWkumVjofHGN1skBVZm2FokPrfIq3UZyJNDd94qAfd4soulIrYpXJUxXvHOB3aV2XqybEaQX9PEG2I9teb6Nq5seIRDIRjFerOjNHwMvo2xtGIlGKntsdVZ66MzcYy8RkQk4QiMS+MuLZHWykJRkIj2lef+ICvFxximlpIi5tID+xc0oMBaHVKUfRTnhF8Lh7NHcWlNP1zcZCMlUnXS5reW1Cv17g2q+yP9pyqwFwpSj3lF1QCW4VHY1zhavjJJPS3J7VpJYRjaSlMuMR8JXKrtOx3bpI1qqtLw9pGRyPUfQkv/AGYIuS82UuzavNPSF3/1jso+2FqLVcvR/QQY5MIjbJsi0x2Ly90ydjWO8pEQ2lK5ZtKl8rfFthwoiMiUc5mhC4UZfN/cm6E/MEO7aop0ymRaRJOJsS9Pn0d1yl8sz6MpW+lcqc2A6tyJAoy9VFRp03Lek4EO6Q9oU+7njZCY8wl2tq5SdQ4MYo2lekMiIpJHJVhpa2NS0Q2xdEpal6m6OvToaRzmbG70V40dfi5q3CS9bdAKjiZRQlzNj9/yYLPNrItmGMV5/wCtxggzVwo2uiPxF6AwwXGeu1rz7B2xIY9rQSjG+VZfTntC7t5VJC8oOidKRD9VSOBRW+2Qwai7ctvVpCO7UhDOJRu08qFqHB0iWm5O+RodjmpDqTIZkV3a7Sjni5i+KmzK21TYaRdrSWqStKUYqJdfiMeygG8yujIhId21TarrauruMhIiGJDp7S3ktOVW0/SFKVr49n74qEZfM9va1SV/6vsodcdcfIYiTESHtLPP6acHjJQKrKCGRDIuHIh7ya9jPi3IdXLyq2VZxMrS9sIdNv8AJO09Q1qiMfjLgvh7+N3FaCrf4cokUdVyXTZq4UZNkPpaf2qaqXGhIdIjuGJKuVTxce3QXL/PBK5VUm0iOZDKJF6SJZzKO6UdXKhaKsbcEW+GJOSiUhUq3QsFK0RjEZR1ftVynozjnQ23D2Yl/P8AMiKbMpCXaLm2pNTlLUm9IxLaKQWXSIilp0qonQ1vMRESlJNY5uIjbIlGVVMUSGVyinKV8dOrTajZdVjczgY3WqKrM5Hn1bbVBVhOiJfFUJUVBJbHVYK/ORjq1KDqMxG67dzKJfkQ6k086ItxjIuZRs9QxmFeRWiW5RmGKkG8vIhIiKJaoodmnkS0jPOW1IMUzjoiIjLdy2q25DlvC9sEU30WoBiUiGQjzRVpYysjuErRuUWKgfMWCIRFqN2rmH0lC1uVVbTTkXLSafkMtQxVvZouHdqlH0UVmbQ+xn9pcCp0xL3Ik5PKcvpVendcQ5o+4MvONx07YrmWYDJ0juH0Y/xXtXM+h2X1JCRNNkWn2sdKqXSDq1y0pRERIuwMl34yx87nve3kk8CJwdRK15a7FoZAXoyL+C7rlXVbQ2iQjq08MVbKTqyy0RG3/wCtv6FV3FSvMVe+RDERcG3TElWHm3ZkUS9IV7MLq5y3lHwCsLq3y3lHwN/Is+2R6eQKdg+GNpCpKjaiMibGUdW5eqvxcUPKMe4KUHVvl9pRH4Nv6ES2nNPOfQpsn8ypgcAuHKREQyXqanqWhERI5CLYjqEdv7lGU/QCkacEgiJD/wAtv6EYHReJe2+qP0KpsV5t62K4n81ciJRFwhGN3Kqs+RiIxEl6lrOrqkdc4jhDL/4xL+KaPqzoyKVtv/LGKradR5cp6lwdpIkBccjbHwr00fVfRx0jd/y0n8WFJtiP/bS7F0eZH6HmGX3/AMkjCnIbYlHbqXp/Dqzpox82XeFaLq1pOUSLu2p9j6vKzrRCSQGDshjLWO0l6ld6q6YuUf8AtpGPVYwJam/RHcjsXXSV6POCOVNiW1gpF/2v8l5W6RNE5UukI+6F9/yr19T5GTTBNSERII2y0qoO9VjBFLiBdIitL5MUbDzfl4kIxJKeppaR09leig6q2B0uNfBl9KQfVS2Ol1q7/lufSjsrq8540kR0xWNtl4l6Dc6rJRi616Tbn0oeo6o3C0usesP8cUtjTz84BSinCaJduLqgfuKTFva/mmR6on4yk1LvCjY04dUtFLUUor1f1SvS6P0Oq1vd3lzU+qOpIhHiNDdyy/gur9E8lcyrKxpnSFwmi26btKWWRyLm3pHurl/Xgx/6Zp3cLojL0V0LL69twREbSEfWVY6ycvKry98REScEZDLmUYHZZ9vP7DjYkRbkS7USHtKuVjNSw4TZNlId0pfwW6cnzIYtldLcruRdRNdUvyttLmGRer/goeqxq2/OEZFIuUVZss6PVlScbhHVKX0q0UvV0+VxOCQ7Zfy/Mp76VMXJK6qqYjcQ8pR/kjcpqak4jdp5iXaMv6tmruKQ+kMlM/2dy+jacPhBIB5Yy/ejvs+vlyJzKXSbFxy3m7SFp6QZjIRiMZaVY+kmfi5EeHwx2iJKpY13eiS4+Tnetw+j7Y7d06I0OXuMCQtNyGMvNlLSrTgTTAWjEd0V57yPPzpiGJFGQ2yVrznpaT7TfDIhL3S4R/iieon1UZeisy8Bekplx33KQrRISjzKkFmFY2+W5twhiMitVpoKsZOcS6QjFSQsNORskUbSIRtUW7deONk1ULSU3HGLlpFbHV6yar8icuHi2jtkMVYyZb06e6KWzSjzSjtLclY0cydyx+mIiEjltIZF/FO5TnNXSOALnFcEnLpStXSHmRiNo+FI4FMI3A2Uk4SHzbpS2LZcMXJbbUJl3SXiMSISIiLUVqlMaGmcMSFvduG3xYKZHLmhEYiPdEVYVluvcddEokIx+/51Y6JwYRKMiLdFPtUYDLT3ULUUQk5LagguYutXDwxl2RJRmORMHdECldu+TH8il36AiIolGXN/JLp6Em9USt7SRbVap6JtxI5BEfS+VRzXRsuz6Su9dSk4MZRuHVJOYZcQjEnB0jtRINqRU5DbzFpjEiQTGTGLnN6OlXnChJspcTwyUa9mDYO6SiVsrudMtq+FDUgRCOki0kp/KHHWyKdwkI7rVMng3baifNRjEZDuQQTCsKVwpOciRUz5CUfMP+HhEl5iyJCJDESuj9/2KvZm+XAdkWlpwbbtpIn2V+npEa0kBXVRSLuoMHJSitPaV6ssv0+bv2coK+JCPaU6OZ+FVlvCPKigKSLFLGNeKcxrhVfDElv2TGSXUrVhwrh3aU5jWN6RlHmkuE9dPSl+jJhpgoi4LZFcXa+qqj0N6b1blTTN8UiInbhIiu/kl1Ldensa0ZRTwPy3KtE+RFpuiJF4UXg6jqvdTxHpWYOd1QHELmSuKi4msWBrMSVd4y3i6SXQLE2SWJEqd0hzj2HTcXlluXP6Hrek5GI28xJXAbdwwJYOKovQ/pe3mYmTYiJBaV0lYcKhLonaXIh8Sy2W5Q2D2lOey1UkG0ot44qL9kJbdSKejuQ7DHspQoLGsiswrR3fFS6l2HeTurPKgfZA9rwpfsoe0l1OZixiPKq50zqeGxbbIh3KabqRIh7pfFXMesvO+FwG+YrvWU5ySO70HHOTlkpvHOCaIYlH0leKKtF+mEyISEhjHUvPeadJxkQirl1UdIOKRNEUoiRXERLDgzm9Pf8A5X+NmHH2wi/VPRqjdLiE23q5VjXRykArWmx9G1adzJuRDdL4qewftG4pRXT4fJ408FMw3cINjHltWpDK1Ck5atAJFcpuO1yluuiREMpRu2qjdYtdEWm2+bm+hXOpm2JFEbdS5T0hk/U8UnBERIojJYZ3q6PT4XLJUuluEX24jaTUvSko5nLXHSkIlG0RttVjqswpmhj7YRS1RKKiq3pNtaERXn5+a+l4/E0YzDLiYERu9ZBcaJcvKmajMjdKRFL0kPi5Jc9l2rKJzL6wRckUhGQxjH5f7leaIxK4StIeyub0DZO2jar3kIONDFweJLSS6sJ4c+eOhLzJXXEMdUkhtwh0lLsin6gSckPMo3DKHBcEhct3LWshpm47pKP3/ilDRlbKRRuSqZgh7XoooakQlIh26lMAV6RRiJeqjKOW62Mi9FYdQMrvmikE1cJWlO7VK1VE6opvCVw7k5gPpIUHIjEdqcaqZWo2R/DBKwwGOqW0k0eMij/SmTetjEbSjaSAU7iJboxJIcOW5aB4o3CmX3C4bhCIyigH3qaW7V4lH1OUtEVzY29lCUNS6RcMhjHdKSj6t2padJxwxiRW7hQEq7TkJS2yGIxTpMkRWoeizAiHzhCQiMpDFF0tc077WRW6tMvD/gjYJffMRFuIqr58wQiXLw3SLwkrdg+Mv6Uzm7LZUzpS0sP/APiJCcvpfsszFt8RJuN3dRnH1CuG9AelpUzgtOXCRFcUiXYsvqm32+K2QkJbR1Lr4ObfivK9T6frdxIYYpQkQoZt6UtvZW+LIYj4l2WuIew72tSdxcUa26ngcSORHdLejLGaCJH7YIiI2iURVcybq5apqlupIrmylHhCPyK9C5/Stk6g9HweEikPKnsCUU3UiRF2U/g8g9D8HCWcVBC6lcVOAZg6KVgaB4i2LiCqA6z6Z1+hFtoS1F7XLk/yXnh7Iq4C9rcGRcpL1XiUhiUYoWopGCEZD8VTU7U3qPyl2mYqeIJCRxISiQyt/wA10V2pEUK2QtjEbRWHiKUTodTuyGScw7yAYNOiSoaGYktYEhscO0lSip2oRisxxQ+Dq3NMdT8lvBxDYmswxTh9RQPXeiQrkHWldWNN7uGJdnSS6tJcZ60Xv/d2oy9qEVz+qusNvU/ibJzeXN64C4jmrUStfVhWcCubEoxckJS1XKOeYkRFcN2ntIqiYJp1p8drglyyXgYeq1k+29ZnhycNn/w7RUMDxLS1EixtTVEfEabc3EIl6qzE/EvoOK7kr835JrI/JF0dW0wRE4Q2jpJRwlcKo/WpVuNNNk3bdd2lpldQTaV6edK2xpHOFqIoiIkuF5lnLhf1Kbqq/iN3FstHtKpFhLVzLzubPb2/Q8fjZs3CLUSyKeFkSJFDSjGV3oxIVyXJ6kBtEnAxWnMbiit4KaVqTyhzhPtFtIhXSscJjaURlLUqBRsDw2yL+n6FYsvqRiIi7IuWa048mWd2nniut0pJv7Vv2REJWyiWokMdaItE6Q3N7eZa2smPuEIuRlKO0ex/kqDmbxMOOuGbpSGI7f4q6nmYuDIW4l9/8FCZiw2/GRd4bfW8qUNT6npI+QiMtO6RSJTnRzpW5Y24IlGUSlH+Kj8yyljSxIiEroxIf3KIbpDbdEdMS1FIZKtzXgOmYZ4REJcMRl2lJNVDbo6olqKKqNCDvDEuGRCPugiRD3Zf4oD8IVLT5NN2iQy03D/JTBpb8vrHCqya026pakc1UkL4tEIxH3SRSL0f8VVxzuLZFESdFu6NooSnze2ZODKNwyEv4flVp6uhC4V0lHVlfwuW3m3KuUle4+EhcEY9pB1jT7giRHK+67b+vFGxpZMM0Fy0RHxJtwBNuLkbuW0rUHSYNiJCIjIu7JOutFESHaNw6VPYjuVsCx5pspC7zbVp3KhFwiEiEpEMRHV9CCxbccGQui0UhISPb2E9l9bIyk2RW3FzJ7hHaNh1gi1FLmRuccQaR0huk05LwEgM2ry0sDqIhItRD9CLcqn/AGM62cY+xn5EUtXCL96cLL6QlR0fFohISlzafCp/o/nTtGPC4XFEpFKRIKor24xbKTkhkMSSmHx2+kOpZ4240cmHZZcelzcouNC12rSR9H0oYL3Rv4o/yVMrm+KP+kZKGqOj4kJRJxdPyLPDiy9Fv6djo82YcH2wPRIS+VGjWtc4rgo0D7UYulbbFbCuzACiLhFy2yT+UznoMtu+YVbfvg+qsdrmBH29sfSXAqnOswEouFItOmKZwzl/SUpbhU31TSehruWX5q2++62JCXD0x1EpocLdQrz3Q9JnKZ0SGIkWqSnS6bGWkhl3oq8fUJvorHaW8C7ywxJcdb6XVI6SEfS+hOtdLqwiGUSHvfzVT1KZ6OuvYAXiSykPpLmdT06fARtGXeUM/wBYdTKPD095VPUxPwsnZcDScDu1LmOX9Y7gj52ml6UVPZN01pnyITi1ujxEe/Kzy9JlFvcdjq3JxkSc0igMuJp/S4JDKWpTx5rQ0Y3ONyjzirmbPLi6n6fLjIZJ7DLTVad60aFu2Yd2Yp3Lus7L3yiJN/CCtNs+qafpiFMt4ju+cpmir6asEYkN20SEv4JrM8uERkMkbGtIcSuWYkmccJFIpS8K15eXSKdoFSSSIhTGJrDJEpCJWri3WRIs5CNxCLRRXZMCtLukuJ9YFfws8aIokAtjLbt3eVZ89lx8u3+O/wDUOgAkNw6fv/gkG2PDK3SMu1qWN5/QiREJDLvDqUXUdKWouCLYyiUbhXi4+mx7dntcnJzbsk/xdeyB+VJTF2BFHeQlDdDnhcyqkc3EMo6lKv1ANCRObV7XHZMXg8mO8tFt4xJU7rPY49NGUY3CUooaj6aC++TUW4yKJCQ3Ivpy+0NIJEWqMbtSnPkli8eCyxx+sGIt26hHchHG+USUjUOtCVw+sgirxHSI+kvLyvl7vDhrEtmmkMit9FLccbAYiUvSQFVXG52UOLktXxlHV0SHSO6SksvouLHagqZqV21XDL2BERIR5VOVZ5UJS5aZWiVrZaokj6SjJp/bG3al45pw4iIyIiKUR0qXpXBIZFGRR2j8i0wY5GsG3HTISiIxtQlUyQge4uGUYqXwwG0o7u0hMweFsRiNw7pfNWiUJltU+VIJE1EpEMeyiWcvEhLVcSNp2CdGUu8IiOpFg3G1MkSzkwyIikI8wxUdneRiVwkW31Vb6ohARiWpDkAOjFGw5zUZu6w37GiUdxSK5Q3slxx23UQx9FdNrsoaIh/0qDq+jbcvNykJff8AMn2Ur1biLQA2PthDI4kolhgiOI6i+KrZjkvnZFIo7SkphvL2xEXBGJRG7UjsPCuZY1whi7IR2y3LdbmYCURbKO0uIrg7lnEGRS8Kh3cukXDIbR0lH7+VTcjiEoa+RhEt3MrNWZq22MiLUMSjqH9qCLKhbEibauG4bvm4ISoyBwwk6JARRLcX+ycm0X7PvONk2RjEh1RjqQFFmJNOSlaW0kuloxYGBOSCWmMbe9it5q9TcUeHp7275E7Afp6oRd4kic4pSK2MVN41QuMVJEPuTke9ElX8kNt8iGMhEh0qyP5IAtPmMg8w+XNK0i/UjXlOU8Csae0f9KeBjw7dKJxG3TpEk00OlSo2NOIjFL9jp7ELuyW5MP1AyIZWosPZAtN825O8ANKzARjKSJxjbDlS0NgKqnbjdzII6Nort3dTWfVZCQhEoy5Uiobg0JSlP1Uah7oSqyZt30iKJf7KHrcgIfaylFSFTmosNiA3Od5A5fWC4JCZRK64i1dn86WlVEO0To8xeL5E04JjGJF3fOKzMN223DK4pJb7BORHhkQp7G1VKqdtkXxkS3nDgcpeJTVdSN2jwy03IfDJWyCQyG7SjZ7AlnLkpRH79lMjVkZkWktQx5U5mOVOCVokSBZZMCuEk5YepftcQ6VFRMREpEW2RcqomcdI6l8rjIfSJazXEpXKEJdfF5jyvW+KdOrMtRF4iW2qtwSkJF4iQ2IpYCS224JNuhdEusmpoSb3AJSK76cV6M6EdZVJmQNjIRIreHIbf24LxkIKbyDNnaN0HG9pDzD8qXZftWva+YUhOecbkQ+j8ihniiJD2lUegHWKL7QtuxEiIhuIVb3KhpyRNkrxzmmfSsZKQ9pLHFDUeHEcIRIZIzClJXE5QmXxSJcC62zEs1K6RcMdMl3upx4AuOOaRErl566ftE/mTjjYyEhEhIdK4/VZ68PT/icJ7naqniSRL5yMxy5+MuGS2GTv8pCS4ZY+m5rJPDs/Vk/xMqYj7nIS7siUl00f4VC4fZFVvqmx4TZUzm4SKRc0kL1t53ECpBIboy1SFdk5P8HzHJxW8rmlLW8N0TG64iJTnTLpD7Lbaa22kUblVAKK04creVc/a162HDLGzdkkyWhCXN6yKYoJc3hU3TbroI4JJ6lpDK6Nqm6TKU5VNixaN3NFTc/BEs0LgjtEdUiKKlaeuaKLTZSK37/kUdgBVIjdER1CpbK8rGJ2qZNsskrR0wDdzff+9FtYDKRFHsxQwt2jtEUcUYCS1kZ5FiV0UC4152MpXDt+/wCVFBUNkMhtEUAyc3LdJbiTSmWWBEbbhTuLIiMkw3xI23RTrco3apJwEONSG5RTtL51uMYx3EposYj3kCdPIpXW6UUjjNLEpFclOxErt1qzFyJRJMvP6hjaXiSBx1touWUU4NKIiMfmptpndzDcnal63s6RQCsQjqj9/wDoku0wltUc3XXELsdJENyVR5wJcQXIj3UrBsSNMIlIRiRbhWVtLxBjqTlNVi5cO21OcQi/qVyaCr1WQkW3tRUDV9HLdRXFy6V1AStQp04uF62lEoVDonloscQZSIiG6KttZEaZ8Ze4VMfgiS6bLmxKUtXhWq8R9jPiWoWKmPwRKp9oyqPN8hkIyL1vkUfV1b4jIWiLtRUq/gIylG7xJOBDEo6VmtC0eYVbjgt8Moj7oSlSoplIrdJJnFyJxERGO7s/r/Mi8XyuIdw6Y/OTBY4xbISt26dS1SUxSG60pJMy3CN13L/BENVrbYy0x2oAZzL/ADgkV0S1LVVlkhkJDGUiFZhnDZHGNxEUZEpBx846hKOq2SQVWlyBsqkXCEt1pIXPMmblIW93aVuEyIZW+iMVtxoSGRD6KZ7QGUUXCbERGVsu6iazE4WjEhK6Kk8GBaK0oyEu1FLpKASuIpeHUp0cqm1LL7oyESlHSVqcy99rhxIxlyy3KwV7TYlEdRCRcyhmKMboSEiu+/8AgjSj7ES5f9Khs5JuUeFdpkKJzOldLhxcIY6o9lSNFQNk2JOFIxG63cg3OK/Ln3HCi2RIBzJ3RK4C5l006MeJEZF2ko8sbcltKMVc5LI5uTimd8uYUmVOObYipZrJWxuJwe7uJW/HK+EOmPa1KKraMi0kWnlReW1pjwYxGfglsxKI+ElF1OX8MiEpdlWvK6Fxsd121bzCi4hFbcs/dsbezgquX1LjDgk2RapErvSdLnyaEeJGPLH4qqtRQE3K3amKY42q7ybhX02Dr3QHppwyi64O7UK6Xl+etVMibcGXKK8wMPk3cJFJWDJuk9SxpMpbrlXH6i4sOf0eNnh3bpZmIsUL5lcXDivPz2Zm4RFKNxW8qmc86ZO1bBUxiN2opFpVRko5+T3Ltp6P03tpxqpcIdReqnsSIi1bVHU7zcRu+Ki2XBItS5fLuyTnRtwmnWykQiJXKtdMKg6uudIRkMu7apZutECHTErfv/im3aikEiKWqVsVvjldOXPi3ltXafLzLaSMp8vjqtT9ZnTQ2tCJd61QlRm7hW7e8iY2uiXUTWLrDQ7Skgn8z5Ij6KhReIt2lONYjuT66ZZZeUm3VOlpIi7oil4AZFKJImjJhsbR1bpSL9yNpnGy0lIuaJLKwbuiMoZInNJRV1YZFsRESlzKAp24xtJTAvDaJfGRImiaelHVLtdlNvYx1bluqrybERFsYkofMakiG0roraVGkm9FsCiQxK2Uhj+5Ch7XIZR7KgsmrCbB1pwiKeknLol8iPy9w9JERCOnaPyo2LE7ltYPDjKPZS6uubGN+4bSUW4x5uQlGNygKmTjoi4ZCMtQiJfKjadLm7mbEhGQyt1JwHxLTp8KrLGXNOlHilIC1RHT9KmgoyDzYlbHVuRshpgUpFyyTA4iVzhJbuLhCIyiI26bopbdEQiMilFMiMHiESjcO1IqKiIAUZFKUY7f2rKkhbbIi0ih3czpCaEiiNttyAZJoi4hOadsbvv+RBUeVcSLktRdkk5VVzbhSbck2O0tqLoKxsRLzgtg3pLaRfrQBtM0TVsSEdyjqusdlYMtVtycqqsiblxBLtIamdEW7olzST2SWoa1yF2rcNqNbeIhly8oqHpAldLVtRwVUbY3drlRKWhLThOWl6v8kLmBeaqbS9of5rfNEn28I3SESLchswOLD8i1MVMbdRcIleKaQyQvi5LbpT1C0I23EMlVejebk4QtuWuXauVWRtwmy2xLsioaFVTbZEUZD85QbuZEJE0IyJvcn85rXytbGUo7dqbZYbFvjmRC6RRjtSESbXEcj3dKIdoJCKboCEbS0kJRIVuqzAWmyIiLuxQYXDKxbdGQyjdIf5ImtfIRi2Po7lHtdIuKNsuyMSFIqswfBxgyEpOFErZWxQcTNNVkLYytLlTlTXW2iN3a0qDqXyYudkREUhttH9ih6ytGRE4RREuVAWIcwbEpO3DdHvfqS6PN24uXDaQxVTrs0aICiQjGMtSrrFY6RlBwhl2dSZ9XTnTbdkRbbRt1JFKTDhWxkIxLvKAyqoMmhbcHzg+siadlxsi2kkek04BEPZElomSFu26Q3Dyp2kwH3QroiSewuIuHdt9FI0d7XIiuIhHxJpqUhctiRDuRNczLUWlC07VyD0kqjAXO6g3aRuQ2lEd1yk2BEoiWlPVQiXDiWkruyp0AR0DdtqAzPLiGJDqU1iQjK5Bk6MpakukKWq3U5YRELfMoiryQhMhEdIq011SPE0xEU3CQkXZU6sbS7UyooibihDEhVtp6LiOzjaNqFzfKbiiPatSl/Wn2rRuCmcTRLtGYlpJDHTkP5YktZpUum+NpTrdaQoMsEjyYqusTchTtaRbkyTxEmcQJKAC5VXWRl3rZEkDgRfmRlLTkX1lLt5YPDHUJSuSuUiLki6WgcLbbFSlPQCIXRIi2yUi1TtxEW5SlEiiSc/BlwkIyisrlsvvyDayyV0hGIyG5E1FaIR4QjIdURtJaqmCtjJMOgTESuGW3akpZcvdJ9sTIYly6YlFakJEUiiqw/nzgNcNsoiWrT8ZaoK4yG64h091Gv9hYOkFaY0wxKRXd5RNI8QjJ0olylqRrkiaKXKWlDZfk5u8N0pRHtDqVbKeRlGLcboylbItqlKUdoiPeFV+kyt1x1wjIhFsSIdqYyTOHW3RA7pSG0bkbKrRVYah9ZMY0AvyG63d2lJs4CVzlpEOkrfvikey225DqISRPJGclyXgEThF3fvgptoJCUh5o7kBS5kRFpt3I5uolIto26U2eRQsiPZFJceEtUYjqKSiMyzQolcq9V1rpMFErSLl+cns4suY1VMTRSIe7JVLpAyxwh4erVaUiQOOHFGROjKPd9BRr9WQlE7o6bk9noTlLxNleW3cn6qviLgyIpEJaf5oNs+IIlq7KewoScERG0t0hikNJbJcwYMYuEVukS3KdpWmykXufLpiqp+CDbEXCLT6qOo85IfNFKPNEkFpZMKkW9sRHTuSal4SIHJR7SreZ5p5uLZbu0iqDM2ybESdGURGJIiU7WvuWi0UnPVQ+fPuFTSGROcMhIdojC5IZcuGJCXdJF1pRYflq4T8vgiVS+U5TwpuSYFxRcckV0RIoq2i7EolGMfv+ZQlLTk2IiQ3SlqGPiTtbWCLgjuIe8Mkl6TwmRRiLcea5AZhiUYkIkIlaQ/OQFFnQxcEi06bU0OciREMhOPZSVjE9R1QiMCLUO7b3UDmlZLzRRIdMi1Ku5nmIyEhLlTr1S46I8NsiLcQjJBp1tjhA25EYxS6Or4ok6RSFoto2/fyKOp3HSai7btiVsf2JLNI6LXDaIhkUtQxL9yAedqPZlSQXcIR9b6FuvytuMbiKUu8icuYFpseJ7YRXFzfsRRYMSIiduiURkmNKfW0IkJRG6Q2ltRn9nuBEuIMijEZKZby0ikd0S5k3mFE6UYlIe19/zIMVlFCWpzV6pCmKusEXf9PKiPZBMNjLVHvF2lCvVok5ErZbo3KdmnKdwrT08RStOQtkIuFEdUiiKgG6xtgREyERHmGSYzDMyftaIYj6KexpLV9YNxiXmxIvv/0WNvNi1xBIpFzaVGtZZUuNw4YxK6RI0sucFq4hiPzUpQkqWr4mqPo/zTkxESu1F2VXG33RHTqLu2ogHXBGToxlpuH5EHEw9qjK2N0v5Jkm24iVpXW7h/eozNqkhalxBldKMVEZTnTZEIuORESLs3bf+qD0slXTiRDb6QxFMYEIjwyUhR1LTlolIuaKHqLZDEdXetUU41SttxIRlLVKSZJqJERXWx7KUzU8O4rRT2Fa2RSGN33/AL0j2CxZbcES4Qy9VR9VlwlbD4qmHajdaI7UI2ThERXRLchW1aqsnIfc0ikyWe0hU44bsinaO0huReW4iZcMSGUdxCKfaixXs3yttrhiI6h1XEUk5S0LYtW3FpKSsvBbIuG4QEQ+kn8KJspcO4tyNosQlPlLbbYx4npRinSpBuHl7qmMWI2x3eqmK0iEZC3K6MiGKRSI6iauKN0vVU3S0lpDEfnJrLKUikQjp1KWYARtLcjRq7VNC24IkOrmQFdhISEtMojzKd6QMSJsh5SVRfdiUZbpR1IBmppGiEW4lIR2l8ZSGU5eMbdpbtRLeL7WkoiW77/3ogSERk3pIrSFVN6LQp0W4EUiAREuUdqGyXMHDaIRK0bpFqWjAnWybl6w6fpSX8rJprzGkhuKSegPrcyEaYikPEISGUSt7f8A1VeyaicI23xuiUZFpTr1U0VNwDEeLcIlH537VIZDiXCIYxEYyK75Uh4JriqePcZCI7RK3wpprM2AlxSiRFuReaYjwC85p5ZESpwUbrpebkV1pEnilcDzsBGLcXHS0iMRu/b+VTOUG+bY8UeHK6PKqAxQVLbg2jzDu/2Uq10iqaRwROMBlpul/NUmxPZ5TjdErtIj/qVbq8H22+HGQ6bbh8SVmWfFWebaERlG4h3J/K6N9wYmXduEbv8AZGtFA+X5RxBJ0ijEpRttaQONCL7pF7nK2PKpZ2gdbJ0SfkJCVolt/UougfhaUR4dumX+6FJLLsqEbSEo9pSY0bQyK4e7FA02ZWj/AKlL4vNOCJFbJFAb2I262WrtbfvgmxyYSGN2nVbJSjLIDKJXF9/70S1aN1opJqtfgoRIhISIdsub9S1V9FxiJDxJEUiES/kpb2U2LpERDHtCWpPeyHRcGQlw7YxlFNKPoKTgRE5R2kWpEuGT7D5XWtO6d3miRWZME6LRN3COqKSJcNh0SHUw/wD+Ik59lYqblSRXaoy3INnMylFyNsoyHSmKcyES4mr76vIhXmiItO1JpThPNkRERFLUMbRkmHXSIpDbbqEtST7DdIpQKJbo2p+lpo3FIo8sVXiFDLTD7twjKQq6ZGy3RsScIhKOooxl9CHpItNiWmQ3SEdK3T48eTZFxJaR1R7SVp62Qw+4+ZkRCQSKMVM0hCLYluHVyqNpmSaGJDERkUo3JujqyJxwSIiH1UjkHVDhOy22yTdKIyk5dLSKDfrBbIhErVEs5k5Io7e8imuTdYRFHS2OkURgIkQyLSJWqtUtW+7G2OnbFS+Z1bbBNyltu5vpQWialyRXXCO61AVuDfEAgFTD8SESbGQldK2P38qbBxgYi6TYlulwxU6FqCzCjddiLctUtyHy6lJoiFwtJaZKyuvtDcJj2SEhTDQtx4hE2RFzRIlPbQmSWxfkLYiW0e8kug6VgkUdw8qDbqWpBGNsRR2MhkV0i09pPatkZiLbcRK2Q2qE6Svxai2V0Y2o7OKR0okfEIRHmjEUANIJiUiJwtsiknKast1zul2RD2kTl+XiLrRGNpFzfNUy7lgkQk5dH7/rUoVDTEOoScEZDdGKLTNZxWexICIytlaO1RtLm5PuSICEdxS5kjOHx4dxScG24tqbo8Ci3aJERCRR5UgnXsfNWxKRDpTlO2QiIiMkyLxFpEbdsdyLoyIhu1EXdH7+RIEVFCToiMo90opzCkcFrUQ6RG7SXOiqh4WhGZalvGuYiIyu5UaMDhQlEhcKRS28q3S5aI3CPnOzqTj9aOkSu9H5EUxUCXDISHtdpLR7R+GTuNuDUkTmrcledEijaPo3Cic2rhbbtKMS5lCDnhcIo6iK2WkhkmSxsecbIiEhc3IauCQCIykJXRWma9txiRODIhu4dtyj2M1AiNsSIiG3vFzI0SRy6qJoXBtlEfuSOpSIRIyjIlCtlIhkWgZF2u8ttZkT5EI2tj6IoA6tqmyKJKn1VLF3iah02qdpmBd4o8TSURuuinmsubjEpFdqFBItikacKZRuiI6Y/wC6KfohbER3bRSwbEXCGO4U5WARDIRuHvSQas0lfw3HRcLu8yFr87dK1sojpikV7oG4ZFqjG4R2/wDRC0FCTjgW2EVxbVpComheEzEam5sSHbH+Cu9WQ8PzQkIi33RUR+Dx2i3Ee7L0lKlVk2NvLpiJS9FTslRrakhO6W0StlFNZpVkBALUhiPLEkXmFTJxsibkQl5wY/G/zUfUNlxSIh4crhlbEU4BFLmbpFFwe0JR0pWFKLguuOERFK3spll0ZCNokVspKRKicGUSttlyoJGNALDm7m0pfsx8ok2RXFpuRAUgTk7d4pR+hJiIuFEIgJDGSBYk8kYdN8RfIol6yMeycW3TccHzcpWju+lRo58LZkI6uaKJezg3TBspXXR26f4ppSJNtxGz1Uh1sdqJaIStIhjGSAdcIrtIpA0wbolIikIop3Oxlw4pWA+bkcdvLpUYTLbhETbYyEtv80ELffFwZRHmu0qKLpE6NhQIZd6P708YOiNol4UHjTXDxWxER5tSBpbuj1cLjFwiMttwqv5vnBtE6FsSF0bhkVwxtQGGbAwcm7pahiSjc3ruOUxttKXeVyeSyngbUiMiEiGPLpTbjwjHhjpHlld+tacyqpdInNu4pI2j6POEJSG4e0ouooHS1BuEIjKXoxRboRKO4tojpVgy/KxabtERIhuKMkzV5ddJZZZfhopyojESGQiO7chmq9tohIZCUtvKpl7IXHRGJW26lpzotLsx3CW5Vjl+mBrc44jbgyK4bSUc1IbZEJRlJSRdGXJEIlHvFK1Cnk5tFcREWmIqxtG1b0S1ERfOQDTpCVtslMfgV9yRCNvoyTTGSOmURG4S3Wq5qQkzklURsFIrhiKdzMW3Sk7cA6RkswyOrFqIiI80UNS5W75xtwpesShawfhBomBFoYxERHtKGcFw5C433S7KAr8rdYGQuadoptjOyG0ht9FPSanvwX5qUuXd9P5kqmYciIj3dQxkm2+kIttebbKRDEpDtULU5m+6Xm5S7JRU6Gkk2RNuuW8SMpdkh+RNN9IHDKMiHvKKD2SV0ikWq7UpChyghIXC7Pa+X8iPB6STvSEjbjqIpCW1ZlNWIiQEJS23D9CRR4OzIeAEdumQ/fBLN8WiGTd13KlDGk6LjkYkMtUopuooicck04I2xKXMgsKwnXB23R5VYDEGxIglON2q5ByKrU5Y6LjZOE2Qlbb9/wAqk2KSIyFwYiMbfmoCqrHRK4S1bij/ABTZVzpNEI2yJBjuI63cJW7vvinyzArZW6dKhKN5+JadW5TAgRMXxEpcu79qAdzHz7YuCZSEruZR+YMO8VuNokI3cS2UfpWMtEJCIlISu5bljucuCUYy4VulAZFxspEVw6rvi/kS6bM+AUnCIhO0QHUJcxf5KWy95ipYkRXaZFao5/Lm5CRQIhKQjqkkZVFgVSRFEhAi0lcj66gbIeE22PEG1AC+bbgiLcRK61E45gYucUhLdyp7Ko38B1heaEiEm9QytTdA6LBEJiMhKMlJDn7huRaiLjno2/InMsoi47jjoibbhdkruZBFtOyciI2kJalGsGLUhcIh+LH9SkHqgW3XCPSMuHElUcyqyfcKOmWnT9/yo0pY2XeI4PCKIlaXiU0BONEUrh291c9ZqibK0o6e1crnleZibBC6I8QdNupGk0TgYuEXCIhLUUhuko1/MhFsikUhkJCm8czFgmyjKVpWqB6QviTkmxERK4k5PIbi0Qk4Y80U7lT8uG3pECIrRuLs/kWZbUN+xjEhuIi7yhHXCEijaPZ5VfXYXPKa0SErdWki0yS6/F0SErSkNpCQqm0jp6RLtcqlMc4dgNxFHmU3FKx5NRsCDhO+dddlcV0S2qs5vTG05uKRW3SRFFmrnEEikPdRVbXNuOCJbeylAiKejLiSIStJTNdWcMRiUQIbpcyKx4fClw7SuH74KsZtUk4URtEdqr7IfldSLjlxSiVsuVF1cdqrYSEZCWkvWTx1rmmVu5PRUVTi2LpE4QxlyipPiN8QHBH0tKrmDki3eJYLhS7MkaJfgujt7RbuymHqwRKMbVXGa522JcqmKcBdKRd65IJEzEmokQyKMRJJp6LhecLdyqHraZz2SIiVolJWalbcKIlGO1Kp2bcakMpDHuqFrGeIJSKPL3VaSAoxiMUDhl4lLmigbUwcicK6VveSX8sEQIuUS3CW1WippyG1sh5UzUU0WnJe9ul6pK8b5LKrG64MrRtLsrMQG4Rjan8GRK71tyQ4z/VzLHKnCMEksVjgkK1i0RXRUyKOYVAjbp+skG7LcXiWgpvWRAU4j6SvRbOM3DpW8KJuUjbEh5rZJxhoRJFGqKVH+xWxEo7tMbVunpxEpfNTpEOmIp0RQq0jgCW3takrFhv0ttqXhilGQoEyA1+WtOhG2XLElXHOjgy0jy6VcRwWYEI7dqDV0OjzZahGPdRI5U2A2xlHlUjK7lTuONqiqiADKm9Vsk8xTi5tFS4YW95MutkJDH6qVNFxISIR8SDeopCPFIS26VYCiMiIZd1JAmnBiOoht2/7JyhSqml9jCTojJsS7qAfzdxz2yUY2CJDqVsrsncOIuCPDEhIhlqjtQdRkIkdoiIjpERtFPZ7VlrMitEhKKynJ/iCTZWiYy0l/FWPHJxEh82NvZ0rdS0I2jaJXW7kW6DKp9oYkUdMtNv7lHZvm7RNELUpFaOoRTFVlN0iIyHl/wB8UC/lwjcBR0xEtSJTS/RvE4kLvnAIZDzSQtaJcRwhG2RFp2pOWE6Lg3EVpCQp562W64kbOI4KkgK0oxuIRtUzl1YLoxIrtqbfZbiMh1CNyZo6Rsi9K0dKRpCsxITEmylZq3J/LXjfIG5DEtQko6oYclIit06tqapR4TgkOoSkMrpJEtGe5dTF7U3FxuIyG2Rbv1KGoq0RImilaVylMpB2pIi0iRSl8xRfSDInG3eK3G8okmV8Gs5zAWyhwrS7ur9ihsKbilYMbpEpbP6HhONXCREO1IGkJhonClcNvZVQtoiupBGRCUo2+kgctxk4IlLl1RRXswOBEriUZjjK7tLTELTncRFsRutu7Kq1Qd2r5yJZfIZSIiEuZDwuRPBMafIY8qW4+TkiLl2oulpxKJEMRSH6YWyt0kjtAHoiiVyN9huEQkIy8KaYYAtW0lbMmMY8ONqnKhE09K4LYyb901aoqVECIhG3VyiKk6kmB83plH0iW8QbbkIjpHUVxKBaiK/C2263Ttiqw7TXco8qttPEicIiIdvZVeraYZOkMpSJVEtvNNsUwiQyIilpQODAu6Yy7SHddMoiUrbRFGUNM6JDaQ7o7lV8QDMoy7zkXIp7OMti3IRuIkRT1PCujIrbbSTr58fs3SiplNFUeXibjYy73eTdRSPsHujtiVyk2aEWylKRS7SfaBx0huHl7qNlW8hN110iICjHUUbVYiZ97u7pJ9ukIW4y27bU2wyTZaiJG0aKHCI7vjIWrqOFcSNqCGMh5kJiyJWl3kFoxT5q2W26W5LzFoeE6XM07/4iW6hhoRGIxkJaVHOPSF1sSIvNOFdbbwiVY/abit7DVur1ljjHa9GSGnbuTYuXLMz5U/MQpIuRt5Uh2RaUtsPEno9l+VONpgjuitsmntNGgNy1VYR09nlWA4m3XJFpQcNgJSReA2pDRCn8UxszEk5hgP3uSTJbHFAZihqnEh0onAkPUFcPdSPYYCLcSJBxLFkUnBjteqkqUyTpJ3AiLcl8IUrhI0vYeqbIhIdMhig6WihG6XeipPHArk3FIi2Qlb/SnBbFNBilY6kaONu04lyxURW08bYjbptU5hpUfUEPeRYe0IDBOCRFzaU7gI6iGPZijTbKPL3U3gyfZU1SOfYArhbES7v0oU6Ps6lO4N7Y/f8AUtE3cNqW1RUqpktMtKVTU8buJpIZKeqqSUrVAVIEEhtS2eh7lYREI2x3ERRl+1M1OXjc58XtJlghIm2yukXKPL/mpXHHzgCOjdtVwj+TYcICkMRIdVqKzWvb4QkEijqlKIoPMHR/7Yjqlu7qjX3It3FIS0j9/wA6RWC8XhqZGTftYlES+aozpEbrjTYt9mQo7LcBKUrRiPZ/imaurbcLzYx4ZeL9qqUtKYeXuyEY6tIp+koCIiE7YqwhmLbRiTo6rh++CeewEyIrRb1dpV3thaV5yjLUI28qcao7RkMbhiRKVbwu4YjIR0ykPpfnUrV5bxBbuHulaP71FyVpA1jURERtlyonHLx4XMRDqjpJIIx4gjH3RFO1UbR0o2QWly62RRt7yfJsWyEhlFB19eXDiMu19/InKcjdaIhcjGNttyYGsm2TjfEGQiVpEVqX0mzARAeEMfOXRLaovF8iG/aQj3f2JvNnBiFuktwkKJ9oDMZiVxbSJEMuTErbvEh6k2jG60pStG2KVTk1FwgIpCMhFWbKCjbcIuIURG6Q3XKRqGXRNspRbjGVopTzDg0w1cR4ZRkPa5k1WZyxwxEW+IXDEZFtU+aTdQy6RlEJCPLzfQpLLKXSRCMtw3CQqP6N1RQc5iKXLuVkYASu3eJING03IRjcnqOpa02CUtRXKGqcwjU8MRLVqRANSckXzfkQFhCpkMYj4rkM88kUzkhIR+LFJNohK7UmDToEReknakOGMtRREoiSRhqEU+40fo7kJyBPZkJC3IY26iiIkSHfaHhvlESk07dqiPCL/BSr9M1GJDLlt0qLqq1tph9qMZNOiPatL96rH7Z1OODpTXDulci3orY4lFToNNGK07jFJI5JoTIijuQZomylIkQGCWDXMsLCIpWlI3gXxUjipnEiTzTZFFG2khbJ3IjBwkzwCHlSscO0n2k+2HJy443Vy6nfKlB2UhoR5hWCXdS7T9KcuH9p/wBw/hhEUyWqSSTnaHxJOJJ7ipyY36yh7B1Zi7JDFhJJxEvqob6FEt4Ykh8ESKoNJGOKdcHdJDjjJAbHBKx7SzydpMv1jbQlIvVRsQ7xCHTpQhDcJKN4hP8Ateke1FSdG0QxkpuUbahRFtTFQZRiKKdxTGOKyuQ0DpwKRERIpsrZS095axbSY7VG1RhSc5YqvZ7RlISEe1LmVh4Ue6scZEh1S7yNnpVadshiRCNtyIcqScKIqVfpx4Zjujaolqk4YiRDKIqpQFxfi25KRFIeVOu4iIyJBVrjZkLcYiOqMrUy8yZaZRVbhCcM7CJDG7tKPazHhlplK7upY5WW6PoxJPDSiPo8win3iATtWLo3Dut730JFwkJSKIxtuUvT5SQlKIxIvVU6zl7ZagEu9FFzCv1eYOSkAxEhiVpfIhzzh+0ZS27lZ84baYbblGJFbwxkQqLcbbIbLSLcQxTl2EKOJNkMtxFFGNjw4nIe0mCpSmIxulutFaecuJsbfrJELYwF8SHT4flUR5TadER5rvvgn2RJvml2U8zTxLicu0tyueAaqqmIx3SG7ak1L5EMtUkupxJwhKMY6kQZC42MR3XClske+0MQkUiIVmXNSKWkR1DzJ+sbGQFtHVbqR1fVtexOGLdxEMSj9CsCcwzYCo2mxGMbY+lyqrjdL0k81TERCJeipinyjuxuS3pOhORU48MeYlMY5gLcmxHTuiX0JDFIItCXKP3/ADKCdaIiIpRuIvRUbLSRYAXDJwi9VSYuiJXXSjG0lAs5O6PnPFd9C21ROC5xCujHUSZrXTMk3dJP44S2qvVGYkJNiREI7vvgjanMhIR4UijtjGSCHCHnOVEOHbEhl9VVahzDiEQ6oyLUUZI0c8bG1wSlGNskFUyLHDES1CQkonOMoadbcdjpafL1SWMZs1cJEVpR0p3MMzb4DsSKJNOCNpbhJXj9xll9JrAB0kUvCnMMNqZZw9JLLEZXeFSZuF2pJwC7vEicBHuprhDLUgm8MPvJaIVo2i2plxspKaez44CnsMUK0BS1IrDGKR7azHCTUdMiG4lEYt9pTwtk+QtNiJERWiXdR1P0dDh8WqdpaYSIhHju8CRDsHy/3rzPX8nXOR8P/wAmyy97CY/io8NNvNK3HkuVSiWZ0wf8wKmlfb/jhjgiT6PULTZOfhXL3WhtkFSzcXJHy6vIuPvlPOr/ANPmv/G1vy5u+2jOjAedc7nzv81bavo/lpShmtDHtVNPL+Kj6egYpnS4VSxUkQXcJwXYjzF5F2+m5LcpLK9n+G5M/k8cuxGLPa+Kti2nMCS/IvWfoxIgtGHaL1UsFs9XMg6ZPG0hQzIxRpghyw1CMUwa8lpXXIN9niDdtTgU5TlJEm2pyOItoOHpHVzKRHEiH+pJwZEkQ0AiMZEsbFygSAhK5acLso3EfSSIdn1UtNIBmSVjIkQ5h2UzjiQqVkYtkt4B2k5g4XZWzO3akZnVb8ZNvDbH4qUtJbLSIOhk7KPpCllTiMiEZFp3KVbHmituNDdd81PY0hKNkhK4bSHtKQClEtvqoyYjaI+JIwciSmxNhLIC3t226k7hcWlJF9ONODzJp0iMzy8nRttES8SZxoYgIy9X7/lU+RCkRHcrlVpCv05ORjISHuqFq6eR3D3oq3uBIUEdIMt3hVbQqFSJCMbrbttwrZOnGQ6R9JW38HiZcT1YrTeUt3CIld4VXYKydcZNR4fpbiQzDjjfNdzK6jlQbhitO5O32i7qNjUVZ/HiRiRRjzDqSacSIhlpH1lYyywR0y5dNqbwpx0pyppuioWyJuWovCpbGkJuIx1d5CBQEJSlHajAEoxIrtSaRVTTiTUbVEjl0R+ttR+BRHctOYi4MdMkGigfIbR+NqSWHCIStEZbil98Vs6ARkXEu5SJMv0/DbG+661EJvOuELcrpej8qh8tqxHSlV+J1LYiAxESu7SAbpX2ytVSEncBAS4sRa3FEYy/aiSbYdk6NwxjG3VzfkUHUU9S4IiQlzblM9HKMmxIS72pBVD+w3OLIbh7IyUpno8OmG7aXLyqYHHhlGIlqQWZ0guNuFKUWnbeWxPG+YzyWoTj/pTLuN0opb46YpvEVIKAxiNqUBXJrAU7whinojmGK3hgK0Ip/ARijRmcBW8W06lRQBeRNSrKYR3OfNJAf8RIcPLMv7VY/wD+AVJZD/8AmMf/ACfMJWvpV0Upc8pWGKo3wFh0ng4GLYkREMYlLDH8i8b1XqMOH1nFln/5ZK+O/m+XHj9ZxXP8rgfRutBro9WE7lxVOBv8EauTeAsOkAxl+SW3FFdSuWYvFmFTjiB4Cz7E4VpOiThCXsiOP5hw4eOHl/xLyLr2TdEKPL6Z3L2hdfpnycmNRwyIpBH+7Dyfk/IXooXIOgWX5K/7Mpn6wnYk3wjJgmHhIdJDgHlxHDHyY94U+X+W4csOfHHeNv1f3/8AHkcn8jxZY8uMvW36rh3Wlhi3nVSP5BjwAtERGQsNiVuH5PL5fKuxZ9lrTBsEACHEpGBKICIl5oStLD86RnHVvl9dVcd92q9kulxXRDFkW2h/Xh5fL5EZ0hfIxG4iaB0m2pcrdOI/wgieu4+e+nxwv+WP3/1HX/Hesw5eb0+ON/yx+0QOMdqXg4m8MFvy91ey+/OTiK21juSRLVpuS8ECmn6iO1C8a6UUus1JjFK1WJ3A7dMZLfliKawSscdqVpxsCTk0xgW1LxxUrkO4HbKJLBMdSZcxt9JD4ulplaotXoU6YkXzu0mSAUtl0YjctkMlCyMA7SbJq7V4RtRGDJc0STDsuaSSjeLfaW8Wu0sx73hW8S7JJFtvCn7SQ5T9pZxSTmGPEtuSgDYtjpScQT7jabxFUZrAUrH0U+2Eh2peFP3UyDinRwWEN2mSR5OylsqfKPKl4Rut9VMYkl4FK1XInRIEMS0pbQjqW8WdK2YxuTqdlxWcMZJImsF3mRoq2VPbGSjKinJsxG4rZd1SuDwoasOQkQ6k5SC1DxFaMhKKQ04Q2kMilqHSkOPRGREMkzhUkWkYq0jcT5vCo6ozRtsiHlLmu/3TjrxGQt8PTu3aUzWZY2TbjhEPEkPLJAFUzjcuIWnaNsiUfUsi6ZFKIbR1fKhWW3NJFLvff86O9jlHuilABEm2xIRIZW2pVO6RREh9IUw1QuOGJlddpEUbm4+xmvNlKQ7SG3veT8yrYbfzIScFsboxuUrGIi7HVaqRkoydIi3aVZmcwu4ZaeblQmwRWEIxIiIZeFA1Alw3b5Dw3C9UkdWUjbojcNt1xDcgatqInGNrTunTpJOfacp4WlzGS35VnkSmhuSSwRSsRT+AxSccFcBxrBOYprEloSQWyztW+JakuFLakHGI9rakaXyPD/1lMO0nPmkuiUbRFbcPaiS5LIoyHSOlOC+RWkXrEvL9b/Gz1Oczt11fNfzP8Nn6zlw5MbrrNOy1jODQTiUiiJYCN0u1/hgoioYICF1y50tAco/IuZE6W0il3iWG6XMXiXLn/C45f+55nJ/xbky8zKOhHSkMhEZOu6yjpHlJVrpY3HgCOgOK2JcxbiVe4hS3eIlmJq/S/wATODkmfbem38d/xzl9N6jDlyynXEvDDupJjas8topw9K9p9mQGKcwxTOGP3kkuupgjjCRLToR3LfAFbxEVNXCASCwuRMRjah3AKX9KSm22pJeDUVrBtYOEVNqpTD+FyaK1GWkmnGxlap02CeS5FtOEIimsQiRJ0OWSk23CMk2bRIkcB5iS4paCPwwJbxcLSjsQTLrUkaK0imak4Mk8LEf6VpkSEvRRjQkQp3FOwDwEgzUw60o19lTpcpkCinwdGKZglYYo0Vp8buytlj9xSME7gtJEZUjAE5hFKwaWQtRUGje07uZDVOZtCN2kdVpfQiTbQr2UtOykTl2q6IihNFUb7boybuFYdPIrVlDSNtDw29KKMREZFajake82Q8qGckI80kY7UbbUFVHp7qsILMHiF0bdQ9lO4EXDl8VJq5EekrVo8CERES7W5MhNC8Uh3bdUk/mGIuW6dOm3+CjssJ/jiMR3XR+lSr1GJSKREfLbFAVphwxcjcV24pKxg7Z2rUOzl5CV0R1apJ9qlKQkJae8QoJuR8JyOqNsR3KDpaV/hGJiRcRWOnfi7ErbpStTtSAkXDEtRJhRjYKmuHuqPCqKYykQyV6JkSImyG0dMoqBzLKhlIRL1UShuudIeEQCUdV1ykHH4sHzE278VBccm2AEolG0RKXyIcKgiaOXK7HVyon2mukYCtFglN4+ikP2kqYleXmSiIYofDV3VsnIpg8IrWAfe1M4mnQJBnMMdKYwwkRd6SU6RJA3eklQfwPUPMtYY3JOGNqTNJRw4kkEkzkscw0ig5G5cqS6cdSRhhciOHxB1DakZAOSj2U6HLJBt1DcybEht5UTiJfeKYOYtikQuWxxSsbkJZiIyWOxWxxWFjdpSXs0IrMATmIpl61TSlaKUkkkvAZdlIdwSafRylFPuNW/1JhnG30kW3hamvYDELitSeGjSwl/SmTC0lFPZjHvJeJDFMONlHupo8LUtDYyVtpSSTIhQrJFaiMUxttp1Etux3Ch8QSoXI2Wz4uDdIhTT+AlpJNvNxTTeCWjmRl3UttrbzREkYtkke9nScTmLqaitY4KtpFi8tYmmW4peJI2VOYkkYmIpJYoZ1FpCiqCTNSZGMe0hwJLxxS2ZLLBDcnsSHSXdimeIkcGRS1Ft7KNhlRTjLmEtqjswkLkY2jHlUuGBRimzp5K5kQBq7sEiwwiQxIS3Fama2lKI3fFFCsNOXXR5kwlMXBld2lEvVZNjG66UU6/U8O23l9JRbhOFqGXooTsaLQyEpWrWNYTdzd11sUUy95vTdaIiofMSIZEJR9FMxdJXjOTlvZJFykUmxkJFpUa1SiYyJz1YqVYNsR1S0jyoSX+DhcjMYj8VQua0otCYjylHwqSr3SGMdyh81YcdAiutEpR5YJ4/YROHWVX+80fgqPtEk+sivL3Kk8D32ixYvU9rH8fPe/n+tYdY1d71SeB77Rax6xa73qk8D32ixYj2sfwe/n+s/GJXe9Uvge+0Sx6x67D3Cj+De+0W1iPax/B8jP9ax6yK/3ql8D32iwOsivH3Kl8D32i0sS9rH8HyM/0rHrJr/eaP4N77RIx6xa73ql8D32ixYj2sfwfI5P1mHWLXe9Uvge+0W/xjV3vNL4HvtFixHtY/g+RyfrWPWNXe80vge+0TeHWBW3eap7uy99dYsR7WP4Pkcn6Fp+mNS2Uxappf5i9j89SH4yK73mk+De+0WliPax/D+RyfrePWRXe80vgqPtFr8Y9d71S/BvfaLFiPax/B8jP9bHrIrsPcqXwPfaJWPWVX+80fge+0WLEe1j+D5Gf6z8Zlf7xR/B1H2i1j1k13vFH8HUfaLFifs4fg+Rn+tfjKrveKPwVH2iTj1i13vVL4HvtFixHs4fh/J5f71sesauw9xpfA99onh6zswwwjwKPwVH2ixYj2cPwvlcv9mses7MPeKPwVH2iRj1k13vFJ8G99osWJezh+K+Vy/3pH4xa73mk+De+0WvxiVvvVL4HvtFixHs4fg+Vy/3rWHWFWe80vwb32i1j1hV3vVN4HfrrFiPZw/B8rl/vS8Osau96pfA99olYdZFd71S/BvfaLSxHs4fhfJ5f71hdY9dj7lS/BvfaJOHWLXSlwqXwPfaLFiPZw/D+Vy/3rPxiV3vVL4HvtFr8Ylb71S+B77RbWI9nD8HyuX+9Jx6wq33ml8D311mPWDWe9Uvge+usWI9jD8HyuX+9a/GDWe9Uvge+ulYdYld71S/BvfXWliPZw/C+Vy/3rf4xK33ql8D32iTj1g1nvVL4HvrrFiPZw/B8rl/vWv7fVnvVN4HfrrC6fVnvVN4HfrrFiXsYfh/K5f71r+3tZ71TeB7662PT6sH3Km8Dv11ixHsYfg+Vy/3rf4waz3ql8D311vHrBrPeqXwPfXWLE/Zw/B8rl/vTZ9PawhjwqbwvfXTGPTOq5GPC59daWJ+zh+F8rl/vSD6WVBY+UgZx/L5fzOfWTp9M6ohjw2PA55fjLFiPYw/B8nk/TYdLqkRjw2PAX0pT3S+oMYkzTl2ovS+OtLEezh+F8rk/sZw6T1H9wtYeif0rGulFQO1ou8JfStrEezh+H8nk/T/9sar+8GC7wufWW3OmdUQkMGMMCHEcYgeFviWLEvax/B8nk/X/2Q==\n", "text/html": [ "\n", " \n", " " ], "text/plain": [ "" ] }, "execution_count": 1, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from IPython.display import YouTubeVideo\n", "YouTubeVideo(\"xEMWh6EjJoY\",560,315,rel=0)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "![hungary_sludge_grap_930321a.jpg](https://raw.github.com/jckantor/CBE20255/master/notebooks/figures/hungary_sludge_grap_930321a.jpg)\n", "\n", "'Red mud' (or 'red sludge') is the by-product of processing of [bauxite ore](http://en.wikipedia.org/wiki/Bauxite) to produce alumina. Alumina has a wide variety of commercial uses including plastic fillers, catalysts, textile dyes, sunscreens, ceramics, and as the raw material for the production of aluminum metal. The [Bayer process](http://en.wikipedia.org/wiki/Bayer_process) is the most commonly used process for the refining of bauxite ore.\n", "\n", "In the Bayer process, the aluminum oxide (Al2O3) present in the raw bauxite ore is digested in a hot solution of caustic soda (NaOH) at 175 °C. The caustic solution reacts selectively with the aluminum oxide to produce aluminum hydroxide that dissolves in solution. The other components in the bauxite ore do not dissolve and are mechanically separated from the liquid stream. The stream of solid waste and entrained processing liquid forms the 'red mud' that was stored behind the earthen dam at the Hungarian facility. \n", "\n", "News reports suggest the red mud produced at this facility was stored with no further processing. This would be inconsistent with modern practice in plants of this type, but if true then the processing liquid entrained in the red mud would be highly caustic and constitute a significant health and environmental hazard.\n", " \n", "The purpose of this project is to use your chemical engineering knowledge to gain a deeper understanding of the issues associated with this tragedy. In particular, the news reports provide little hard information about the actual health or environmental hazards posed by the red mud. Some reports refer to the trace metal content, others cite the caustic components with pH ranging from 10 to 13. \n", " \n", "We approach this problem starting with the basic chemistry involved in the processing bauxite to produced refined alumina. From this we develop a simple process flowsheet and determine the basic operating constraints on the process. This provides a basis for a worst case estimate of the caustic component of the red mud and assessing possible remediation strategies." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Exercises" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Generation/Consumption Analysis\n", "\n", "The portion of the Bayer process we're considering consists of two reactions. The digestion reaction results in the dissolution of aluminum oxide Al2O3 under caustic conditions \n", "\n", "Al2O3 (s) + 2 OH- (aq) + 3 H2O (l) → 2 Al(OH)-4 (aq)\n", "\n", "Following separation of the undissolved solids, the remaining liquor is cooled under controlled conditions in a crystallizer to precipitate aluminum hydroxide Al(OH)3. \n", "\n", "Al(OH)-4 (aq) → Al(OH)3 (s) + OH- (aq)\n", "\n", "The Al(OH)3 precipitate is subsequently heated in kiln to produce the purified alumina product of the process. We omit those downstream operations from our analysis because they do not contribute to the production of red mud. So considering only the upstream operations, prepare a generation/consumption analysis to determine the overall process stoichiometry for the production of the Al(OH)3 precipitate." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Estimating Caustic Soda\n", "\n", "Typical bauxite contains about 40-60 wt% Al2O3 with the remainder comprised of iron oxide, silica, and other minerals. Initial news reports on the Hungarian situation provided information from which we can estimate the approximate composition of the processed bauxite ore as shown in this table.\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
Component wt %
Al2O355%
Fe2O335%
SiO210%
\n", "\n", "Following digestion and separation, assume the solid components of the red mud have a density of 4.0 g/cc, and that 35 vol% of the red mud is comprised of liquid retained in the void spaces. Based on this information, sketch an input/output diagram for the production of Al(OH)3(s) and estimate the minimum amount of water required to process 1000 kg/hr of bauxite ore. \n", "\n", "From this information, can you estimate the amount of caustic soda (NaOH) required to process the bauxite ore?" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Free Energy of Reaction\n", "\n", "The dissolution process involves reaction with hydroxide ions (OH-). Hydroxide ions are formed by the dissociation of water\n", "\n", "H2O (l) → H+ (aq) + OH- (aq)\n", "\n", "The dissociation constant at 25 °C is\n", "\n", "$K_{H_{2}O} = \\frac{a_{H^{+}}a_{OH^{-}}}{a_{H_{2}O}}=1.2\\times10^{-14}$\n", "\n", "The molar enthalpy and molar Gibbs free energy of formation of H+(aq) is normally taken to be zero at STP (1 atm, 25 °C). Using this information and additional data from Table B.3 of your textbook, estimate ΔĜ°f,OH- for OH- (aq)." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Estimating Equilibrium Constants\n", "\n", "The standard free energies of formation for selected compounds of aluminum at 25 °C. are given in the following table. Data is from [Parks, George, American Mineralogist, Vol. 57, pp. 1163-1189 (1972)](http://www.minsocam.org.proxy.library.nd.edu/ammin/AM57/AM57_1163.pdf). \n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
Component ΔĜ°f kcal/gmol
αAl2O3(s) -378.2
Al(OH)-4(aq) -311.0
αAl(OH)3 -275.3
\n", "\n", "In the Bayer process, alumina is dissolved by digesting Bauxite ore with sodium hydroxide (NaOH) solution at 175 °C. The digestion selectively dissolves the alumina by forming Al(OH)-4(aq) according to the reaction\n", " \n", "R1: Al2O3 (s) + 2 OH- (aq) + 3 H2O (l) → 2 Al(OH)-4(aq)\n", "\n", "(a) Using the results of problem 1 and the other available data, estimate the equilibrium constant KR1 of this reaction at standard pressure and temperature.\n", "\n", "(b) Because we don't have enough information to compute ΔĤ°rxn, use the Gibbs free energy as a crude estimate for ΔĤ°rxn. Estimate the equilibrium constant KR1 at 175 °C.\n", "\n", "(c) Write an equation for the relationship between equilibrium concentrations of OH- and Al(OH)-4 in the digester." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Estimate Equilibrium Yield\n", "\n", "Following digestion and separation of the undissolved solids, the liquid from the digester is cooled and seeded to precipitate Al(OH)3 by the reaction\n", " \n", "R2: Al(OH)4- (aq) → Al(OH)3 (s) + OH- (aq)\n", " \n", "The Al(OH)3 precipitate is a white, fluffy powder. \n", "\n", "(a) Estimate the standard molar free energy of reaction ΔĜ°R2 and the associated equilibrium constant KR2.\n", "\n", "(b) Write an equation for the equilibrium relationship between the liquid phase concentrations of Al(OH)4- and OH- in the crystallizer." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Flowsheet Analysis\n", "\n", "Prepare a mass balance of the Bayer process flowsheet shown in the accompanying figure assuming a basis of 1000 kg/hr of bauxite ore. \n", "\n", "![Ajka_alumina_plant.pn](https://raw.github.com/jckantor/CBE20255/master/images/Ajka_alumina_plant.png)\n", "\n", "Assume that the dissolution of Al2O3 reaches equilibrium in the digester at an operating temperature of 175 °C, and that the precipitation reaction reaches equilibrium in the crystallizer at an operating temperature of 25 °C. The spent solution from the crystallizer is mixed with makeup NaOH and water before recycling to the digester for reuse. The portion of the liquor from the digester lost through entrainment in the red mud is equal to 35% by volume of the red mud. The density of solids is 4 g/cc and the density of liquids is 1 g/cc. Adjust the flow of caustic soda to maximize the recovery of Al(OH)3(s). \n", "\n", "(a) What is the maximum recovery of Al(OH)3(s)? \n", "\n", "(b) What is the pH of the red mud when recovery of Al(OH)3(s) is maximized?" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### How much acid should be shipped to the accident site?\n", "\n", "There have been several proposals for remediating the environmental impacts of spilled red mud. At the time of the spill, one village poured vinegar into a river contaminated by the red mud. Suppose 10,000 cubic meters of sludge entered a creek. Assuming vinegar is approximately 5 wt% acetic acid, how much vinegar would be needed to completely neutralize the effects of the red mud? Roughly how many railroad tank cars would be required to transport the vinegar?" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Long Term Remediation\n", "\n", "There have been several proposals for a long term strategy to remediate the red mud from alumina facilities through [CO2 sequestration](http://www.epa.gov/climatechange/ccs/). Global alumina production is estimated to be 80,000,000 metric tons per year of Al2O3.\n", "\n", "Use the results of your calculations above to estimate how much CO2 would be required to neutralize the red mud produced as a by-product of alumina production through the reaction\n", "\n", "CO2 + 2 OH- → H2O + CO3-2 \n", "\n", "[Laboratory measurements](http://www.ncbi.nlm.nih.gov/pubmed/20036053) show that about 5.3 g of CO2 can be sequestered per 100 g of red mud. How does this compare to your estimate? Where would the necessary CO2 come from?" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] }, { "cell_type": "markdown", "metadata": {}, "source": [ "\n", "< [West Virginia Chemical Spill](http://nbviewer.jupyter.org/github/jckantor/CBE20255/blob/master/notebooks/B.01-West-Virginia-Chemical-Spill.ipynb) | [Contents](toc.ipynb) |

\"Open" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.3" } }, "nbformat": 4, "nbformat_minor": 2 }