A new productivity model incorporating the whole domain fracture closure shape under non-uniform proppant placement: insights for proppant placement optimization and cost-effective ceramic substitution