We study a model of a lipid bilayer membrane described by two order parameters: the chemical composition described using the Gaussian model and the spatial configuration described with the elastic deformation model of a membrane with a finite thickness or, equivalently, for an adherent membrane. We assume and explain on physical grounds the linear coupling between the two order parameters. Using the exact solution, we calculate the correlation functions and order parameter profiles. We also study the domains that form around inclusions on the membrane. We propose and compare six distinct ways to quantify the size of such domains. Despite its simplicity, the model has many interesting features like the Fisher-Widom line and two distinct critical regions.